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This Journal is a publication of the Ameri. 
can Rocket Society devoted to the advange. 
ment of the field of jet propulsion through 
the dissemination of original papers disclosing 
new knowledge or new developments. As 
used herein, the term “jet propulsion” 
embraces all engines that develop tirust by 
rearward discharge of a jet through a nozzle 
or duct; and thus it includes air-c: nsuming 
engines and underwater systems as well as 
rockets. JET PROPULSION is open to cop- 
tributions dealing not only with p:opulsion 
but with other — of jet-propel! d flight, 
such as flight mechanics, guidance, t: \emetep- 
ing, and research instrumentat sn, [n- 
creasing emphasis will be given to the cientifie 
problems of extraterrestrial flight. 
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such financial sponsorship appears as 8 
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amuel Herrick is Professor of Astronomy at the Uni- 
versity of California, Los Angeles. There he instituted 
courses and research in celestial mechanics and observa- 
tion theory that are now a part of the astrodynamics pro- 
gram of the University. This program can claim the first 


university course in astronautics. Dr. Herrick’s ‘‘Astro- 
dynamies and Rocket Navigation”’ was instituted in 1946, 
and now the first Ph.D. in the field is his co-author, Dr. 
Baker. Prof. Herrick has held two Guggenheim Fellow- 
ships for his space navigation studies, and is currently con- 
sultant for such corporations as Aeronutronic Systems, 
Rand, Martin-Denver, Librascope, General Precision 
Equipment, General Electric, Systems Corporation of 
America and Republic Aviation. He is chairman of the 
Executive and Space Navigation Committees of the In- 
stitute of Navigation, and is an honorary member of that 
society and its British counterpart. His B. A. is from Wil- 
liams College and his Ph.D. from the University of Cal- 
ifornia, Berkeley. 

Robert M. L. Baker Jr., presently a Staff Member, Aero- 
nutronic Systems, Ine. and an instructor of courses in 
Astrodynamics at the University of California, Los Angeles, 
obtained a B.A. and M.A. in Physies, and a Ph.D. in 
Engineering (and Astronomy). His specialties are in the 
fields of celestial mechanics and rarefield gas-dynamics. 
Dr. Baker is currently engaged in studies concerning space- 
vehicle orbits, meteoritics, low density drag and sputter- 
ing with particular emphasis on high-speed re-entry 
from interplanetary orbits. 


1 Introduction 


PHENOMENON of the post-Sputnik world has been a 

spreading of the realization of a few that celestial me- 
chanics has become an engineering subject. But the classical 
field of celestial mechanics, as delimited by many of its 
present-day practitioners, is at once too broad and too narrow. 
It is too broad in that it extends itself into areas such as the 
structure and motions of the galaxies and even of the stars and 
the nebulae that are found in them. _ It is too narrow in that it 
has not traditionally included forces that are important to the 
study of trajectories of objects such as those that we are now 
sending out to space. Accordingly, the new term “astro- 
dynamics” has come to be used at least partly in its stead. 

Astrodynamics is defined as including those parts of celestial 
mechanics, geophysics, aerodynamics, electromagnetic theory, 
exterior propulsion theory and observation theory that bear 
on the trajectories of planets, comets, meteorites and artificial 
vehicles above the denser portions of the Earth’s atmosphere. 
It may be noted that astrodynamics bears the same relation- 
ship to astronautics as does aerodynamics to aeronautics. 

It is evident that certain of the areas of astrodynamics have 
been overworked while others have been neglected and await 
intensive research. In the former category, we list optimiza- 
tion of orbits, satellite stay time, equatorial bulge perturba- 
tions, ete. The investigation of these topics has actually been 
hindered by the variety of notation introduced and the lack of 
active coordination among the numerous contributors. (It 
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is estimated that there have been published over one hundred 
unrelated and often redundant intercompany memoranda, 
technical notes and papers on the subject of equatorial bulge 
perturbations alone.) On the other hand, articles dealing with 
precision orbits, with precision study of nongravitational 
forces, etc., have been conspicuous by their absence from 
most astronautical literature; and it is in these areas that 
we will look for the future advances in astrodynamics. 

In this connection, one of the most significamt of the recent 
advances in astrodynamics has been the spread of the realiza- 
tion that the subject is a complex one, that not only is there a 
need for a thoroughgoing knowledge of what has been done in 
the whole field of celestial mechanics, but also that the advent 
of artificial vehicles in space is presenting us with new prob- 
lems that are worthy successors to those solved by such men as 
Clairaut, Lagrange, Laplace and Gauss. This realization is 
focused in part on the distinction that we can make between 
“feasibility orbit work’’ and ‘‘precision orbit work.” 

Many excellent feasibility studies have been made in recent 
years. In this country some of the pioneering work may be 
found in the publications of Rand Corporation, and in this 
country and abroad in the publications of the various rocket 
societies, especially the AMERICAN Rocket Society and the 
British Interplanetary Society. A significant amount of work 
has been done also in proposals for study contracts submitted 
to military and scientific research organizations and in inter- 
and intracompany memoranda and technical notes. Feasi- 
bility studies are characterized by the search for optimal orbits 
and minimum requirements upon fuel and guidance. Notable 
progress has been made in error analyses, but they have tended 
to lull people who have only a cursory acquaintance with the 
trajectory field into the comfortable delusion that it is, after all, 
asimple one. This is because these studies may make use of 
approximate theories, whether based on the two-body prob- 
lem or taking into account the perturbations of a third or a 
fourth body or even of the Earth’s bulge. They may be 
based upon approximate coordinates, such as ones based on the 
assumption that the moon is moving uniformly in a perfect 
circle or is fixed ina rotating framework. And finally they may 
be based upon approximate constants, such as the laboratory 
value of G the constant of gravitation, which is good to only 
three significant figures as determined in terrestrial labora- 
tories. 

In precision orbits, on the other hand, one must leave no 
stone unturned in the quest for accuracy in the basic theoreti- 
cal developments of perturbation theory, in the consideration 
of nongravitational forces, in the coordinates and reference 
frameworks, in the constants, in observation and in correction 
theory. The more recent advances in each of these areas will 
be explored in what follows. -—. 


2 Perturbation Theory 

The term ‘“‘perturbation’’ is used loosely. Sometimes it is 
used to designate the sum total of forces acting on an object 
over and above a principal force that is responsible for some 
kind of analytic reference orbit or for any one of the compo- 
nents that goes to make up the sum total. It may be used 
also for the resulting accelerations, or for the departures from 
a fixed reference orbit or for the variations of parameters that 
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define a varying reference orbit. It is evident that these are 
capable of more specific reference as “perturbative forces,” 
“perturbative accelerations,” ‘perturbative departures” and 
“perturbative variations.’”’ By extension the term is applied 
also to a process that is not strictly perturbative, but involves 
numerical integrations of components of the total acceleration, 
without the use of any reference orbit. The latter process, 
which has been widely employed in recent years, sometimes 
without realization of its relationship to the historical develop- 
ment of perturbation theory, is known as Cowell’s method. 
The term ‘“‘perturbation’’ also appears in the terms “special 
perturbations” and “general perturbations,” the former refer- 
ring to the handling of perturbation problems by numerical 
integration and the latter to the development into series, as of 
sines and cosines, that can be integrated term by term. 

In special perturbations the integration of the departures 
from a fixed reference orbit is known as Encke’s method. The 
method of variation of parameters, which has been simplified 
by the work of Herrick (16)? and Musen (28), utilizes a ref- 
erence orbit that varies in such a way as always to agree in 
position and velocity with the actual path. These two meth- 
ods and Cowell’s method indicate the three basic approaches 
to the perturbation problem, but it is possible to find many 
variations or combinations of them, a notable example being 
Hansen’s method. Some criteria for choice between methods 
of special perturbations are well known; others are emerging; 
still others are needed. The best established criteria are for 
choice between Cowell’s method and the other two. 

Cowell’s method is to be preferred when the perturbations 
are very large. Such, for example, are the perturbations of 
the sun on a distant satellite, such as the ninth and twelfth 
satellites of Jupiter; such also might be the case for large 
thrust perturbations during the modification of an inter- 
planetary orbit. In these circumstances the ‘perturbative 
accelerations” would be as large or as unmanageable as the 
principal acceleration term, on which the reference orbit is 
based. Cowell’s method would accordingly be handled with 
as short an interval of integration as the more strictly pertur- 
bative methods, and its equations, not involving the relatively 
complex developments of the analytic reference orbit, would 
be much simpler to compute and program. 

When the perturbation terms are very small compared with 
the principal acceleration, however, Encke’s method and the 
method of variation of parameters have a distinct edge over 
Cowell’s method in that longer intervals may be used in the in- 
tegration, thus reducing both the amount of calculation and 
the inevitable accumulation of end-figure errors in the integral. 
The orbit of the minor planet Icarus, which comes closest of all 
planets to the sun and so is subject to very large primary ac- 
celerations, is a very good example of an integration in which 
Cowell’s method cannot or should not be used (17). Even 
more questionable is the use of Cowell’s method in the integra- 
tion of the trajectories of lunar vehicles. For these trajec- 
tories the perturbations of the sun and moon are almost negli- 
gible at the start, and the many steps, at short intervals of 
time, that are necessary to integrate a Cowell trajectory, be- 
cause of the rapid variation in the Earth’s gravitational attrac- 
tion, may be avoided by the use of a reference orbit based 
simply on the two-body problem. A very large number of 
Cowell integrations of lunar trajectories have nevertheless been 
ground out of modern computing equipment. Such calcula- 
tions may be justified in feasibility work when the amount of 
calculation is of no importance, but when precision requires 
that the accumulation of error in numerical integration be kept 
to a minimum, a Cowell integration scheme should be called 
into question. In one such calculation it is estimated that the 
accumulation of error amounts to nine significant figures! 
It constitutes an advance in astrodynamics then, that several 
organizations are experimenting with Encke integrations in 
the realm of lunar trajectories. 


5 Numbers in parentheses indicate References at end of paper. 


There are several devices by which one may reduce the effect 
of the perturbations. Doing so may make it possible to ignore 
the remaining perturbations entirely in feasibility or other 
approximate studies. One of these is the adjustment of the 
central mass, or of the constant of gravitation, in such a way 
as to transfer an average effect of the perturbations into the 
principal or two-body term. Only the departures from the 
average effects then remain as perturbations. The «verage 
effect of the equatorial bulge, for example, is to augm«nt the 
central mass by an amount depending upon the inclina' ion of 
the orbit. In the n-body problem, the average effect of ob- 
jects inside the orbit under study is to increase the : entral 
mass by an amount somewhat in excess of their own 1:asses, 
Exterior objects decrease the average pull of the «entral 
mass by amounts considerably less than their own masses, 

Another possible way to reduce the perturbations is to ad- 
just or alter the analytical reference orbit. Acting «1 this 
suggestion Baker (2) has shown that a rectilinear gravi: y-free 
drag orbit is preferable as a reference orbit in certain rc-entry 
problems. Sterne (34) and Garfinkel (12) have stud ed in- 
teresting solutions to the satellite problem, when there are 
bulge perturbations present, by employing an exactly -oluble 
Hamiltonian. Unpublished work indicates that there are 
special or approximate solutions of the three-body prol)iem of 
the same character. For feasibility work, as indicated ::bove, 
these solutions may provide us with useful information ‘hat is 
relatively complete. For precision work, however, it must be 
remembered that at best they are approximations to tlie real 
circumstances. They may supply a useful reference orbit for 
Encke-type perturbations or for the variation of parameters, 
but additional perturbation developments, as carried out in 
(2), must be undertaken for the full solution of the problem. 
Accordingly the reduction in the complexity of the perturba- 
tion work must be balanced against the complexity of the ref- 
erence orbit. If the reference orbit is exceedingly complex, in- 
volving functions that are difficult to compute or tabulate, the 
new reference orbit may involve more complications than it 
eliminates in the perturbations. The three-body problem has 
been disappointing in this respect; see (23 and 26). Ana- 
lytical solutions are very lengthy and would probably involve 
impractical calculation and an inordinate number of terms. 
The simple integral named for Jacobi has been used effectively 
in checking perturbation integrations and may in the future 
serve a useful purpose in the variation of parameters. 

The satellite problem has proved to be immeasurably more 
complex than it was generally conceded to be before Sputnik. 
Numerical integration is unsatisfactory because of the accu- 
mulation oferror. Series integration is unsatisfactory because 
of the excessive number of terms that must be carried. Fur- 
ther experimentation is to be highly desired. It should in- 
clude such matters as the investigation of Tchebycheff and 
other polynomials in series integration, modification of the 
equations of motion, modifications of the reference orbit, and 
combinations of series integration with numerical integration. 

Ion or other low-thrust propulsion offers an interesting chal- 
lenge to perturbation theory. The continuous thrust is easy 
enough to handle theoretically, see for example (35), but it 
may be considerably smaller than other perturbative forces, 
and so may become lost in the uncertainties of the latter. 
Thus for a take-off from a circular satellite orbit the design 
trajectory may be a uniform spiral, whereas the actual! path 
will probably develop a sort of “elliptic resonance,” developing 
into a series of intertwined loops from which the “break-away 
direction will be highly uncertain. Paradoxically, this ci 
cumstance will be less troublesome for interplanetary voyaging 
than for trips to the moon. For the latter it is desirable, to 
put it mildly, to go somewhat in the direction of the moot 

rather than a contrary direction. For the former, when the 
vehicle starts to recede from the Earth it will begin to slow dow? 
so markedly that, whatever the direction, it will be tray eling 
on a heliocentric orbit with approximately the same velocity 
and very nearly the same shape as the orbit of the Earth. Ac 
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erdingly, subsequent low thrust will follow a pattern very 
guch the same, whatever the direction of the vehicle from the 
farth. Once free from the Earth, however, a heliocentric 
diiptical resonance may prove troublesome. 


3 Nongravitational Forces 


y recently have investigators been stimulated to look 
dosely at nongravitational forces in space-vehicle motion. 
(Unfortunately the investigations stem more from attempts to 
aplai away erroneous results of purely gravitational orbit 
determ:ination than from basic interest in the subject for its 
gmsake. The principal nongravitational forces that appear 
ignifi:antly to affect an object in space are aerodynamic drag, 
dectrial and electromagnetic forces, meteoritic drag, radia- 
tion pressure and the relativistic effects. (The last effects, 
ithouzh strictly speaking not nongravitational, are placed in 
the nongravitational category for convenience.) 

Altl ough aerodynamic drag is small for an orbit above the 
denser atmospheric levels of the Earth and planets, it must be 
treate'| with considerable care as a perturbation. Recent in- 
vestig:: tions of the drag coefficient at high altitudes (4, 27 and 
5) have all shown that the drag coefficient is significantly less 
than the so often employed value of 2. In fact, the transi- 
tional parameter (the reciprocal of a modified Knudsen 
number, a or B) has the value of approximately 0.4 fora satel- 
lite at a perigee altitude of 200 km (and with a 1 m radius, 
surface temperature of 300 K and an orbital semimajor axis of 
?Earth radii) and, consequently, reduces the drag coefficient 
by approximately 5 to 10 per cent. This recent advance in the 
serodynamie side of astrodynamics should soon find employ- 
ment in satellite orbit determination. 

Once in orbit, the space vehicle will acquire an electrical 
tharge that may interact with the Earth’s or sun’s magnetic 
field and perhaps with a highly ionized atmosphere, to produce 
woth dissipative and nondissipative forces. The factors which 
ifluence the satellite’s charge are the photoelectric effect, ion 
jickup and ion conduction. The resulting charge appears to 
vequite small, but the forces (e.g., Lorentz, induced drag, etc.) 
rsulting from this charge could become important as pertur- 
ations. Pioneering investigation of these electrical forces 
ave been carried out by Howard Chang and Myron Smith of 


the Rand Corporation and by the authors listed in (14 and 21). 
Inaddition to such forces, one might conceive of a magneto- 
tydrodynamie “shock wave” that would affect the motion of a 
space vehicle. More investigation of these areas is certainly 
indicated. 

Meteoritie drag, occasioned by the momentum transfer of 
uiero-meteorites impacting on the surface of a space vehicle, 
vould be quite similar in action to aerodynamic drag. The 
juestion of the density of meteoritic material in space must be 
wsolved, therefore, before adequate conclusions can be drawn 
irom satellite observations. A large amount of theoretical 
work in this regard can be found in the numerous papers of 
Whipple and of Kallmann. Typical of the analytical work 
involved in the interpretation of micro-meteorite detectors on 
board a space vehicle is the less well known work of Bauer (5). 
In this connection the problem of surface pitting, not only oc- 
asioned by micro-meteorites but also by impinging molecules 
sputtering) needs careful attention. The phenomenon of 
molecular sputtering as studied by Wehner, Henschke, Honig 
ind others has interesting effects, not only on surface deterio- 
ration, but also on drag and heat transfer, as noted in (2). 
Radiation pressure will be particularly influential in the 
inalysis of the orbital motion of low mass and large area 
bodies, such as dust particles, balloons, etc. The analysis of 
weh bodies ranges from the classical work of Poynting and 
of Robertson (32) (done in connection with the study of inter- 
jlanetary dust) to the more recent investigations of Garwin 
13). 

Relativistic effects are not only of importance in the motion 
of the perihelion of Mercury (11, 37), but also in orbital 
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studies of artificial asteroids (38) and satellites of Earth 
(25). The problem of the efficient incorporation of Einstein 
terms into perturbation theory is of great importance and 
deserves careful future investigation. 


4 Differential Relationships 


Differential formulae are useful in the improvement of an 
approximation to an orbit into a better approximation, in the 
correction by thrust of a physical orbit that does not intercept 
an objective into one that does, and in error analyses or studies 
of the propagation of uncertainties in initial conditions to un- 
certainties in terminal conditions. In astronomical work, 
heretofore, it has been customary to obtain the necessary for- 
mulae by differentiation of the basie two-body formulae, 
usually without any attempt to determine the effects of these 
changes on the perturbations. Modern electronic calculating 
equipment, however, has given us a powerful new tool for such 
analyses. Instead of differentiating the basic formulae, we 
are now inclined to repeat the basic calculations once for each 
of the constants that is to be altered in the correction or whose 
uncertainties are to be studied. From the differences in the 
end results that correspond to arbitrary alterations in the 
basic constants, made one at a time, the differential coefficients 
are determined. These can as well as not include the effect 
on the perturbations. 

The differential correction of orbits is not, as some have sup- 
posed, a matter of determining from observation at an instant 
all of the components of position and velocity. It involves in- 
stead the calculation of what an observation would be if the 
adopted elements were correct. These elements may be the 
six customary elements of a two-body orbit, or the six com- 
ponents of an initial position and velocity, or some other ap- 
propriate set of constants. The discrepancies between calcu- 
lation and observation are called residuals; it is these that are 
related differentially to the desired corrections of the starting 
elements. It follows, then, that if there are six constants to be 
corrected, a minimum of six residuals is needed. These re- 
siduals may be of any accurately observed quantities or of any 
combination of observed quantities, subject of course to special 
circumstances that may lead to indeterminateness in the 
solution. 

Rapidly moving Sputniks offer a special challenge to dif- 
ferential correction theory, especially in view of the difficulty 
of getting an adequate perturbation theory. For finding 
ephemerides, at least, it may prove to be preferable to adopt 
a program of continuous correction, based on observation (1), 
instead of developing a thoroughgoing perturbation theory. 
Without perturbations or with an insufficient perturbation 
theory, however, there is a limit to the usefulness of the earlier 
observations as compared with the later observations. Since 
the early observations have an important bearing on the accu- 
rate determination of the period and semimajor axes, how- 
ever, they cannot be completely discarded. It is necessary to 
develop a system in which their effect is properly weighted. 
An interesting essay in this direction, of a preliminary and 
theoretical rather than practical character, was presented by 
Swerling of the Rand Corporation to the Los Angeles “‘Astro- 
dynamics Colloquium” in May 1958. 


5 The Determination of Orbits 


In classical astronomical theory it has been customary to 
determine circular orbits from two observations of direction 
only or from one observation of the two direction angles and 
their velocities, and to determine elliptical orbits from three 
observations of the two direction angles. The circular orbit 
procedures appear to be directly applicable to the determina- 
tion of Sputnik orbits, but the elliptical orbit determination 
procedures are unlikely to offer much to the problem, because 
the slight departure from great-circle motion usually en- 
countered requires accuracy of observation beyond that which 
seems likely in satellite observations. 
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The loss of the classical methods for the determination of 
elliptical orbits, however, seems to be made up for by many 
new techniques such as those employed in (20 and 36) as well 
as those discussed in the references compiled by Benton (6). 
Most of these are approximate in character, ranging from 
rough estimates as to when the satellite first passes over an 
expected point to determinations from radio Doppler obser- 
vations. The state of orbit determination, however, appears 
to be somewhat chaotic, and it is questionable whether it will 
function well for an unexpected and unannounced satellite. 


6 Requirements on Observation 


The comparative precision of the various types of observa- 
tion used with satellite and space-vehicle orbit determination 
is understood somewhat better now than it was a few months 
ago. The accuracy to be sought in visual observations should 
preferably be not less than 2 sec (one part in 105) for close, fast 
moving satellites, and should range down to one part in 10° 
or even better for an object at 5 Earth radii. The 2 sec 
accuracy has been claimed for the ballistic cameras of the 
Cape Canaveral chain, see (8), and has probably been 
achieved also with wide angle photographic telescopes at astro- 
nomical observatories. The stituation with the yet untested 
and long delayed Baker-Nunn cameras is not clear; it may be 
that 1 min is a more realistic estimate for these and perhaps 
even for the other cameras mentioned. 

Two seconds of are (10~5 radians) nevertheless stands as a 
challenge to the electronic methods of observation. Observa+ 
tions of lesser accuracy will be unlikely to have much value 
in the improvement of our knowledge of the gravitational en- 
vironment. For the present, accordingly, we shall tend to 
discount radio Doppler observations (22) and to hope that 
transponder-type observations of high accuracy will be 
planned increasingly. 

The question of observation is often associated with the 
availability of a stable satellite platform in space. In this 
area of astrodynamics the work of Roberson, Klemperer and 
others (30, 31, 24, and 3) has been particularly interesting and 
useful, not only with application to satellite stability, but also 
with respect to lunar librations. 

It is important to distinguish between the processes of ob- 
servation reduction and orbit determination. Specifically, an 
observer or observation group should not consider it necessary 
or even desirable to combine incomplete observations to de- 
termine a position. By incomplete observations we mean 
measures of one or two components of direction without dis- 
tance, or of distance without direction. With such observa- 
tions, observers unacquainted with astronomical techniques 
tend to believe that it is necessary to triangulate in order to 
obtain a position. Such a position, however, may multiply 
the uncertainties of the individual observations to the point 
of vitiating their value. It is far better to supply the orbit 
analyst with the individual observations, not necessarily in 
groups capable of triangulation. He will then be able to 
“represent”’ the observations, i.e., calculate what they would 
be if this assumed orbit were correct. The residuals or dis- 
crepancies between the observed and the computed values 
then supply him with a basis for the correction of his assumed 
orbit (compare section 4). 

Accordingly, we may see that any observed datum, a dis- 
tance (or‘ “slant range’’), a direction angle, a radial velocity 
(or “slant'range rate’’), etc., can be used by the orbit analyst. 

A corollary to the foregoing discussion is that the observer 
need not attempt to find a zero point, as on a Doppler curve, 
from the gross values of points on the curve. The orbit 
analyst instead can compute residuals, and either determine 
the zero point to their lesser range, or leave it as an unknown 
for his final solution. 

The proper province for observation reduction then is the 
correction of raw data for known systematic errors, instru- 
mental ones, such as level and azimuth and collimation, or 
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environmental ones, such as refraction and aberration. Even 
geocentric parallax, especially since it is so large, should he left 
to the orbit analyst. 


7 Coordinate Systems 


The observations themselves do not have to be reduce: to a 
common reference system, but the systems in which thcy are 
made must be closely tied to one. Observed coordinat:s are 
necessarily “topocentric,” i.e., referred to the place (“topo~”) of 
observation; their frame of reference may be tied to the local 
horizon or to the equator. In either circumstance th» co- 
ordinates of the observer (latitude, longitude, height abo. e sea 
level) and the time must be accurately known and reduce: toa 
specified reference spheroid. The connection betwee the 
observation frame of reference and the geocentric dyna:nical 
frame of reference is then handled partly in the reduction of 
observational data (as in topocentric coordinates or dir tion 
cosines) and partly in the orbit calculation. 

The dynamical frame of reference is usually not an invrtial 
one, especially for geocentric orbits, since relative motio:, can 
be handled as easily as inertial motion. But the refe:ence 
axes should be unaccelerated; the logical choice for such: axes 
is equatorial and equinoctial, i.e., with the plane of the eqiiator 
adopted for the zy-plane, and the z-axis directed to the v-rnal 
equinox. Even in these circumstances, because of prece=sion 
and nutation, we must choose between several equators and 
equinoxes: 

1 Fixed to such a timeless frame as the mean equator and 
equinox of 1950.0. 

2 Fixed to the mean equator and equinox of the beginning 
of the year. 

3 Fixed to the mean equator and equinox of some ar- 
bitrary “epoch,” perhaps a date during or shortly before 
flight. 

4 Moving with the mean equator and equinox of date. 

5 Moving with the true equator and equinox of date. 

6 Fixed to the true equator and equinox of some arbitrary 

epoch. 
The first three, all inertial, differ by a set amount of the 
“precession of the equinoxes”’; the fourth and fifth are non- 
inertial, being affected by progressive precession 3 to 4, and 
nutation, 4 to 5. System 6 is also inertial, differing from sys- 
tem 3 by a fixed amount of nutation. The orbital work, then, 
would have to be done in such a frame as 1, 2, 3 or 6; but 
unfortunately the observer and possibly his reference axes are 
necessarily in system 5. For a lunar vehicle restricted to a 
few days in orbit, it may be preferable to work in system 6, 
neglecting its departures from system 5, but for a long term 
satellite it will be necessary to consider the problems of pre- 
cession and nutation with great care. 

When observations are made with altazimuth instruments, 
such as ballistic cameras, and the angular coordinates are read 
from graduated circles, the full effects of refraction, precession 
and nutation must be taken into account. When they are 
made against a stellar background, by measurement of a pho- 
tographic plate, some of these effects are absorbed by the star- 
place reductions, but at least the geocentric coordinates of the 
observer must be corrected for precession and nutation. 


8 Astronomical and Gravitational Constants 


The problems associated with accurate constants may be 
illustrated in lunar-flight trajectories by the fact that the best 
values of the geocentric gravitational constant and the lunar 
parallax as found by Herrick, Baker and Hilton (18) lead to 
uncertainties of the order of 100 miles at lunar impact. It is 
clear that the constants will have to be improved for future 
spot landings designed to supply depots or expeditions. It is 
also clear that precise observations of early flights should con- 
tribute to the improvement of the basic constants. 

For flights to Mars and Venus even greater uncertain(ies 
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grise because of the inexactness of the solar parallax, which is 
essentially the ratio of the Earth’s equatorial radius to the dis- 
tance of the sun (the astronomical unit). It is best to lo- 
calize this effect in the conversion of the burnout velocity from 
ft per sec to astronomical units per day; it is easy then to show 
that the uncertainty in a purely ballistic trajectory at Mars 
or Venus would be of the order of several hundred thousand 
miles. Having made the shift to the heliocentric system of 
ynits and constants, however, it is possible to adjust the ballis- 
tie trajectory and cancel out much of the error in the constants 
by a program of observation and corrective thrust. This cor- 
rection can be carried out, because the observations are al- 
ready, ‘1 effect, in the heliocentric system. 

The asic data on which the gravitational constant is based 
were acded to critically and effectively by the recent achieve- 
ment 0: the Army Map Service enclosing two latitude ares in 
this hemisphere and in the other (9). From these measures 
is determined an excellent value of the Earth’s equatorial 
radius. Joined with the equatorial radius in the determina- 
tion of the constant of gravity is the acceleration of gravity at 
theequator. The accuracy of the latter requires both a care- 
ful stuc'y of the statistical determinations of the acceleration 
of gravity at the equator from relative observations made 
around the Earth, and of the basic determination of the acceler- 
ation oi gravity in Potsdam in 1906. The latter determination 
has been subject to correction by later determinations at 
Washington (19) and at Teddington (10). Other accurate 
determinations are known to be in process. 

The distance of the moon has not been determined by 
parallax measures since the series of observations made in the 
early part of this century at Greenwich and the Cape of Good 
Hope; a more recent determination is reported by O’Keefe 
and Anderson. These observational determinations are some- 
what inferior to the determination dynamically, that is, from 
the constant of gravitation and the period of the moon. 

The critical determination of the solar parallax has been 
subject to a number of painstaking investigations in recent 
decades. The best value to date seems to be that of Rabe (29). 

Recent years have seen a revolution in the accurate deter- 
mination of the fundamental time variable. The irregu- 
larities in the rotation of the Earth have been tied down with 
increasing accuracy (7), and the result has been the establish- 
ment of a new kind of time, more nearly Newtonian, which is 
to be called Ephemeris Time. It will be the basis of the 
tables of the sun, moon and planets to be published in the 
National Ephemerides beginning with 1960. The correction 
of Greenwich time to Ephemeris time is well presented in 
Sadler (33). 
9 The Future 

The current great upsurge of interest in astrodynamics is 
attested by the many groups being founded in the missile in- 
dustry to study trajectories. The number of rockets sent 
theoretically to the moon is beginning to rival the number of 
persons killed in mystery literature. In some instances the 
calculation of a single trajectory to the moon fills the heart of 
its calculator with perhaps too much of a sense of achievement 
and knowledge. But in the great majority of instances the 
young men who are assigned to this kind of work soon dis- 
cover the vastness of the field and tasks before them. They 
are the men who will make great contributions in the field of 
celestial mechanics in the future. Some of them have already 
begun to do so. 

The universities are beginning to rise to the challenge pre- 
sented to them by the increase of interest in this field. We 
venture to describe the U.C.L.A. program, in part because of 
our familiarity with it, and in greater part because it can claim 
priority in the field. Formally, the program may be said to 
have started with the introduction of the course “Astrody- 
namics and Rocket Navigation” in 1946, although older 
courses in the determination of orbits and the reduction of 
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observations are now an integral part of the program. To 
these have been added four graduate courses. The first 
master’s degrees have been granted under the program this 
year and, in conjunction with the Department of Engineering, 
the doctorate. Its program of research dates back to 1937. 
Through the Institute of Navigation, it sponsors the Astro- 
dynamics Colloquium, which every three weeks brings to- 


gether about fifty persons from local research institutions and 


the missile industry for discussions on matters concerned with 
precision trajectory work. Cooperation with industry is re- 


garded as an essential part of the program, and it is recognized 


that the future of celestial mechanics is in the hands of a much 
broader group than was ever the case previously. In return, 
industry has come to recognize astrodynamics as an essential 


specialization, no longer to be touched upon lightly without 
appreciation of its basic complexity. 
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The paper reports the results of some experiments con- 
ducted for investigating the influence of the length of the 
combustion chamber, the steady-state chamber pressure, 
the shape of the nozzle, and the equivalence ratio on high- 
frequency combustion pressure oscillations. 


mixed gases, thereby simplifying the experimental work 
and eliminating the effects of such variables as atomiza- 
tion, vaporization, mixing, etc. The premixed gaseous 
propellants used were superheated propane and _ air, 
ethane and air, methane and air, hydrogen and air, and 
ethylene and air. For the longitudinal mode of combus- 
tion pressure oscillation the frequencies ranged from 570 
to 1750 cycles per sec. The results indicate that a relation- 
ship exists between the amplitude and frequency of the 
combustion pressure oscillations, the combustion cham- 
ber geometry, and the burning rate of the propellants. 


Introduction 


N EXPERIMENTAL study of some of the factors affect- 
ing combustion pressure oscillations has been initiated 
at the Jet Propulsion Center, Purdue University, with the 
object of determining the effects of certain pertinent variables 
upon those oscillations. The research is in its initial phases 
and is being conducted with a rocket motor burning premixed 
gases as the propellants. In this manner several of the vari- 
ables which ordinarily enter into the combustion process for 
a liquid propellant rocket motor, such as atomization, phase 
change, mixing, etc., are eliminated. 

The initial phase of the experimental program is concerned 
with the studies of the influence of the following variables 
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Satellite 


when burning a given set of premixed gaseous propellants: 
Length of the combustion chamber for a fixed diameter. 
Shape of the converging portion of the exhaust nozzle. 
Steady-state combustion pressure. 

Equivalence ratio for the propellant combination. 
Diameter of the combustion chamber for a fixed length. 
Location and method for injecting the propellants. 

The discussions in the present paper are limited to the ex- 
periments concerned with studies of the influence of the fol- 
lowing upon the frequency and amplitude of high-frequency 
combustion pressure oscillations: 

1 The length of the combustion chamber for a fixed di- 
ameter. 

2 The shape of the subsonic part of the nozzle. 

3 Variation in the steady-state combustion pressure. 

4 Variation in the equivalence ratio. 

To date, investigations have been conducted with the fol- 
lowing propellant combinations: Premixed superheated pro- 
pane and air, premixed ethane and air, premixed methane and 
air, premixed hydrogen and air, and premixed ethylene and 
air. 


43 
Experimental Apparatus) 


Fig. 1 presents a cross-sectional drawing of the research 
rocket motor. The internal diameter of the cylindrical com- 
bustion chamber is 33 in. Its length was varied by changing 
the number of uncooled cylindrical sections. The gaseous fuel 
and air were introduced into the mixing chamber through con- 
centric holes; the fuel entered through the central tube and 
the air through the annulus surrounding the latter. The fuel 
and air were thoroughly mixed by causing the combustible 
mixture to flow through a screen and around several baffles 
located in the mixing chamber. The injection system was 
designed so that it operated with a pressure drop of «p- 
proximately 150 psi thereby suppressing any tendency for 
low-frequency combustion pressure oscillations to occur. 

The pressures in the combustion chamber were measured 
with water-cooled capacitance type pressure transducers 
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Photocon Model 342) and their associated electrical equip- 
ment (Dynagage Photocon Model DG-101 and Power Supply 
P§102). The electrical outputs of the transducers were re- 
corded on a six-channel, cathode-ray oscillograph (Hathaway 
Model SC-16B). 

The Photocon transducers were calibrated statically before 
each experiment. Photocon transducers which were to be 
operated with water-cooling during an experiment were also 
water-cooled during their calibrations. The calibration pres- 
sures and the combustion pressures measured during an ex- 
periment were recorded on the same oscillograph record. The 
dynamic characteristics of each of the Photocon transducers 
employed in the experiments and also their associated elec- 
trical equipment were determined by the shock tube tech- 
nique described in (1 and 2).4 


Experimental Results 


The results are presented under the subheadings: Classifica- 
tion of the Oscillations, Effect of the Chamber Length, Effect 
of the Nozzle Shape, and Effect of the Steady-State Combus- 
tion Pressure and Equivalence Ratio. 


Classification of the Oscillations 


The combustion pressure oscillations of a longitudinal mode 
have been classified into two principal categories called for 
convenience the “shock” type and the “sinusoidal” type. 
Other intermediate categories exist. The shock and sinu- 
soidal types correspond to the oscillations having large and 
smal! amplitudes respectively. 

Fig. 2 is an oscillograph record illustrating the shock type 
of oscillation, so termed because the extreme rapid rise in 
pressure at the beginning of each cycle is quite similar to that 
fora normal shock wave (1). 

Fix. 3 illustrates one cycle of the shock type for a gaseous 
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bipropellant rocket motor burning superheated propane and 
air. It can be seen that the overall pressure rise comprised 
two rapid increases in pressure. The first occurred when a 
longitudinal pressure wave moving toward the injector passed 
the pressure sensing diaphragm of the Photocon transducer. 
The second pressure rise was due to the same pressure wave 
passing the pressure transducer after being reflected from the 
face of the injector. The extremely short time in which the 
pressure increases occurred indicates that the pressure wave 
was substantially a shock wave. 
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Fig. 4 is an oscillograph record illustrating the sinusoidal 
type of high-frequency combustion pressure oscillations ob- 
tained with the gaseous bipropellant rocket motor. 

Fig. 5 illustrates a transformation from steady burning to 
the shock type of oscillation when burning premixed super- 
heated propane and air. Similar transformations were ob- 
served for the other propellant combinations investigated. 
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Fig.4 Oscillograph record illustrating oscillations of a sinusoidal 
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Effect of the Chamber Length 


Table 1 presents typical data obtained from experiments 
with superheated propane and air as the propellants. The 
data illustrate the effect of changing the combustion cliamber 
length, the only variable in the series of experiments. 


Table 1 Typical reduced data from the study on the | 
effect of chamber length for superheated propane anc air | 


Chamber Type of Frequency Ampli' ude 
length oscillation of the of the 
in. oscillations oscillations 
eps psi 
28 shock (L)* 570 22 
22.3 shock (L) 685 22 
16.3 shock (L) 948 17 
12.6 sinusoidal (L) 1219 4 


Values for the variables: 
Equivalence ratio, 1.2 + 0.08 
Propane line pressure drop, 158 + 8 psi 
Air line pressure drop, 150 + 10 psi 
Mean combustion pressure, 45 psia + 2 psi 


* L indicates a longitudinal mode of oscillation. 


The data obtained from experiments with other pro)ellant 
combinations indicated that varying the length of the com- 
busion chamber influences the amplitude and frequency of the 
combustion pressure oscillations in a similar manner. 

Gaseous bipropellant rocket motors burning premixed 
superheated propane and air and having combustion cham- 
ber lengths shorter than 7 in. did not exhibit high-fre- 
quency combustion pressure oscillations of a longitudinal 
mode. For those premixed gaseous propellants, a length of 7 
in. may be considered to be the critical length below which no 
high-frequency combustion pressure oscillations of a longi- 
tudinal mode occur. Table 2 presents the measured critical 
lengths for four of the propellant combinations that were in- 
vestigated. 


Table 2 Tabulated data of critical lengths* 
Oxidizer Fuel Critical length 
in. 
air superheated propane C;Hs 7 
air ethane C2H, 8 
air methane CH, el 11 
air ethylene C,H, 6 


* The parameters held constant were: 
combustion pressure (45 psia), 
alence ratio. 


(1) steady-state 
(2) nozzle shape, (3) equiv- 


Effect of the Nozzle Shape 


Fig. 6 illustrates the shapes of the three nozzles employed 
for determining the influence of the shape of the entrance 
section (subsonic portion) of the nozzle upon the amplitude 
and frequency of the combustion pressure oscillations. Based 
upon their axial length the three nozzles will be designated as 
the 1 in., 2 in., and 4 in. nozzle, respectively. 

Figs. 7 and 8 present typical portions of the oscillograph 
records of the combustion pressure oscillations obtained i 
experiments in which the only variable was the length of the 
subsoni@ portion of the nozzle. The propellant combination 
was premixed methane and air. 

Fig. 7 is an oscillograph record obtained with the 1 in. noz- 


ale. “Channel 2 presents the static pressure as measured by 8 
Bs ‘ Photocon transducer located 2 in. upstream from the entrance 
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tothe nozzle. Channels 1, 3 and 4 present the static pressures 
in the plane 2 in. downstream from the injector face. All of 
the statie pressure records (channels 1, 2, 3 and 4) indicate 
that the combustion pressure oscillations were of the shock 
type. 

Fig. 8 is an oscillograph record obtained from an experiment 
with the 4 in. nozzle. The wave shape of the static pressure 
record obtained from channel 2 (the Photocon transducer lo- 
cated near the nozzle entrance) differs markedly from that of 
the shock type. In addition, the cor responding pressure 
records obtained from channels 1, 3 and 4 (the Photocon 
transducers located near the injector face) also have wave 
shapes that are different from the shock type. The difference 
may be attributed to the fact that the entrance section of the 
nozzle no longer affords a surface from which pressure waves 
can be reflected in the longitudinal direction with full in- 
tensity but causes the reflected pressure waves to interfere 
with each other. 


The Effect of the Steady-State Combustion Pressure and 
Equivalence Ratio 


In the investigation concerned with the effect of the steady- 
state combustion pressure the latter was varied from 25 to 130 
psia, approximately. Fig. 7 is a portion of the oscillograph 
record obtained with a steady-state combustion pressure of 
52 psia when burning premixed methane and air at an 
equivalence ratio of 0.88. The record indicates that the com- 
bustion pressure oscillations were of the shock type; their 
frequency is 920 eps, and the peak to peak amplitude is 20 
psi. 

For the same premixed propellants burned at the same 
equivalence ratio but with a steady-state combustion pressure 
of approximately 130 psia the combustion pressure oscillations 
were of the shock type; their frequency was 1000 eps, and 
the peak to peak amplitude increased to 35 psi. 

From additional data similar to those presented in Fig. 7 it 
is concluded that increasing the steady-state combustion 
pressure increases the frequency of the combustion pressure 
oscillations and also the peak to peak amplitude of the oscilla- 
tions. The aforementioned effects were noted for all of the 
gaseous bipropellant combinations investigated. 

During the course of the investigation of the effect of the 
stea(!y-state combustion pressure on the combustion pressure 
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oscillations, no combustion pressure oscillations were ob- 
served when the premixed methane-air propellant combina- 
tion was burned at equivalence ratios above 1.3, regardless of 
the magnitude of the combustion pressure, over the range 25 
to 130 psia. It was decided, therefore, to investigate the 
behavior of the other premixed gaseous propellant combina- 
tions at different equivalence ratios. Although that  in- 
vestigation has not been completed, the results obtained to 
date are of sufficient interest to be presented at this time. 
Fig. 9 presents regions in which high-frequency combustion 
pressure oscillations occurred for the premixed propellants 
methane and air, ethane and air, and ethylene and air. In the 
latter figure the curves do not represent sharp lines of de- 
marcation between stable and unstable burning, but rather 
the loci of midpoints of zones wherein the combustion pressure 
oscillations changed from a high amplitude shock type to a 
low amplitude sinusoidal type, and eventually into smooth 
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combustion free from combustion pressure oscillations. For 
a given steady-state combustion pressure it was found that 
the aforementioned zones could spread over 0.10 units of 
equivalence ratio. 

Another observation made during the investigation of the 
effect of the steady-state combustion pressure on the combus- 
tion pressure oscillations was that increasing the steady-state 
combustion pressure caused the mode of the pressure oscilla- 
tion to shift from a predominately longitudinal mode to a 
transverse mode. Furthermore, the results indicate that the 
transverse modes of oscillation close to the injector are of 
much larger amplitude than those near the entrance to the 
exhaust nozzle. Of the premixed gaseous propellant combina- 
tions investigated, hydrogen and air gave transverse modes of 
pressure oscillation having the largest amplitude. On the 
other hand, the combination of premixed superheated pro- 
pane and air gave the lowest amplitude for the transverse 
mode. Since this phase of the investigation is still in progress, 
the results presented should be regarded as being incomplete. 


Discussion of Results 


The results of the investigations presented in the foregoing 
are discussed under the subheadings: The Effect of the Com- 
bustion Chamber Length, The Effect of the Nozzle Shape, 
and Effect of the Steady-State Combustion Pressure and 
Equivalence Ratio. 


The Effect of the Combustion Chamber Length 


From the investigation of the effect of combustion chamber 


length it is apparent that the amplitude of high-frequency 
combustion pressure oscillations ranges from a low amplitude 
sinusoidal type illustrated in Fig. 4 to the high amplitude 
shock type illustrated in Fig. 3. For intermediate ampli- 
tudes the wave shape is substantially triangular with the 
changes in pressure, increasing and decreasing, occurring in 
approximately the same time interval. 

Examination of the combustion pressure records for the pre- 
mixed superheated propane and air propellants shows that the 
shock type oscillations predominated when the latter occurred 
in motors longer than 16.3 in. (see Table 1). In motors having 
combustion chambers shorter than 12.6 in., the pressure os- 
cillations were predominantly of the sinusoidal type; and 
their amplitude was only 4 psi, measured from peak to peak, 
compared to 22 psi for the shock type occurring in the longer 
motors. Finally, with the shortest combustion chamber (7 in. 
long) no longitudinal combustion pressure oscillations of 
either type were observed. 

The foregoing results indicate that the length of the com- 
bustion chamber has a decisive influence upon the magnitude 
and type of the combustion pressure oscillations of the longi- 
tudinal mode. Moreover, from the observations one can 
postulate a mechanism for explaining the influence of the 
length of the combustion chamber upon the magnitude of the 
combustion pressure oscillations of a longitudinal mode. 

Consider a small pressure disturbance traveling longi- 
tudinally downstream from the injector toward the nozzle. 
During the time that the pressure disturbance is moving in 
the downstream direction, the combustion process in the 
vicinity of the injector proceeds at a normal rate. After the 
pressure disturbance has been reflected upstream from the 
nozzle and just prior to its arrival at the injector face, it must 
traverse recently ignited and unburned propellants. The 
pressure disturbance compresses the unburned and burning 
mixture thereby increasing its temperature and consequently 
its burning rate. The increased burning rate, in turn, ampli- 
fies the pressure disturbance, because additional hot combus- 
tion gases are evolved at a rate faster than that for the steady- 
state combustion pressure. The amplified pressure disturb- 
ance is reflected from the injector face and again passes 
through the unburned mixture further increasing its burning 
rate and thereby further amplifying the pressure disturbance. 
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It may, therefore, be expected that in the case of rocket 
motors having a short tubular combustion chamber, so that the 
time for the pressure disturbance to travel twice the length 
of the combustion chamber is short, there will be only a small 
quantity of combustible mixture injected in the time interval] 
for the reflections of two successive pressure waves frorn the 
face of the injector. Because of the smaller quantities of 
unburned and burning propellants traversed by the pressure 
disturbance the smaller will be the energy release, an con- 
sequently the combustion pressure oscillations, if they \ccur, 
will be of small amplitude. It is also conceivable that if th. 
combustion chamber is short enough the amplitude «f th 
combustion pressure oscillations of a longitudinal mod: may 
be too small to be significant. The length of combustion : ham- 
ber corresponding to the latter condition has been termed th 
critical combustion chamber length. 

A long combustion chamber affords the pressure di-turb- 
ance a relatively long time for traversing twice the lens th of 
the combustion chamber, and likewise increases the an) ounts 
of unburned and burning propellants between successi\e re- 
flections of the pressure wave. Consequently, a large pres- 
sure amplification may be expected because of the ‘arger 
amount of thermochemical energy released by the con:pres- 
sion of the large quantities of unburned and burning propel- 
lants by the pressure disturbance. 

From Table 2 it is apparent that the critical length of thi 


fen chamber varies with the propellant combin:ttion. 
7 


There is reason for assuming that the critical length of the 


~combustion chamber is related to the flame speed of the pro- 


pellant combination (the ethylene-air combination has the 
highest flame speed and methane-air the lowest for the pro- 
pellants investigated) (3). 

Experimenters using liquid bipropellant rocket motors have 
also noted that long combustion chambers tended to initiate 
combustion pressure oscillations of the longitudinal mode 
more readily than do short ones (4, 5). The mechanism 
postulated in the foregoing, while based on experiments with 
rocket motors burning premixed gases, also explains the 
greater tendency of rocket motors burning liquids to be more 
susceptible to high-frequency combustion pressure oscillations 
when the combustion chamber is long. The same mechanism 
can be employed for explaining why large liquid propellant 
rocket motors, since they have small aspect ratios, are apt to 
excite the transverse and tangential modes of combustion 
pressure oscillation (6, 7). 


The Effect of Nozzle Shape Study 


The results of the study of the influence of the shape of the 
nozzle (see Figs. 7 and 8) appear to confirm the mechanism oi 
reverberating shock waves described in the foregoing as an 
explanation of the mechanism for sustaining the longitudinal 
mode of combustion pressure oscillation. Fig. 8 indicates 
that when the converging portion of the nozzle does not 
furnish a surface from which pressure waves can be reflected 
with their full intensity, the wave shape of the longitudinal 
mode of combustion pressure oscillation will not be the shock 
type. 


Steady-State Combustion Pressure and Equivalence 
Ratio 


The data indicate that increasing the steady-state combus- 
tion pressure from 52 to 130 psia caused the peak to peak 
amplitude of the shock type of oscillation to increase from 20 
to 35 psi. A larger amplitude is to be expected at the higher 
steady-state pressure because of the greater mass of unburned 
and slow-burning propellants that become available for com- 
pression by a pressure wave in a combustion chamber of a 
given volume. The resulting shock wave, therefore, released a 
greater amount of energy thereby feeding back more pressure 
energy into the wave. 

For all of the premixed gaseous propellant combinat ons 
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amployed, but more pronounced for the ethane-air and the 
hydrogen-air combinations, it was noted that increasing the 
steady-state combustion pressure increased the amplitude of 
4 transverse mode. If it is postulated that the mechanism 
by which the amplitude of a transverse mode is increased 
should be similar to that by which the amplitude of the longi- 
tudinu! mode is increased, the transverse mode must also con- 
st of shock fronted pressure waves. Future experimental 
work with larger diameter motors has been planned for a 
thorough investigation of the factors influencing the trans- 
verse modes. 


Conclusions 


The mechanism which initiates the high-frequency com- 
bustion pressure oscillations observed in the experiments is 
not completely understood. It appears quite likely that the 
combustion process introduces momentary pressure gradients 
which are capable of initiating small pressure disturbances, 
and the latter initiate the combustion pressure oscillations if 
certain conditions, all of which have as yet not been investi- 
gated. are favorable for their amplification. Factors having 
a bearing on the character of the oscillations that have been 
investigated at Purdue University are the combustion cham- 
ber length, nozzle shape, equivalence ratio, combustion pres- 
sure, and propellant combination. 

It nay be concluded from the subject experiments that the 
shock fronted pressure waves are a type of detonation wave 
traveling longitudinally in the combustion chamber. The 
wave is sustained by the extremely rapid release of chemical 
energy immediately behind the shock front. Thus, the shock 
front and the energy release are interrelated in that the shock 
front causes the rapid energy release, and the rapid energy re- 
lease serves as the driving force for maintaining the shock 
wave. 

The sinusoidal shaped waves discussed in the foregoing 
appear to be a type of low-amplitude pressure wave in which 
the energy release occurs at a relatively slow rate. It would 


seem reasonable to assume that the rate of energy release is 
related to the thermochemical properties of the unburned 
premixed gaseous propellants, the dimensions of the combus- 
tion chamber, and the pressure change effected by the travel- 
ing pressure front. 

Experiments are either in progress or planned at the Jet 
Propulsion Center, Purdue University, in which the effects of 
changes in the combustion chamber diameter and other com- 
bustion chamber conditions are to be investigated. 
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This analysis is concerned with the essential features of 
the chemical interaction between boundary 
layer and sublimed or ablated material from a body sur- 
face or gas injected through a porous wall. These effects 
are in addition to the “‘blowing”’ effect of the gas added to 
the boundary layer. The models studied are those for 
which the diffusion coefficients of all species are equal, the 
** Prandtl and Lewis moduli are taken as unity, 
Solutions are 


a gaseous 


“frozen 
and reactions proceed at an infinite rate. 
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Boundary Layers With Chemical Reactions Due to 
Mass Addition 
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given for three problems: (a) A diffusion-limited irrey ers- 
ible reaction, typified by the oxidation of carbon to carbon 
monoxide, (b) a temperature-limited irreversible reaction, 
such as dissociation without recombination and (c) a re- 
versible reaction including complete dissociation and re- 


combination within the boundary layer. 
Nomenclature 


modified blowing parameter, 2( pv) pelt. 


cy = skin friction coefficient, 27,./ peu.” 

Mu( OH ‘Oy 
pete(H. — Hw) 
concentration by weight of species 7, p:/p 
Cp = specific heat, 0h/O7 


heat transfer coefficient, 


659 


| 
7 
} 
| 
‘ 
us- 
ALK 
a 
ns 


oh oc; 
= frozen specific heat, Cicp, = hs or 
Cp, = specific heat of a single species 
D = mass diffusion coefficient 
E = internal energy mm, 
h = enthalpy of mixture, Cihi 
h; = enthalpy of species 7, E: + (p/p): 
H_ = total enthalpy of mixture, + (u?/2) 
k = thermal conductivity 
i = latent heat of sublimation 
Le; = Lewis modulus, 7 
m; = mass source of species 7 per unit area 
M = product of molecular weight and stoichiometric molar 
reaction coefficient 
n = stoichiometric molar reaction coefficient 
p = pressure 
P = product species of a chemical reaction oe = 
Pr = Prandtl modulus, | 
q = heat transfer rate per unit area 
R = reactant of a chemical reaction 
Sc; = Schmidt modulus, p/pD; 
S = sublimate or injected gas 
7) = temperature 
u = velocity in z-direction 
ua = velocity ratio, u/u. 
v = velocity in y-direction 
w = mass source rate per unit volume 
x, = longitudinal coordinate along surface 
y = coordinate normal to surface 
\ = fraction of sublimate or injected gas entering into reaction 
a = dynamic viscosity 
p = mass density 
7 = shear stress, u(Ou/dy) 
+ = shear ratio, 7/7. 
Subscripts 
4 
C = carbon 
CO = carbon monoxide _ 
crit = critical 
e = boundary layer outer edge, ‘‘external’’ 
g = gas phase 
i = species identification index 
j = source (or sink) location index 
O = oxygen 
P = product 
R = reactant 
ref = evaluated at reference a 
S = sublimate or injected gas 
w = evaluated at wall 
* = at the reaction plane, in the case valk a single reaction 
plane 
1,2 = refers to reaction planes 


HE recent advances in boundary layer theory which in- 

clude effects of dissociation, ionization, recombination and 
other chemical phenomena have almost entirely been re- 
stricted to the chemistry of air. (An excellent review of such 
treatments is given by Fay and Riddell (1).)4 However, 
with increasing interest in systems utilizing melting or sub- 
liming surfaces or injected gases, it is necessary to consider 
the fact that significant quantities of surface or injected ma- 
terial may enter the boundary layer and react with it. The 
case of surface reactions has been treated by Emmons (2) 
and others. 

The purpose of this paper is to present in a brief manner 
the essential effects of chemical reactions involving subli- 
mated or injected gases within the boundary layer upon the 
amount of heat transferred to the ablating or porous surface. 
In this treatment, the diffusion coefficients of all species are 
assumed equal and the “frozen” Prandtl and Lewis moduli 
are taken to be unity. This allows considerable mathemati- 


* Numbers in parentheses indicate References at end of paper. 
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cal simplification and yet it is believed that the essential! fea. 
tures of the phenomenon are retained, as Lees (3) and others 
have demonstrated for nonreacting boundary layers. Thy 
present paper is restricted to the case of laminar flow so that 
the number of assumptions may be kept to a minimum and g 
mathematical analysis may be carried out. Since there is n 
apparent reason why the chemical reactions treated |vereiy 
introduce new differences between laminar and _ tur) ulent 
flow, it is likely that the utility of the enthalpy-velocity and 


_ similar relations can be extended to turbulent flows with th 


same degree of effectiveness as has been done for non react- 
ing flows. 
In treating the chemistry of reacting boundary layers, one 


can attempt to use realistic approximations for reaction rates 


This requires either exceedingly complicated mathen atica 
procedures, or subsequent limiting postulates. For ex: mple 
Dooley (4) assumed the concentrations of species at t!« sur- 
face to be known and obtained some very informati.e re- 
sults. Another procedure is to assume infinite reactio:. rates 
(sharp flame fronts) so that concentration need not bv arbi- 
trarily specified and to proceed on that basis. Such ‘s thi 
treatment of this paper. 

Before one can treat a particular example to determiiie the 
heat input (and rate of ablation, if such is the case) it is neces- 
sary to know the skin friction coefficient c; (and consequently 
the Stanton number c,) which is applicable to that case. This 
wall shear is strongly affected by the blowing action t the 


wall in addition to the effects of chemical reaction treated 
herein. The former influence, however, has been we!! dis- 
cussed in the literature ((5) for example) and so will not be 
treated. In this paper the emphasis will be on the new as- 
pects introduced by the presence of chemically active species 

For certain problems, where species concentrations ar 
known at the wall and at the boundary layer edge, a mor 
concise approach may be used to obtain the heat trans/er to 
the surface; this is done, for example, by Bromberg and 
Lipkis in (6). The method of the present paper allows more 
complete knowledge of the state of the boundary layer to be 
obtained. 


Analysis 


_ The applicable equations for a steady laminar compressible 


boundary layer are 


ou ov ou Op 
pu + pu | 
i 


In Equations [8, 4] it has been postulated that the diffu- 
sion of mass due to thermal and pressure gradients is negli- 
gible. Gaseous radiation has not been considered in this 
analysis. 

The energy equation presented above does not include a 
term for “extraneous” energy sources, such as electrica! dis- 
sipation, since such processes are not considered here. [How- 
ever, it should be noted that the internal energy, E;, of each 
species does include all forms of energy which that species 
can potentially release by any process. Therefore, it includes 
the energy (“heat of formation’’) which is released in chemical 
reactions (species transformations) evaluated at some (:irbi- 
trary) reference condition. Although it may not be immedi- 
ately apparent that the heat of reaction is contained implicitly 
in the energy equation, it can be made to appear explicitly 
by introducing a specific reference enthalpy. Then the ditfer- 
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+ Batiation of 2 Cih; yields terms which can be combined with 


thers 
The @ieappropriate mass rate of formation w;, from Equation [4], 
yield a term wyh;, ret which is the heat of reaction at the 
t 
— werence condition (6). For the present paper, the energy 
, : quation in the form of Equation [3] is employed. 
~ Returning to Equations [1 through 4], it proves useful to 
a entify therein the Prandtl, Lewis and Schmidt moduli 
and 
h th 
eact- 
One 
rates 
atical 
npk 
Sur- 
Te- 
rates 
arbi- 
the 
oC; oc; re) OC; 
C:CeS- 
ntly §ltis apparent that considerable mathematical simplification 
This J rsults from postulating Pr = 1 and Le; = 1. Lees (3) and 
the futhers have also shown that these are not unreasonable phys- 
ated ff alapproximations, at least for nonreacting boundary layers. 
dis- Jn the presence of chemical reactions, the effects of depar- 
t be Btures from these postulates may be significant, but have not 
as- J een considered in the present analysis. Recalling that 
r to 
and -quations [5, 6] become 
Lore oH OH _ ( oH 
be pus pv dy 7 
oc; oc; oC; 
pus 2( [3] 
ible 


e 


Then, following Crocco, the independent variables may be 
hanged from (x, y) to (#, u). In that event Equations [2, 


7,8] become 
nf 20 (2 ac; cat (2) oc; 
ox \ Ox | ag” du) Ou ou 
[11] 


lf Equation [9] is now substituted into [10 and 11], the sys- 
tem becomes 


oc; oc; Op oC; 
pur —p ae =T + ee [13] 


OeroBE 


u 0 0 T 


These equations can be further simplified by postulating H 


and C; to depend upon velocity only; that is, 0H/oz = 


0C;/d% = 0. This postulate is often given the name “simi- 
larity.” It may be shown to be self-consistent although the 


question of uniqueness has not been demonstrated. If, in 
addition, one considers the so-called “flat-plate’’ problem 
for which the pressure is taken as a constant, there results 


[14] 
ou? 

ou? 

ac; 

D;{ 1 — |... [5 


Two facts are immediately apparent from this system: 

1 The total enthalpy profile is linear with velocity, inde- 
pendent of whether or not chemical reactions occur (pro- 
vided no “external” energy sources are present; e.g., electrical 
dissipation). This will enable the effect of reactions upon 
heat transfer to be simply accounted for by means of the 
species composition at the surface. 

2 The concentration of a given species is linear in velocity 
except in reaction zones.° 

Equation [15] may be integrated twice with respect to u to 


get 
du 0 


where the primes indicate dummy variables. If (0C;/Ou). is 
determined by evaluating Equation [16] at the outer edge, 
u = u,, and is resubstituted into [16] - 


[16] 


The inside integral of the double quadrature can be written 


0 0 (1 0 
T 


oy 


The latter form suggests an approximation for the reaction 
zone to be a sharp thin flame front, or an ensemble of such 
fronts. In that event, Equation [18] becomes 


u 

where m;; is the mass source rate of species 7 per unit area at 
the jth position between the surface and the coordinate posi- 
tion u. It should be noted that m = O everywhere but at these 
reaction planes so that the shear stress need only be associated 


* Implicit in this statement and in the development which 
follows is that the quantity (uw;/r?) is a function of u only, a 
condition which is assured by the postulate of infinite reaction 
rates. Only for this model have the conditions for separation of 
variables been considered; the general conditions for separ- 
ability have not been examined. 
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with these locations. The u at the upper limit of the summa- 
tion indicates that only sources between the surface and the 
position where the velocity has the value u are to be con- 
sidered. 


With this notation convention, the double integral can be 
written 


(dove ° T Ti 
and Equation [16] becomes 
oc; mij 
C; Ci, ) (u uj) [21] 


The concentration gradients needed for Equation [23] can be 
obtained from Equation [21 or 22] 


(2). 


yielding, for the heat transfer at the wall 


Mi; 


Cie — Cet [25 


0 
where = is the summation over all species and S is the 
0 
eacininiian over all reaction planes in the interval (: to u, 
The first term on the right, expressed in the usual {rm, js 
Crpele(H, — H.) or, more precisely’ 
oH 
(u = Ch Pelle [(Cihi)e — (Cihi)w] + 
oy gw a y. 
In the second term on the right of [26], 7» may be 1 placed 
by (c,/2)pu.2, which for Pr = 1 equals cip.u.?. On intro- 


ducing these relations, Equation [26] becomes® 


oT ij 


0 


= (ue — w)| + (2 - 


The concentration relations defined in - preceding section 
are now applied to determine the resulting effect on heat 
transfer. Since 


and Equation [17] becomes 


Ci Ci, w = 


Heat 


*) 
le} 9 


Tj 


The term in square brackets deserves some particular atten- 
tion. It is not the usual difference between external and wall 
enthalpy, but is the enthalpy difference from the stream to 
the surface conditions, of the external species, regardless of 
whether they actually exist at the surface. The second term 
is the reaction contribution, which will be shown to be equal 


- ar ac to the heat of reaction appropriate to the particular example. 
= E hi 5 The mass of a species passing any plane parallel to the wall 
is 
one may write 

oc; oc; d 
i (= hi m: = Cipv — pD; = Cipv — T [28] 

oy gw Cp oy w oy gw 


or, for Pr = 1 


ar oH ou aCi 2 


For a subliming wall recall that*® 


oT oT 
( oy ). w ( oy w (ev) 


6 This equation is correct if there are no chemical reactions at 
the wall. In the general case which may include surface reac- 
tions 


(3 — (pivshido, w = w w 


so that an additional term would appear on the right of Equation 
[24], namely 


[24] 


(pivihi),, (pv) solid.w hs, 
See also footnote pertaining to Equation [27]. Radiation from 
the wall, omitted here, could be included by adding the usual 
«oT term on the right of nv [24]. 


Tw oc; 


Inasmuch as concentration is a linear function of u, the above 
expression for m; may be rewritten in terms of species con- 
centrations at known locations. It is therefore of interest to 
determine these concentrations. 

The general solution for the concentration of a species for 


7 Note that c, is defined here in terms of the actual species 
present at the boundary layer edge and at the wall. 

8 In cases where reactions occur on the surface, the surface 
reaction terms are not to be included in the double summation 
of Equation [27], which applies to gas phase reactions only. 
The surface reaction terms are already accounted for, as noted 
in the footnote to Equation [24]. Another intuitively satisi ying 
approach, which yields the same result, is to conceive o! the 
— reaction as a gas phase reaction just off the surface (where 

= 0) and to treat it in the general manner described herein. 
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the condition of several planes of infinitely fast reactions 
gecurring in the boundary layer is given in Equation [22]. 
The diffusion boundary condition for the sublimate or in- 
ff jected gas is 
oc, OCs 
(1 Cs, w) (pV) w (2) ... [29] 
oy Jw oy Jw 
(26) § and, for each of the other species, is 
oc, oc; 
w(pv)w — pD; ( ) = Ci, — ( . - [30] 
1s the Oy u Oy w 
to u. f These boundary conditions arise from considering that, for 
rm, is § the su limate or injected gas, the mass rate of the species is 
the to'al mass rate at the surface, and, for all other species, 
the convection outward is just balanced by the diffusion in- 
ward 'o yield zero mass flow at the wall ((5), for example). 
From ‘he general solution we have 
laced 
ntro- ac; 1 /ou 1 
_fo7) § Comb ning [29 and 31] yields the concentration at the wall for 
the sublimate or injected gas 
1 — a; 
there B’ is defined in terms of the actual skin friction coeffi- 
cient rather than the more usual no-blowing skin friction 
coefficient; A; is the fraction of the sublimate produced at the 
wall which is consumed at each reaction plane j, that is 
—Msg, j 
= 
ten- 
wall # The wall concentration of all other species is found from 
1 to § Equations [30, 31] 
Ol 
orm 1 M; oe 1 — i; 
erm i 
Ci.0 = 1 Cie - | ——— |... [183 
jual Ms ( 7; )] [33] 
ple. 
vall § Where M; is the product of molecular weight and stoichio- 
metric molar reaction coefficient of species 7. Oo 
For a single plane reaction represented by 
28) nsS + npP....... 
the concentrations of the sublimate, reactant and product of 
— § the reaction and of any inert gas present, are 
B’ 1 — 
Cs, = | 1 — Aw [85 
1 Mr 1 — 
Cr, e — A+B’ .. 
B’ Mp 1 — tx 
to 
Cinert e 
Ci o = ee 38 
or +1 
es § Since in the case of a single reaction plane all of the subli- 
mate or injected gas which reaches that plane is consumed, 
4 an additional expression for C's, » can be obtained by setting 
y. | Cs.* = Oin Equation [22] and solving for C's, ». This yields 
d 
Bt 
Equating [35 and 39] yields an expression for the fraction of 


>. ae from [22 and 36]. For a single reaction plane, on 


the sublimed or injected material which reaches the (single) 
reaction plane 


Tx 


1 + Bix 


The concentration of the reactant at the reaction plane is ob- 


combining with [40] 


= ( ) 


Similarly, the concentration of the product at the reaction 
plane, from [22, 37], is 


Bi Me 


B’+1Ms 


1+ B’ 


(1 — te)... [41] 


i B’ Mp 
( tix ) [42] 
2° Ti 


It should be noted that in obtaining the expressions for the 
concentration of the reaction product ([37, 42]), the concen- 
tration of the product at the outer edge of the boundary layer 
has been taken to be zero. This is not true in those cases 
where the product exists in the outer stream (such as oxygen 
or nitrogen), but even in these cases it is permissible to “‘tag”’ 
the molecules which are products of the reaction, whose con- 
centration does go to zero at the outer edge, and to consider 
the remaining outer stream molecules of the same species as 
inerts which do not enter the reaction. The expressions above 
can, therefore, be used even for products which also exist in 
the outer stream. 

It is now possible to apply these relations to certain idealized 
reactions. 


Example 1. Diffusion-Limited Irreversible Reaction © 
Consider a combustion reaction typified by 


with infinite reaction rate in the forward direction and zero 
rate in the reverse direction, for the case in which the idealized 
critical temperature for the reaction 7’ is below the wall tem- 
perature: Terit < 7... The critical temperature is that be- 
low which it is postulated that the reaction cannot occur and 
above which it must occur. The concentration profiles in 
the a-coordinate can be either the profiles shown in Fig. 1 or 
those shown in Fig. 2. 

Fig. 1 represents the case of reaction at the surface. There 
are two fundamentally different situations which might lead 
to this: 

1 The wall temperature is appreciably below the “subli- 


mation temperature,” so that a solid phase reaction takes 


2 The sublimation rate of C is less than the stoichiometric 
quantity compared with the rate of supply of O by diffusion, 
so that all the sublimed carbon is immediately consumed at 
the wall by the reaction. 

Fig. 2 represents the case where the mass rate of carbon 
sublimed is greater than can be combined with all the oxygen 
atoms reaching the surface per unit time, so that the reaction 
plane moves out into the boundary layer to a point where 
stoichiometric ratios exist and both C and O concentrations 
go to zero. 
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From [28] the mass of O consumed at the reaction plane is 


Co, 


mo,* = Co, « w= Te 
Ue — Ux 


since the concentration of O at the reaction plane is zero 
in thisexample. In addition 


REACTION 
PLANE 


T 


Fig. 2 Concentration distribution for example 1 (gas-phase 
combustion) 


(OT /Oy)],, (Equation [24]), the heat of formation a pplic- 
able to Equation [44] is the gas phase heat of formation, re- 
gardless of whether the reaction occurs in the gas phase or on 
the surface; for the reaction considered, Equation [43), this 
is 29,000 Btu per pound of oxygen. Further, the value of 
the heat transfer potential is the same whether the reaction 
is C +O— CO orC + 4 0.— CO because, while the + 0 
reaction is more energetic, by an amount equal to the heat 
of dissociation of Oz, the enthalpy difference across the bound- 
ary layer for oxygen (the summation term applicable to 


Mc 12 7 oxygen in Equation [44]) is reduced by that same quantity 
ee ae a _ for the part of the oxygen which enters the boundary layer 
dissociated. 
a 28 Using the concentration relations previously derived, it is 
possible to calculate the blowing rate necessary to cause the 
reaction to occur within the boundary layer rather than on the 
For a single plane reaction, the second summation of [27] wall. The critical blowing rate, below which the reaction of 
becomes this example will be on the wall and above which the reaction 
o. ; will occur in the boundary layer, is obtained from [36 or 
Ma, * 41]. Itis 
and, for the present example, [27 ] becomes? [ Me 
= 
oT _ 12 28 oT 
(: ) = ChPe Ue |= CE e (hi, hj, w) + Co, e (io w + he, heo, .)| (: ) 


Thus, the heat transfer potential to the surface is the usual 
potential of the outer species plus the heat of formation of CO 
(calculated at the wall temperature) per pound of air. The 
solution is independent of the location of the reaction plane, 
whether at the surface or within the boundary layer, and is 
independent of the concentration at the wall of any species 
involved in the reaction. It is apparent that this same re- 
sult would apply if the reaction occurred over a reaction zone 
of finite width made up of many reaction planes rather than 
on the single plane treated here. 

It should be noted that, in view of the expression for 


® A similar result has recently been obtained by M. R. Denison 
and D. A. Dooley of Aeronutronic Systems, Inc. 
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2 
+ Co, e AB reaction, w, [44] 


2 


Chpetle [= C;, Ahi, hi, w) + 


which, for the carbon combustion described by Equation 
[43], is 


Berit = 0.17 


By setting Cr, x = 0 in Equation [41], the blowing rate neces- 
sary to push the reaction plane out to any position in the 
boundary layer is, in terms of the coordinate a 


= B 
1 — (1 + Berit) 


[46] 
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fig. 3. Concentration distribution for example 2(a) (excess of Fig. 4 Concentration distribution for example 2(a) (excess of 
reactant at Terit) sublimate at T eit) 

lic- 
re- . . > 
ron § Equation [46] indicates that the condition 
this 
er — 1 , 1 - move out to the point where the sublimate and reactant reach 
on 1 + Berit stoichiometric proportions, at which point the concentrations 


of both will go to zero (Fig. 4). 

Equation [27] still applies, but now allowance is made for 
the possibility that the reactant will not be used up at the 
reaction plane; the amount consumed will be the difference 
between that entering and that leaving the reaction plane 


| () § corresponds to B’ = © (orc; = 0). Since the principal pos- 
eat & tulate of boundary layer theory, 0/Oy >> 0/0x, becomes ques- 
nd- § tionable under the circumstance (Ou/Ody). — 0, the present 

to § malyses must be limited to reaction somewhat inside of the 
tity § plane where 


yer 1 —mr,* = — Rl = 


bis Ue 1 + Berit 
pDr — | | — 
the oy oy 


Example 2. Temperature-Limited Irreversible Reaction 


dew - 


the 
1 of In example 1 the idealized critical reaction temperature was For Pr = Le = 1, pD) = pw. On replacing 0C/dy by (0C/du)- 
ion § below the wall temperature so that the reaction would occur (Ou/dy) 
or § whenever an S molecule encountered an R molecule. The sii : 
whole process was limited only by diffusion and convection Te (=) ] 
activity. In the present example the critical temperature for , Ou Just Ou Jus 
45] teaction is above the wall temperature so that the reaction 
annot occur in the region between the wall and the idealized The concentration profiles are linear in u, thus 
ritical temperature profile. In other respects the postulates r ic — tee 
of the problem are unchanged: There is a single reaction = | (Cr Cae Ce 
plane; the forward reaction rate is infinite, and the reverse oo 
rate is zero. Two somewhat similar vases are considered : On substituting Cp, » from Equations [36, 40], and for Cr. « 
Chemical combination of the sublimate or injected gas with from [41], the mass of reactant consumed at the reaction plane 
44) the reactant from the outer stream, or dissociation of the sub- hecomes 
imate involving no chemical reaction with the outer stream. yo wy. 
(a) Chemical combination of sublimate with reactant. The 7 Mr Bil 
oncentration profiles in the @-coordinate are shown in Figs. : Tx Ms 
sand 4. As before, the sublimate disappears at the reaction =) 1 + Bae 
jlane. The reaction product goes to zero at the edge of the 
10 § boundary layer (by definition of the edge) but is finite at the As before 
wall. The reactant from the outer stream (usually oxygen) _ ee 
vill still be present at the wall if there is an excess of it at the ® } tise Ms ee Mp mp. * 
reaction plane (Fig. 3). If there is an excess of sublimate at ee Mr Mr 
es- § the critical reaction temperature, the plane of reaction will Oo 
the The heat transfer is 
(i [= ue ( Ms _ Mp — tx) 
= crpette| Ci, hi, e — hi, w) +> + w + As. w 
16 2 Mp Mr Mg, (1 + OF 


Ue? — tx) 
ele C; e h; hy w ion, w |... [47 
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where AH reaction, v, g is the heat of reaction per pound of R where there may be one or more dissociation products of the 


evaluated at wall conditions. This equation requires knowl- sublimate or injected gas. This is identical to the case con- 
edge of the location of the reaction a. An alternative form sidered above except that there are no reactants, only sub- 
requiring, instead, knowledge of the wall concentration of the limate and products. The mass of P; produced at the reac- 


reactant can be obtained by employing Equations [36, 40] e 


T ie? 
(: = C:, € (hi, hi, w) + + AF reaction w, RICr, _ (B’ + 1)Cr, a [48] 
ow a 2 


This is similar to the diffusion limited case of example 1 ex- 
cept that here the heat of reaction is modified by a quantity tion plane, following the similar calculation above fo: the 
which is a function of the position of the reaction plane in the 


2 mass of reactant, is 
boundary layer, rather than the constant quantity Cr, . of : 


example 1. If the critical reaction temperature is reduced “4 Mp, Bl — de) 

so that the reaction moves to the wall, the solution should : t Msg Kidioaas 

agree with that of example 1. In Equation [47] as the reac- 1") ae 1 + Bis 

tion moves to the wall, a — 0 and the multiplier of the heat atts 7 
of reaction becomes from which £ 


Tx — 


Ue — Ux 1+ Bix 


The heat transfer, expressed in terms of the sublimate, i- 


oT ue B’ (1 — 
= elle C; h; = hi w ion, w 
(: ) [= A tion, w, + Bie ] 


CePrlle Ce e (h;, hi. + ry + AB w, S [B’ (B’ + 1)Cs, [50] 


where AH reaction, w, s, the heat of reaction or dissociation per 
pound of sublimate evaluated at wall conditions, would nor- 
mally be negative, thus reducing the heat transfer. In the 
a case of a dissociation reaction, no reactant from the outer 
_ edge of the boundary layer is involved, so the remarks in part 


Equation [41] shows that for the conditions of example 1, 
with the reaction on the wall (Cr, « = O and a = 0) 


Me, 
Ms BY = Cp, (a) above regarding the critical blowing rate do not apply; the 
dissociation will occur at the specified temperature plane no 
Thus example 2 contains example 1 as a special case. matter how great the blowing rate. 

When the condition of Fig. 4 is achieved, the temperature An alternate interpretation of the results of examples | 
limitation is removed and diffusional processes govern. The and 2 should be pointed out. In example 1, the effect of the 
solution should then agree with that of example 1 for the con- combustion is to add a constant term to the boundary layer 
dition where the blowing rate (of sublimate or injected gas) enthalpy potential (Equation [44]); the combustion reac- 
is great enough to force the reaction off the wall. The blow- tion, from the point of view of the surface, looks like an in- 
ing rate necessary to cause the reaction to occur to the right crease in that potential (for an exothermic reaction). The 
of the plane of critical temperature (Fig. 4) is precisely that situation in example 2 may be viewed differently. These solu- 


given by Equation [46] (Ce, * = 0). Substituting this in 
Equation [47] leads to 


ar 2 
i 2 


Therefore, as long as the blowing rate is such that the reac- 
tion occurs farther from the surface than the critical temper- 
(1 — Ux) 


Pa oT M 
(: ") = ChpeUe [> C;, e — hi, w) + + (pv)w AF reaction 
9. w 2 


tions (for example, Equations [47, 50]) can be written in the 
form 


Mr 


The chemical reaction is of greatest importance when dx is 
small (the term approaches zero when wd approaches unity), 
in which case the denominator of the bracketed quantity ap- 
proaches unity. For the temperature-limited reaction, wx is 
constant so the last term can be approximated as (pv) ..L’, 
where L’ is the bracketed quantity above and is an only 
slightly varying quantity (defined as positive for exothermic 
reactions). Recall that 


ature plane, regardless of whether the critical temperature is 
greater or less than the wall temperature, the heat trans- 
ferred to the wall will be independent of position of the reac- 
tion plane and will be as given by Equation [44 or 49]. 

(b) Sublimate dissociation or other similar changes of states 
involving a heat of formation. Under the idealized condition 
that the reaction rate of the reverse process (recombination) : 


is zero, such dissociation may be typically represented by ~ ar ar 7 
+ Ps + ( ). w (ev) +( oy w 
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The example 2 solutions can now be put in the form 


oy 


oT 
L + (« = Dy Ci, (hi, — + we + (pv) L’ 
solid, w 


(pv )w (L — L’) + 

so that in this ease the chemical reaction or dissociation now 
looks like an increase (for an endothermic reaction) or de- 
crease (for an exothermic reaction) in the latent heat of the 
wall material. 


Example 3. Reversible Dissociation Reaction Mee 


Consider a simple reversible dissociation reaction, such as 
might be typified by 


where, for convenience, the dissociation is postulated to occur 
at temperature T'aiss within the boundary layer, such that 
H,O molecules cannot exist at a higher temperature, and re- 
combination occurs at a lower temperature Trecom below 
which the dissociation products cannot co-exist. The con- 
centration profiles are as shown in Fig. 5. 

The example remains valid if the two reaction planes are 
moved together so that there is a single temperature profile 
separating the dissociated and combined molecules. 


oy 


If Equation [21] is written for the concentration of He at 
the outer edge, it becomes oe 


MH 

2.7 

— Cu, » = ue — (ue — Uj) 
0 


= Trecom Tdis 
| | 
| 

= 
| 


H,O 


2 
Cj 


WALL | 2 EDGE 


Fig. 5. Concentration distribution: Example 3, reversible dis- 
sociation reaction 
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Ue 
a7) = ChPelle Ci e (hi, hi, w) 
solid, w t 2 


ve mo, mo. no 
(ue — uj) = — (Ue — + (ue — U2) = 0...... [53] 
0 Tj T2 
and 
> (ue — uj) = (ue — 1) + — = 
if 0 Tj T1 T2 


oT Ue i 
(: 27) = ChPeUe b (hi, hj, + = hi, w (Ue 
a Tj 
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Since all of the He is consumed at the recombination plane, by 
the values of Cu,, and are zero. is zero 
because H, is not present in the outer flow. Therefore, the — 
summation term for Hyis zero. Specifically 


Mee, , 2 
—— (Ue — Uj) = (u, — wu) + 
0 i 1 2 


1 


(ue — uz) = 0 


The mass source and sink rates of O2 and H,O at each reaction — 
plane are tied to that of H, by Equation [51] 


Mo 32 
Mx,0 36 
Therefore, it must follow that a 


The heat transfer to the wall (Equation [27]) is 


ue 
Clearly, since the summation mij(u, — u;)/7; of each spe-— 
cies is separately zero (Equations [52, 53, 54]), it follows that — 
mij 


> 


Tj 


uj) = 0 


That is, as long as recombination takes place before the dis- 
sociated products diffusing toward the wall can reach the 
wall, the dissociation has no effect on the heat transfer at the 
wall. This is true even though much of the dissociated prod- 
ucts diffuse outward, thus remaining dissociated and pro- 
viding a net cooling at the outer portion of the boundary 
layer. For simplicity, the example considered a single plane 
of dissociation and a single plane of recombination, but the 
conclusion can be seen to be valid even for distributed regions 
of dissociation and recombination by approximating the dis- 
tributed regions with many small sources and sinks rather 
than just one of each. At each of these planes of reactions, 
the summation term for He, and therefore the corresponding 
terms for the other reacting species, will be zero, so that the 
heat transfer to the surface remains unaffected by the disso- 
ciation and recombination. 


Concluding Remarks 


The analysis presented herein has treated the effect of 
chemical reactions within a laminar boundary layer for which 
it was postulated that the Prandtl and Lewis moduli are | 


» 
of 
)- 
the 
rs 
= 
19] 
3 
) wb’ ; 
Li, 


unity, and that chemical reactions are restricted to a narrow 
gone and proceed at an infinite rate. Under these conditions 
it was found that: 
1 If the reaction is irreversible and the critical reaction 
temperature is less than the wall temperature, the reaction 
is _ limited by the diffusion of reactant from the external 
_ stream, and the effect upon the heat transfer to the surface 
is to change the potential from that of the nonreacting (but 
blowing) case by an amount equal to the heat of reaction per 
unit mass of the external stream, independent of the location 
the reaction plane. 

2 If the reaction is irreversible and the critical temper- 
ature is greater than the wall temperature, and in addition 
the blowing rate is less than a certain critical rate (E quation 
aa, [45]), the reaction is temperature-limited and the effect is to 
é change the heat transfer potential in a manner similar to case 
1 but by a lesser amount, the decrease depending on the dis- 
tance of the reaction from the surface (Equation [47]). For 
_ greater than critical blowing, the reaction is no longer tem- 
a but will be limited by diffusion of the reac- 
tant from the external stream, as in case 1, and the results 
are identical to that of case 1. The solution is similar for an 
_ irreversible dissociation except that in this case there is no 
critical blowing rate. 

3 If the reaction is reversible but recombination takes 
place such that no reaction products reach the wall, there is 
no effect of the reaction on the heat transfer potential, inde- 
pendent of the location of the dissociation and recombination 
the boundary layer. 

_ The present analysis does not examine the effects of depar- 
ture from Pr = 1 and/or Le; = 1. One method for deter- 


mining these effects would be to expand the pertinent equa- 
tions (e.g., [5, 6]) in terms of (1 — Pr) and (1 — Le;) so that 
one could calculate perturbation solutions. In view of the 
possible import of these departures and the inability at pres- 
ent to carry over the conclusion from nonreacting boundary 
layer theory that these departures are not significant, it is 
apparent that the more generalized solution merits further 
effort. It is of equal importance to define the role played by 
the rates of chemical processes involved; they must, in actual- 
ity, be finite. The ultimate goal of an analysis would be to 
include these effects as well as those of realistic values of 
Prandtl and Lewis moduli. 


References 


1 Fay, J. A. and Riddell, F. R., 
Point Heat Transfer in Dissociated Air, 
Aveo Research Laboratory, June 1956. 

2 Emmons, H. W., “The Burning of Carbon,’’ The Ramo- 
Wooldridge Corp., Internal Memo., Aug. 1954. 

3 Lees, L., “Laminar Heat Transfer over Blunt-’osed 
Bodies at Hypersonic Flight Speeds,’’ Jer PROPULSION, voi. 26, 
April 1956, pp. 259-269. 

4 Dooley, D. A., ‘Ignition on the Laminar Boundary Layer 
of a Heated Plate,’’ Heat Transfer and Fluid Mechanics [nsti- 
tute, Stanford University Press, June 1957. 

5 Eckert, E. R. G. et al., ‘““Mass Transfer Cooling of a 
Laminar Boundary Layer by Injection of a Light-Weight 
Frozen Gas,’’ University of Minnesota Institute of Technology, 
Technical Note 14 (ASTIA AD-132 395). 

6 Bromberg, R. and Lipkis, R. P., ‘‘A Critical Examination 
of the Boundary Layer Equations Including Chemical Reac- 
tions Due to MassAddition,’’ Space Technology Laboratories, The 
Ramo-Wooldridge Corporation, GM-TR-0165-—00288, Jan. 1958 


“Theory of Stagnation 
Research rt 1, 


A review is presented of the essential features of the dif- 
ferential equations applicable in the boundary layer with 
mass transfer at the wall, including consideration of chem- 
ical reactions. A brief derivation is given, and the equa- 
tions are presented in the alternate forms useful for 
various applications. Special discussion is given of the 
wall boundary conditions because of possible confusion 
arising from effects of mass transfer at the wall. A gen- 
eral formal solution to these equations is obtained which 
makes it possible to determine the heat transfer potential, 
subject only to the requirement that Lewis number equal 
unity, for all possible boundary layer chemical reactions; 
no limitation is placed on reaction rates or other details 
of the reaction, or on the Prandtl number. In the limit- 
ing case where reactions go to completion, the potential 
may be written directly. As examples, solutions are pre- 
sented for the case of combustion of carbon to CO and for 
the dissociation and recombination within the nonenay 
layer of an injected material. : 
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R. BROMBERG! and R. P. 


Transfer in Boundary Layers With Chemical 
| Reactions Due to Mass Addition 


LIPKIS? 


Nomenclature 
C; = concentration by mass of the ith species, 
pi/p 
Cp = specific heat at constant pressure, (0h/07')» 

tp = “frozen” specific heat, Cicp, 


oh oC; 
AC; = — Cis 


mass diffusion coefficient of the ith species 
in the total gas mixture 


E = internal energy 

f = function defined by Equation [25] 

h = enthalpy of mixture, YCihi 

h’ = = — hi.w) 

Ah, = heat of reaction 

Ahy.j.1.0 = heat of reaction of chemical reaction j, 
per unit mass of species 1 of reaction j, 
evaluated at wall conditions 

Ahs = heat of sublimation 

AH = heat transfer potential, defined by Equ:tion 


“frozen’”’ thermal conductivity 
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L = 


= operator, defined by Equation [19] 


Le = Lewis number of the 7th species, pDiCp/k 
= molar coefficient in balanced chemical 
equation 
0, CO, C, ete. |= chemical species 
= pressure 
Pr = “frozen’’ Prandtl number, ué,/k 
Qe = heat transfer rate per unit area at the wall, 
oT Oh’ 
oy Pr oY — 
= heat transferred 
= recovery factor ao) 
= time 
= velocity components in the x and y direc- 
tions 
= volume 
= mass source rate per unit volume of the 7th 
species 


= longitudinal and normal coordinate 

= chemical species entering reaction 

= reciprocal of formula weight (product of 
molecular weight and molar coefficient ) 

= dynamic viscosity 

= mass density 

= dissipation functic tion 


Or 
Subscripts 
e = boundary layer edge 
q 


a = gas phase 
a = refers to chemical species or energy mode of 
a single species 


I 7 7 = refers to nonreacting species 

j - = identifies the chemical reaction 

0 = oxygen, or per unit mass of oxygen 

= reaction 4, 

8 = surface material, or per unit mass of surface 
material 

w Pr = at the wall, or at the wall condition 

1, 2, etc. = wh seg number identifies species 


3; 9, first number identifies the chemical reac- 
tion, second identifies species in the reac- 


OME schemes for protection against high rates of heat 
transfer for high performance airplanes and missiles in- 
volve injecting foreign gas into the boundary layer or ablation 
of the surface material. The opportunity is great for this 
added material to react with the boundary layer gases, or to 
undergo dissociation or other changes involving release or 
absorption of energy, with the possibility of appreciably af- 
fecting the heat transfer to the surface. Several authors 
have considered various aspects of the chemically reacting 
boundary layer, for example Hirschfelder (1)? and Dooley 
(2); the emphasis in this paper is on those reactions involving 
gas injected through or vaporized from the surface. Refer- 
ence 3 also treats this problem but with a different and more 
restrictive approach. 

If the equation for heat transfer to a surface is written 


Introduction 


it is quite useful to examine the Stanton number c, and heat 
transfer potential AH, separately. It has been demon- 
strated by many authors that c, has a relatively weak de- 
pendence on the details of the boundary layer profile (provid- 
ing that some approximate account is taken of the property 
variations, as for example by evaluation at an appropriate 
reference enthalpy). Therefore if, as shown later, the po- 
tential itself can be determined in a rather general form, and 


Numbers in parentheses indicate References at end of paper. 
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if c, can be determined to a sufficiently good approximation 
utilizing existing solutions for the nonreacting situation (in- 
cluding the blowing effect of the gas entering the boundary 
layer from the surface, now available from solutions by various 
authors), we have immediately at our disposal useful working 
relations for the calculation of heat transfer with chemical 
reactions and gas injection. 

In (3) the heat transfer potential with chemical reactions is 
obtained, subject to the restrictions of Prandtl number equal 
to unity and zero pressure gradient, by determining the de- 
tailed boundary layer profiles in terms of the velocity, using 
the Crocco transformations. In this paper a more general 
approach is used to obtain the potential for particular cases. 
The Lewis number is restricted to unity, but the Prandtl 
number and pressure gradient are unrestricted. 

In addition to examining certain features of the governing 
agent considerable attention is given, as noted above, 
to the determination of the heat transfer potential with 
chemical reactions and surface gas injection. 

It is believed that the derivations leading to the final for- 
mulation, Equation [36], for the potential are new; however, 
except for the specia] formulations given, the specific results 
of the examples presented are not completely original, see 
for example, (3).4 : 
Brief Derivations of Appropriate Boundary Layer 
Equations 


Energy Equation 


The energy equation describing the compressible bound- 
ary layer in the absence of chemical reactions and cross diffu- 
oh 


sion effects is® 
oh ra) oT 
pu — + pp— =—(k 
oy 


ox oy oY 
where we have dropped certain terms as in the usual boundary 
layer approximations, and @ is the energy addition per unit 
volume by viscous dissipation. It will be convenient in the 
remainder of the paper to combine the last two terms into a 
single term, ¢’. For the steady-state boundary layer, the 
equation of continuity is 


Now, keeping in mind that k(07'/dy) refers to heat transfer 
by molecular exchange processes and does not include any 
energy transfer accompanying mass transfer (even of a single 
species in the case of a multicomponent gas), we can im- 
mediately write the energy equation for a gas composed of 
many species 


pins + hy = 


where w; is the mass rate of chemical formation per unit 
volume of species 7, and where represents the rate at 
which enthalpy is “soaked up” by such formation.’ Note 
that this term acts like the term convecting enthalpy out and 
is thus added to the terms on the left. Any effects of chemi- 
cal reactions will be associated with the species conversions 


4 Similar results have recently been obtained by L. Lees in 
“Convective Heat Transfer with Mass Addition and Chemical 
Reactions,’’? Third Combustion and Propulsion Colloquium 
AGARD, NATO, Palermo, Sicily, March 1958. 


5 The equations here are written for a flat plate; the extension 
to a curved cylindrical surface or body of revolution is 
made. 

6 The subscript 7 is normally taken to indicate the chemical 
species, but more generally can be taken, without change in what 
follows, to distinguish between various energy states of a given 
species as well as between species. 
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(or mass production) represented by the w; terms. When ac- 
count is taken of all reactions only Ah’s of formation (+ and 
—) will appear in the term Zw,hy. The k used here is the 


“frozen” value for the gas mixture (4). p,u; and p,v,; are the 
existing mass fluxes in the x and y directions for the species 7, 
whatever might be the causes of such flow (convection cr dif- 


fusion). 
The species continuity equation is aait, 
ox oy 


and £ w; = 0, so that Equation [3] still applies. 


It is convenient at this point to introduce a particular 
reference value for the enthalpy of each species. We define 


a quantity 
where hi, is the enthalpy of species 7 at wall conditions’ taken 


to be constant in this analysis, and substitute this into Equa- 
tion [4] yielding 


S44 


= hi — 


- no apparent absolute enthalpies appear. 


or, on subtracting the total (global) continuity Equation [3] 


re) oT oc; 
+ 2p > )+o 


It is important to note that by introducing the reference 
enthalpies in this manner a source term, ¥ Wii», appears, 


thus introducing the heats of formation in the gas pliase 

evaluated at reference conditions. Also, it is quite clear that 

all enthalpies enter as differences, taken separately for cach 

species. Inasmuch as the concept of enthalpy is meaningful 

only in terms of enthalpy differences, species by species, it is 
useful to have equations (such as [9 and 9a]) where tliese 
enthalpy differences are demonstrated explicitly and w ere 
The heat of forma- 
tion is, of course, itself a difference of enthalpies. 


Diffusion Equation 


If Equations [8] are substituted into the species contin ity 


_ Equation [5], the expression 


OF 
i dy 


re) re) 
— C; = D 
(Cipu) + (Cipv) 2( ) +u 10) 


is obtained, dropping the x derivatives of the diffusion term, 


oh;’ oh;’ of. or 
i ox oy oy 


Multiplying [5] by h,;’, summing over species 7, and adding the 
resulting Equation to [6] yields 


re) Oo. re) oT 
— (piushi’) + — (pivshi’) + Wihiw = — (: =) 
Oz Oy oy 


We now write the species mass velocity (flux density) in 
terms of the mean mass flow plus a diffusion term® (relative 
to the mean flow) 


oc; 
piu; = Ci(pu) — pD; — 
ox 


oy 


where C; is the concentration (mass fraction) of the ith spe- 


cies. Substituting in Equation [7] we obtain 
— (pudy Cis’) +— (pv Cihi!) — (x pb, 
ox oy i Ox ‘ ox oy 


Transposing terms, noting that =h’, and neglecting x 


derivatives of diffusion terms, we have 


re) re) re) oT 
— (puh') + — (pvh’) + wihiw = — + 
ox oy oy 


7The choice of reference value is a matter of convenience. 
Choosing wall values simplifies the statement of boundary con- 
ditions because h;’ = 0 for all species at the wall. 

5 We have neglected diffusion due to all gradients except species 
concentration gradients (thermal diffusion, pressure diffusion, 
etc.). Also, note that the D,’s are the same as the “‘number’’ 
diffusion coefficient used in kinetic theory. See for example 
(6). The D; for a particular species is the diffusion coefficient 
for species 7 in a background of all other species present. It can 
be computed by resorting to complex analytical techniques, and 
in principle can be measured experimentally (by suitable sam- 
pling, plus chemical or mass spectrographic techniques) in the ac- 
tual flow situation. 


670 


ac; A ac; 
= Ci(pv) — — 


This may be combined with the total continuity Equation 


[3] to yield the alternate form 


In general, in addition to the diffusion equations, the chemical 


re) oT 
pD; i’) + Wihi,w ar) = 
oy i oy 


equations and rate reaction equations are necessary to deter- 


mine the w; terms. = 
4 
Momentum Equation 


The momentum equation for the boundary layer is not al- 
tered by mass diffusion, except insofar as the density and 
viscosity are affected by the concentration changes, and re- 
mains as 
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Pr oy oy solid, w 
Energy Equation 
ti t When writing the edge boundary condition, it should be re- 
ow iy valled that the enthalpy for a particular species at the edge 
erence speck f of the boundary layer is not constant; rather, the total en- 
pears, ff understood by use of Equations [4 and 5]. aquation [5] thalpy including any chemical energy due to species conver-_ 
is multiplied by h; and summed over 7 and then combined sions is constant at the edge 
plase ff with [4], thus eliminating the “source” term, then the equa- ay 
that § tion applicable for the region near the wall can be written : 
r cach (because all uj. = 0) as Diffusion Equation 
ingful 
For the surface material 
3, 1 fo) re) oT 
these oy\ oy\ oy (pivi)ims,w = (Pr) w 
orma- Keeping in mind that a coordinate system has been chosen From Equation [8], we have at the wall 
(implicitly) in which the surface is fixed and mass flows ee ; 
through it (including the wall material flowing through the 
surface) Equation [13] can most conveniently be integrated ; Oy Jw 
across the interface by employing a fictitious boundary a small nathan oan 
muity distance € out from the wall in the gas stream, such that no 
reaciions occur between the boundary and the wall. Clearly, , dCs 
for any reactions not involving the wall material, for ex- (pr)w(1 — Cs) = — | pDs Oy ers Senne) 
10) ample recombination of atomic species from the external flow, ; 
: the temperature and property variations are so small over the For all other species, excluding those reacting at the wall, 
Vert interval that the contribution to the wall heat transfer will (p:)w = 0 and oe. 
— be the same whether the process occurs on the wall or a aC 
distance away from it. (The kinetics of the reaction, on (Cipv)» = (oo. 
the other hand, may well depend on whether it is a surface OY Jw 
[6] reaction or not.) Similarly, reactions involving the wall ma- 
terisi] behave the same thermodynamically whether the reac- Development of Alternate Form 
tion is at the fictitious boundary or on the wall. That is, in ar 
one case vaporization occurs first, and the reaction takes place For some purposes it is desirable to write Equations [9 and 
in the gas phase at the distance ¢ from the wall; in the other 11 ] in different forms. In particular, the introduction of 
case the reaction occurs on the wall in the solid phase, and the ratios of diffusivities (the Lewis and Prandt] numbers) is use- 
7 gaseous product leaves the wall and crosses the fictitious ful. We note that" 
boundary; the end result is the same regardless of the path of oe gh. ae ac, 
the process as long as the end temperatures are the same.® — 4 
In summary, all surface reactions behave as gas phase reac- or oy 
bade tions very near the surface and can be so treated, at least as oT > OC; 
far as heat transfer considerations are concerned. Cy oy + hi 
Now letting the fictitious boundary approach the wall as ’ ; 
¢— 0 and integrating across the interface yields where ¢, is the “frozen” heat capacity (that is, the heat capac- 
se i ity summed over species 7 and not accounting for species pro- 
ll (= puts) — (pv)whsotid, w = (: 7) = (i =) duction or conversion due to chemical reactions). Then the 
i Oy g.w Oy solid, w first term on the right in Equation [9] can be written 
(the contribution of @’ is zero when integrated across the o/, o7 ofk or 
ical interface). Because we have in effect transferred the surface by 4 ay) dy ue? oy 
reactions to the gas region slightly away from the surface, the * & : P 
only mass flow at the interface is that of the surface material. (J oC; i) 
The left side of the above equation then becomes ( Oy\ Mp OY bly i (OY 
(h h ) wh I h Ipy diff he I oe 
— Nsotia, Where the enthalpy difference is the heat o ci 
sublimation, Ahs. From the following section (Equation 
[18}) it can be seen that at the wall, where hj, = 0, Ofu oh’ , 
( on) dh’ oy (Pr oy Pr oy 
k—) =(= — 
oy J g.w Pr OY where we have defined ‘frozen’? Prandtl number in terms of 
9] It is this quantity which will be referred to as the wall heat a pe tg > the regular Prandtl 
transfer rate, g,; in the discussions to follow, special emphasis 
number, using the complete c, and & (including mass dif- 
will be placed on its evaluation. We can now write the fusion-carried energy) as is done, for example, in (1). How- 
boundary condition in a form more convenient for use here By 
as!0 : ever, use of the frozen properties has the virtue of allowing 
easy identification of the phenomena involved in the diffusing 
and reacting bounday layer. — 
* A more direct, but somewhat more involved, approach for It will be helpful to adopt the notation : 
handling surface reactions is to treat them as they actually occur, OS 
|. on the surface, in which case the result of integrating Equation av ) av ) ofu av ) 
d [13] across the interface is to yield terms of the form (pv); Ah,.i.w LC ) = pu— + pv = 2(£ oe ... [19] 
: by the methods of the section on “Determination of the Heat : Ox oy oy\Pr oy 
Transfer Potential,’’ for both recombinations and reactions in- : 
volving the surface material. The factor (pv); is then evaluated The diffusion Equation [11] now becomes (after replacing 
using Equation [8] and the conditions imposed in the above 
named section to yield terms for surface reactions identical to — 
,] those for gas phase reactions (see Equation [36] ). -« 11 Procedurally, we are following the scheme used by Lees (4) 
“ Radiation effects have been neglected throughout. in using ¢, and the ‘‘frozen’’ Prandtl number. 
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C; by Ci — Ce = AC, where C;,» is a constant) 


oACc; oAC; ofu Ac; re) 

L(AC;) = pu + po (Le; — 1) +w (20) 

and the energy Equation [9a] becomes 


The source term may now be eliminated between Equations 
[20 and 21] by multiplying [20] through by A;,., summing 
over all 7’s and adding to Equation [21]. There results!* 


re) oAC; 
L(h’ +20 ACihivw) (Les — 1) + 
OyLPr oy 


[22] 


As before, enthalpies appear as differences only (it will be 
shown later that the term involving 2h; in Equation [22] 


represents the heats of the chemical reactions which occur). 


Determination of the Heat Transfer Potential 
(Case of Le; = 1) 


Without completely solving the various equations it is 
possible to obtain the form of the heat transfer potential, and 
this will now be done.!* Because of the difficulty involved in 
handling, in a simple way at least, the case of Le, ¥ 1, only 
the case where this parameter is unity will be presented. Al- 
though one can argue, from previous experience with the non- 
diffusing and nonreacting boundary layer, that the effects of 
Pr ~ 1 are not appreciable and can be approximately ac- 
counted for by simple formulations, it is not clear that this 
will be true for Le; # 1. The question of the importance of 
such deviations must be pursued, possibly by perturbation 
methods—or possibly by machine calculations. 

In the examples to be given ¢’ is taken as zero with the idea 
that flow work and dissipation effects can be added to the 
heat transfer potential by use of an appropriate recovery fac- 
tor times u,2/2. On taking Le; = 1 the combined energy 
Equation [22] becomes 


Lh! + = 
The argument of this equation may be normalized so that the 


limit values are zero at the wall and unity at the edge of the 
boundary layer. Then we have 


h’ + AC 
= [24] 


The solution to this is some func- 


for the energy equation. 
tion f(x, y) where 


and f is zero at the wall and unity at the edge. We now have 


h’ + = + [26] 


12 The term @¢’ is carried in this form and not combined at this 
point with h’ to form h’ + u?/2, because it is simpler in the 
— analysis to add the appropriate terms to the solution 
ater. 

13 Some of the concepts used in the analysis presented here are 
similar to those used by Zeldovich (5) in an analysis of the flame 
combustion problem, although details of the analysis here are 
different. The applicability of these concepts to the boundary 
layer problem was pointed out to the authors by Lester Lees. 
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Differentiating with respect to y and evaluating at the wall 


yields 
oh’ C; 

Oy J OY] i of 
It is fruitful at this point to introduce the chemical rea ‘tion 

equation, or set of equations, representing the reaction «r re- 

actions occurring in the boundary layer. These equa‘ ions 
are of the form 


Nj = 0 
P 


where X;,, refers to particular chemical species (as O, CO. Os, 
Nz, ete.) and n;,; are the molar coefficients in the cheniical 
equations. (The first subscript 7 identifies the chemica! re- 
action; the second 7 identifies the species of the jth reaction.) 
From these equations can be computed values of the coef- 
ficients for the volume mass production rates, such that 


= 2Wj,2 = Aj,3Wj,3 = ... 


where a is the reciprocal of the formula weight (that is, the 
reciprocal of the product of molecular weight and the molar 
coefficient). For each chemical reaction, the volume mass 
production rates of each species of the reaction can be written 
in terms of the volume mass production rate of one particular 
species of each reaction 


q 

Qj,2 Qj,3 


We are now ready to use these relations. Equation [20] 


can be written, for Le; = 1, as 


= 


W;.1, ete. 


Qj,2 


The first and last forms can be subtracted, yielding 


Aci.) 

for which the solution in normalized form is a gain f, as in 

Equations [24 and 25], and 


Qa 
L(ACi,2) = = 
Q@ 


sothat 


Qj,2 


Differentiating with respect to f and rearranging terms we 
have at the wall 


0C;,1 Q;j,2 0C;,2 
— ! = ACj2. — .. [28 
( of ). a ( of 


etc., so that the quantity required in Equation [27] for any 
species of reaction j can be expressed in terms of one particu!ar 
species (¢ = 1) of that reaction. 
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The quantity (OC;,/Of),. will later (Equation [32]) be shown 
to be proportional to wall concentration for any nonsurface 
gecies; then the terms for any species which does not appear 
at the wall or the edge will drop out of the summation of 
Equation [27]. If then the species entering each reaction 
are grouped together and the reference species (¢ = 1) of each 
reaction chosen as any of the species remaining, Equation 
[27] becomes 


The enthalpy groupings are simply the heats of each reaction 
so Equation [29] becomes, on multiplying by yu/Pr 


From the usual nonreacting boundary layer theory, by defini- 
tion of the Stanton number cp, 


the 
\olar 
nass 
tten 
ular 


20) 


It is clear that to determine c, in the presence of chemical re- 
actions, due account must be taken of the effects of variations 


ah’ of 


OY J ow Pr Oy] of 


(To show these relations, Equations [20 and 25], written for 
an inert material, lead to 


AC,/AC),. = f 


where J denotes an inert species. From the definition of c, 


of = ChPeul 


ac; 
aj i = inert Of Ju 
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species 


thus 


( | 
= ChPelle 
Pr ACre OY /w 


On introducing Equation [16 or 15], setting Le; = 1, and 
manipulating algebraically, Equation [34 or 35] results.) 
The inert species therefore drop out of [33], and it becomes 


oh’ 
qu = (4 ) = expe +r + — (BY + + — (BY + ..- [36] 


Pr oy 


in p and the transport properties to determine the entire flow 
field. 

It remains to evaluate the terms involving (OC/Of),. The 
expanded quantity to be evaluated is 


Pr dy] ow\ Ff Jw Pr oy Jw 


For all species except the surface material this becomes, upon 


employing Equation [16], 


1 
Ve cm), 
CrPelle 


[32] 


where B’ is the modified blowing parameter, defined by B’ = 
(pv) .0/ChPele, and Le; has been taken equal to unity. The 
corresponding term for the surface material is 


CrpeueB’(1 — Cs,w) 


If such terms are reinserted in [30] this equation becomes | 


oh’ 
(4 ~) = Crpelle + Bhs (B’ + DC + Ah,,;,s.w[B’ (B’ + sie] +. 


Pr dy 


where the second summation applies only to those reactions, 
if any, where the surface material is chosen as the reference. 
Note that the quantity ¢’ has now been taken into account 
in Equation [36] by including in the heat transfer potential the 
term r(u,2/2), where r is the recovery factor and is equal 
to unity when Pr = 1. In view of relation [34] the quantity 
(B’ + 1)Cj1, in Equation [36] is the edge of boundary layer 
concentration of the fraction of species 1 which reached the 
wall without reacting. The bracketed term is then the ex- 
pected result: Simply the fraction of species 1 (measured at 
the boundary layer edge) which reacts, with heat release 

Another view of these final equations is gained by taking 
all terms with the factor B’ to the other side of the equation; 
since B’cpp.u, = (pv). these terms combine (recall Equation 
[14]) with the latent heat of sublimation term. The effect, 
in other words, of an outer edge species not reacting com- 
pletely but having a nonzero concentration at the wall, is as 
though the reaction did go to completion but that the heat of 


species 


where S’ represents nonreacting surface material. 

The significance of the terms in square brackets in Equation 
[33] can be seen by noting that for an inert gas the wall con- 
centration is 


[34] 
or, for the surface material when it is inert 
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sublimation of surface material was modified accordingly to 
account for the incompleted part of the reaction. 

We have now determined that the heat transfer potential is 
the usual value, involving only edge species, plus the term or 
terms representing the heats of formation for the reacting 
species, whether the reactions occur in the gas phase or on the 
wall, and whether they involve the wall material or are 
simply recombinations or dissociations of external or product 
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species. In determining the edge species enthalpies, it is 
important to note, as mentioned in the section on Boundary 
Conditions, that the total energy including appropriate 
chemical energy terms is conserved across any shock waves; 
furthermore, it is constant along the boundary layer edge. 

It is important to remember that the expression for the 
potential (bracketed term in [36]) is purely formal in that, 
generally, the concentrations at the wall are not known a 
priori. To determine the necessary wall concentrations re- 
quires (as has been noted earlier) solving the coupled equations 
describing the energy flow, momentum, species continuity 
and reaction rates. This means that due account must be 
taken of the detailed kinetics of each situation; for example, 
in the determination of the wall concentration (as contrasted 
with the heat transfer potential) it makes a difference whether 
reactions actually take place on the surface ora short distance 
off the surface. It is not the purpose of this paper to discuss 
the questions of finding the concentrations. Also it should 
be noted that for the limiting situation of sufficiently fast re- 
action rates (by which is meant here that the reaction even- 
tually goes to completion somewhere in the boundary layer 
or on the wall) certain solutions can be obtained as given in 
the next paragraph. 

The application of the above equations to two simple cases 
can be illustrated: 

1 The irreversible combustion on the surface or within 

the boundary layer of subliming carbon, as given by 


C+0O-—-CO 


under conditions that all the oxygen is consumed: Taking 
Equation [36] to be written for oxygen, then Cj. = Cow = 
0 and C;;,. = Co, The solution for the chemical reaction 
portion of the heat transfer potential is then Co,-Ah,,o,w 
which is the gas phase heat of reaction per unit mass of air 
evaluated at wall conditions. The effect of the reaction is 
simply to increase the usual enthalpy potential by this 
amount. The gas phase heat of reaction is the appropriate 
one here even for surface reactions, because the latent heat 
of sublimation is included on the other side of the equation 
(see Equation [14]). If there are intermediate reactions, 
such as the formation of O, from atomic oxygen and the sub- 
sequent reaction of O2 with carbon to form CO, the result is 
clearly the same as the direct reaction of oxygen atoms with 
carbon. In such cases it is simpler to treat the reactions in 
the latter form, referencing the reaction to the outer edge 
species, atomic oxygen, since its concentration is known at 
both limits, wall and edge. 

2 The dissociation in the boundary layer of a surface or 
injected material under conditions such that no dissociated 
products reach the wall (recombination takes place prior to 
reaching the wall): Recall that Equation [36] does not allow 
using as a reference species one which disappears at both 
limits, wall and edge, as do the dissociated products. Refer- 
encing both reactions (dissociation and recombination) to the 
surface material, the appropriate chemical reaction terms of 
Equation [36] are clearly equal and of opposite sign. It then 
follows that the enthalpy potential is completely unaffected 
by the dissociation and recombination, regardless of the loca- 
tion of these reactions in the boundary layer. If recombina- 
tion is incomplete because of slow recombination rates, a net 
effect will be felt at the wall depending on the quantity of dis- 
sociated products which reaches the wall, which in turn de- 
pends on the reaction rates or detailed spatial description of 
the reactions. If recombination is incomplete because the 
dissociated products do not diffuse to the wall, the effect on 
the heat transfer potential is also not zero, but the require- 
ment of Le; = 1 precludes this occurrence here. 

a It is interesting to note that for the two examples given the 


- it 


results are actually independent of any possible pressure 
gradient effects. Also, they should apply equally well to 
turbulent flow, presuming that the same arguments can be 
made regarding diffusion phenomena as have been shown em- 
pirically to apply to heat transfer effects in the nonreacting 
boundary layer. 

The approach presented here thus enables the heat transfer 
potential in the presence of chemical reactions to be written 
directly, without the necessity of solving the differential ¢qua- 
tions, for the class of problems where one of the reactar ts is 
consumed completely in the reaction. For all other case~ one 
must either know the blowing rate B’ and the wall conceittra- 
tion of the reference species, or the fraction of the edge :on- 
centration of the reference species which has reacted irre ers- 
ibly before diffusing to the wall. In general] then, one req tires 


knowledge of the reaction rates. | 
Le 4 


By examining the differential equations of the boun ary 
layer which allow for diffusion of reacting species and by re- 
taining enthalpy differences species by species and utilizing 
the concept of “frozen” specific heat, it is possible to identify 
the changes in the heat transfer potential due to chemica! re- 
actions of any nature, involving outer flow gas and/or surface 
material. Further, by incorporating species conservation 
relations of the reacting species into the energy equation, 


Conclusions 


— certain details of the heat transfer processes become apparent 
_ without completely solving the differential equations. 


For those cases where one of the reactants is consumed com- 
pletely by the reaction, the added effects of the chemical reac- 
tion on the heat transfer potential may be obtained directly. 
Two examples illustrate the method: (1) For the irreversible 
combustion of carbon to CO, it was shown that the heat trans- 
fer potential is simply the usual one plus the heat evolved by 
the reaction, independent of the locaticn of the reaction, and 
(2) for dissociation and recombination within the boundary 
layer of an injected material, it is shown that the heat transfer 
potential is completely unaffected. For all other cases, either 
the blowing rate and reference species wall concentration must 
be known, or the fraction of the edge concentration of a react- 
ing species which reacts irreversibly must be known. 

The general analysis leading to Equation [22] applies to the 
flat plate in laminar flow without restriction on Lewis number; 
the two examples cited follow from Equation [36] and strictly 
apply only to the case of Ze; = 1, but no restriction has been 


placed on Pr or pressure gradients. Although it has not been 
discussed fully, the general form of the solution should be ap- 
plicable to turbulent as well as laminar flow by extension of 
the present work in the manner in which previous laminar 
analyses have been applied to turbulent flows. 
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The burning program for a rocket powered missile mov- 
ing along a rectilinear path is investigated with the idea 
of minimizing the time interval necessary to reach a given 
point in space. The induced drag is assumed negligible 
with respect to the zero-lift drag; the atmosphere is re- 
garded as isothermal. 
the Calculus of Variations, it is shown that the totality of 
extremal ares is composed of zero-thrust sub-ares, sub- 
ares flown with maximum engine output, and variable 


By using the indirect methods of 


thrust sub-ares. Concerning the latter, and in connection 
with suitable hypotheses for the drag function, closed form 
solutions are obtained for the relationships between alti- 
tude. mass, time and Mach number. The boundary 
value problem is investigated for prescribed initial and 
fnal values of altitude, Mach number and mass. The 
analysis shows that, for given end conditions, two extremiz- 
One which minimizes the time and 
another one which maximizes the time. Attention is also 


given to the case where the final Mach number is free of 


ing trajectories exist: 


choice. A numerical example shows that, assuming the 
initial Mach number is zero, the optimum path includes: 
An initial sub-are flown with maximum engine output, an 
intermediate sub-are flown with variable thrust, and a 


final sub-arce flown by coasting. 


Introduction 


THEN a missile is flying in a vertical plane, two possi- 
bilities exist for modifying its trajectory: One is to 
act on the controls of the elevator; the other is to act on the 
controls of the engine. For a given set of initial conditions, 
there is‘an infinite number of trajectories which are mathe- 
matically and physically possible, depending upon the way 
the elevator or the engine or both are being operated-as a 
function of the time. It is in this sense that a missile, con- 
strained to travel in a vertical plane, may be referred to as a 
machine having two degrees of freedom: One is associated 
with the elevator’s position; the other one is associated with 
the power setting. Of course, it is assumed that the thrust 
direction, relative to the airframe or to the vector velocity, is 
given. 

The problem of the optimum operational performance in a 
vertical plane consists in saturating these two degrees of 
freedom so that the difference AG = G, — G; between the 
final and initial values of an arbitrarily specified function G 
of the generalized coordinates of the missile is a minimum. 

In view of the analytical difficulties associated with the 
previously indicated problem, the attention of most re- 
searchers has been attracted by comparatively simpler 
questions. These questions have been essentially originated 
by freezing one degree of freedom and studying the optimiza- 


tion of the other. For the category of problems of optimum 
burning program, the degree of freedom associated with the 
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elevator’s control is frozen by imposing some geometric or 
aerodynamic constraint on the trajectory (i.e., vertical flight 
or level flight or zero-lift trajectory, etc.); the other degree of 
freedom (the one associated with the maneuver of the engine 
control) must be determined so as to minimize the quantity 
AG. Previous analyses in this area of questions are now 
mentioned. Horizontal paths have been investigated in 
(1 to 3).3 Vertical trajectories have been considered in (4 to 
10). Moreover, near vertical paths have been analyzed in 
(11) in connection with thermal problems for the missile skin. 


Object of the Present Investigation 


The present research falls under the broad domain of the 
problems of optimum burning program. Trajectories which 
are rectilinear and either vertical or near vertical are con- 
sidered. An ideally isothermal atmosphere is assumed, in so 
far as this seems to be the only case, having engineering impor- 
tance, where solutions in a closed form are possible, thus 
leading to clear cut information on the qualitative nature of 
the optimum thrust programming. The brachistochronic 
problem (G@ = £) is investigated: It is the problem of the mini- 
mum time required to reach a given point in space. 

Use is made of a special coordinate system, described by 
Equations [5 to 9]. Such a system embodies a number of 
distinctive features: It enables one to carry out an almost 
straightforward determination of the optimizing condition; 
it is such that all Lagrange multipliers are constant along the 
variable thrust sub-are of the extremal solution; it is specially 
advantageous for the solution of the boundary value problem 
associated with the brachistochronic path; by the use of Green’s 
theorem (13) it enables one to carry out a sufficiency proof 
for the solution, in a case where the classical variational cri- 
teria due to Legendre-Clebsch and Weierstrass fail to 
produce any useful information on the minimal or maximal 
nature of the extremal path. 

The latter feature of this new coordinate system deserves 
to be underlined. As the analysis of (13) shows, for a given 
set of end conditions there are always two trajectories which 
yield a stationary value for the time: One is a solution of 
minimum time; the other one is a solution of maximum time. 
It is, consequently, of engineering interest to have an analyt- 
ical method available by which minimum time paths can be 
discerned from maximum time trajectories. 


Fundamental Hypotheses and Equations of 
Motion 


In this paper use is made of the following assumptions: 

1 The rocket powered vehicle is regarded as a particle of 
mass variable with the time. 

2 The thrust is tangent to the flight path. 

3 The equivalent exit velocity of the rocket engine is re- 
garded as a constant, independent of the operational condition 
of the engine. 

4 The engine is capable of delivering all mass flows 
(8) bounded between a lower value (8 = 0) and an upper 


3’ Numbers in parentheses indicate References at end of paper. 
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_ Fig. 1 Parametric representation of the engine characteristics 


_ § The aerodynamic lag is disregarded, i.e., the aerody- 
namic forces are calculated as in unaccelerated flight. 
The acceleration of gravity is a constant. 
The atmosphere is isothermal. 
_--&-~s« The trajectory is straight and inclined at an angle 6 
__ with repect to a horizontal plane. 
9 The induced drag is negligible* with respect to the zero- 
drag. 
7 In the light of assumptions 1, 2, 3, 8 and 9 the dynamic be- 
“a havior of the rocket powered missile is represented with the 
following set of differential equations 


. D — BV. 

= velocity 

= altitude above sea level 

= mass 


Ore) 


The dot sign denotes derivative with respect to time. 
In view of hypotheses 5, 7 and 9 the drag D is represented as 


acceleration of gravity 
path inclination 

drag 

engine mass flow 

e = equivalent exit velocity of the rocket engine 


D= 2 PoSC po(M)M?2 exp ( — 


where 


the Mach number (/) only. 


Transformation of the Equations of Motion 


The following set of non-dimensional variables is introduced 


go 7= hos [6] 
ie 


* Notice that hypotheses 8 and 9 limit the applicability of the 
present theory to vertical or near vertical trajectories. 
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f = ratio of specific heat at constant pressure to specific 
heat at constant volume (air) ~ 

Po = sea level pressure 

S = reference surface 

a = speed of sound 

Cpo = drag coefficient at zero-lift, assumed to depend on 


These variables are employed to transform Equations [1 to 
3] from the (t, h, V, m, 8) space into the (7, n, M, €. w) space. 
After assuming the Mach number (M) as the new indepe dent 
coordinate, extensive manipulations yield these results 


where primed symbols denote derivatives with respect to 
Mach number 

Tt’ =dr/dM =de/dM n’ =dn/dM 
The g;-functions and the y,-functions (k = 1, 2, 3) ar de- 
fined as 


bi 
= 
1 
= 6)....... [13] 
CpoM? 
yum, 
M 


=M ¥2(M, €) 


CpoM? + € 


exp (—n) = ¢3(M, €, [17] 


M 
= —e€— exp (—yn) = ¥(M, 


where M, = V./ais the ratio of equivalent exit velocity of the 
rocket engine to atmospheric speed of sound (M.,, therefore, 
does not stand for Mach number in the exit section of the 
rocket engine). 

It must be recalled that, in view of hypothesis 4, the 
actual mass flow 8 is a quantity bounded between a lower 
value (8 = 0) and an upper value (8 = Bmax). The non- 
dimensional mass flow yp, therefore, is also bounded between 
a lower value (u = 0) and an upper value (u = pmax). This 
information can be translated into the setting of the vari- 
ational problem by introducing the idea of parametric repre- 
sentation of the engine characteristics, already used by the 
writer in (3, 9 and 11). The non-dimensional mass flow p is 
represented as a function of a parameter a@ having the follow- 
ing properties (Fig. 1) 


<« p=0 
for « < + KM = Mmax 
d 
d« 


With this scheme, a@ is considered the independent parameter 
of the rocket engine and is allowed to vary between — 
and +o. The mass flow uw becomes a dependent quantity, 
varying between 0 and pmax, according to the scheme of Fig. 
1. Notice that the condition du/da = 0 represents either a 
coasting flight or a flight with maximum engine out)ut. 
On the other hand, du/da # 0 represents any other operating 
condition of the engine between the two limiting ones. 
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Variational Formulation 


The time interval Ar = 7; — 7; necessary to transfer the 
rocket powered missile from the initial point to the 
final point is calculated by integration of Equation [10]. 


Assuming that the initial time instant is 7; = 0, one obtains 


Attention is now focused on the differential Equations [11 
and 12], where go, W2, ¢3, W3 are defined by Equations [15 to 
18] and w = u(a). This set involves one independent vari- 
abl: (M) and three dependent variables (n, €, a). One degree 
of freedom is left, and an optimum requirement can therefore 
be imposed. 

Concerning the end conditions, it is assumed that the co- 
ordinates M;, ni, €:, ny, €& are prescribed, the subscript (7) 
referring to initial point and the subscript (f) denoting final 
point.6 With regard to the final Mach number M,, two cases 
are considered: The case where M, is given and the case 
where M, is free of choice. Since a’ is nowhere present in the 
expressions for J, J2, J3, no condition can be imposed on a. 
The end values for a (i.e., for the mass flow 4) must be deter- 
mined from the solution of the variational problem. The 
latter can be treated indifferently within the frame of ques- 
tions of Lagrange type (12) or questions of Mayer type (12). 

If the Lagrange formulation is adopted, the brachisto- 
chronic problem is stated as follows: ‘Among all sets of func- 
tions n(M), e(M) and a(M) which are consistent with Equa- 
tions [11 and 12] and with certain prescribed end conditions, 
to determine the special set which minimizes the integral 
(19)”. 


Euler Equations 


A set of Lagrange multipliers 4; = 1, Ax = Ae(M), A3 = 


\3(.V/) is introduced and thisexpressionformed 
3 


where J; is the integrand function appearing in Equation 
[19]; J. and J; are the first members of the equations repre- 
senting the constraints of the variational problem. 

Since there are three unknown functions, three Euler equa- 
tions must be written 


d (OF oF 
dM (27) 


In the above equations 2; = 7, 22 = €, 2; = a, and F is the so- 
called fundamental function defined by Equation [21]. Sim- 
ple manipulations yield the explicit form for the Euler equa- 
tions 


® In all variational problems where one variable (u) is bounded 
between a lower value (u = 0) and an upper value (ymax), sub- 
arcs » = O and sub-ares » = pmax must be considered in addition 
to those sub-ares which are solutions of a set of Euler equations. 
With the device of parametric representation of the engine char- 
acteristics the sub-ares = 0 and = umax are formally presented 
as Eulerian solutions. In this way Equations [23 to 25] hold for 
all sub-ares forming one extremal arc. 

’ To prescribe the initial and final values of » and ¢ implies to 
specify the initial and final values of the mass m, in view of 
Equations [6 and 8]. 
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Discontinuity of the Eulerian Solution 


As the Euler Equation [25] shows, the extremal arc is dis- 
continuous, being generally composed of sub-ares of equation 


and sub-ares of equation 
d 
d« 


Equation [26] represents a flight condition with a continu- 
ously variable thrust, as shown in the next sections. Equa- 
tion [27] offers alternatives: Either a coasting flight (u = 0) 
or a flight with maximum engine output (u = pmax). 


Corner Conditions 


Owing to the discontinuous character of the solution, the 
important conditions of Erdmann and Weierstrass must be 
satisfied (12) at each corner point, i.e., at each junction of 
different sub-ares. The details of the analysis are omitted 
for brevity and can be found in (13). The conclusions are: 
At each junction point the multipliers A, and A; must be 
continuous, and the latter must also be zero. This circum- 
stance evidently implies that a corner point may be located 
on the variable thrust sub-are of the extremal solution. 


Boundary Conditions 


For the present problem, the boundary conditions include a 
number of fixed end point conditions plus one natural condi- 
tion. This must be deduced from the general transversality 
condition of the Calculus of Variations (12) 


3 7 3 f 
OF oF 


which is to be identically satisfied for all systems of variations 
(6M, 6z,) consistent with the prescribed end conditions. In 
view of the fact that 


6M; = on; = 66; = bn, = =0......... [29] 
oF 

[30] 


the transversality condition reduces to 


oF , OF 
F- 


31 
On’ n ) f [31] 


For the particular case where the final Mach number is free 
of choice (6M; = arbitrary) Equation [31] yields 


oF OF 
a relation equivalent to 
[—¢i + — As(u + = O......... [33] 


The Family of Variable Thrust Sub-Ares 


In this section, the variable thrust sub-are of the extremal 
solution is analyzed in detail, because it plays a fundamental 
role in the solution of the variational problem. 


Since A3 = 0 
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at points of such a sub-arc, the Euler Equation [23] leads to 


After introducing the new constant 
and the definitions 
[37] 


the Euler Equation [24] is rewritten as 


where the subscripts to parentheses denote quantities to be 
kept constant when calculating partial derivatives.’ 


From Equations [43 and 11], the derivative of the non- 
dimensional] altitude 7 with respect to the Mach number (/) 
may be calculated. After extensive manipulations, one 
obtains 

dn M? + 22M.M + — 1) M2 

dM ~ M2 + [yM2 (x — 1) — 1]M — 2M, 
Analogously, the derivative of the non-dimensional] time (r) 
with respect to the Mach number (M) is supplied by 


dM yM.M? + [yM2 — 1) —1]M — 2M. M 
Equations [44 and 45] can then be integrated yielding? 
16} 
17] 


where C, and C; are integration constants, chosen so t at 


A(M) = B(M) = 0 at M = 1. The functions A(M) «nd 
B(M) are defined as 
Ki Ks Ke M Ky Ke KK; M 
where 


Further Transformation of Coordinates 


It is of interest to notice that Equation [38] can also assume 
the elegant form 


4 


a of coordinates from the (e, M) 
plane into the (ex, M) plane, where e, is a new non-dimen- 
sional variable, defined as 


Equation [39] suggests an immediate graphical procedure for 
finding the points of the variable thrust sub-arc: It consists 
of determining the Mach number M for which ¢ is stationary, 
for given values of C and ex. 


Explicit Solutions 
The relationship between Mach number M and non- 
dimensional mass € can be explicitly solved in terms of the 
latter variable if Equations [13 to 16, 36 and 37] are intro- 
duced into Equation [38]. The following result is obtained 


1 
id 


2yM. 


K,; = 


e=M 
1+ CyM, 

For the all supersonic domain the relationship between zero- 

lift drag coefficient and Mach number may be represented in 

the form 


where y and z are appropriate constants. In such a case 
Equation [41] is rewritten as ie 


1+ CyM, 


€ = 


7 It must be emphasized that there is one variable thrust sub- 
are for each value of the constant C. The latter must be deter- 
mined so as to satisfy the end conditions of the variational 


problem. 
8 In practice, z is a number bounded between 1 and 2. 


2yM. w= 
+ Mie + (52 
= x 1 
7M. 
K; = (x — 1)M ‘or [54] 
1 
Remarks 


For given values of x, y and M,, Equation [43] has the form 
e = e(M, C) and therefore embodies a single infinity of curves, 
one for each value of the constant C. After simple trans- 
formations, Equation [43] can also be written as 


€ M Cy 


=2-1+— —-——"___ 
yM* + M, M(1 + CyM.) 


(M — Ky 


With reference to the positive quadrant of the (e, M) plane, 
all curves of the family « = e(M, C) have a common poit, 
the one of coordinates 


® By simple transformations Equations [46 and 47] can be 
reconduced to a form similar to Equations [61 and 62] of (9). 
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Fig. 2 Family of variable thrust sub-arcs 


Its significance can be understood by calculating the instan- 
taneous acceleration along the variable thrust sub-are 
(M — K,) (M — 


= 
Y M? + 22M.M + M2z(z — 1) ie 


_aM 
gsin@ 


and by showing that V is negative for M < K,, zero for M = 
K,, and positive for M > K;. Asa conclusion, after entering 
the variable thrust sub-arc, the missile moves so as never to 
cross the line M = K, of the (e, M) plane. Since in practice 
K, is always less than unity, it seems logical to conclude that 
as far as the variable thrust sub-are is concerned the region 
of the (e, M) plane where M < K, has little or no interest for 
flight operations. 

For the particular case where = 2, M. = 10, Equation 
[43] is plotted in Figs. 2 and 3. In Fig. 2, the neighborhood 
of the point of coordinates (é, M) is indicated. In Fig. 3, 
the region M > K, is represented. As the analysis of (13) 
shows, all the solid lines of Fig. 3 are associated with mini- 
mum time solutions, while all dotted lines of Fig. 3 yield 
maximum time solutions. The special sub-are C = 0 is the 
transition line from minimum time solutions to maximum time 
solutions. 


Constant Thrust Sub-Arcs 


For the pieces of extremal solution which must be flown 
tither with maximum engine output (4 = pmax.) or by coasting 
(u = 0), the equations of motion cannot be integrated, 
in general, into a closed form. Approximate procedures are 
usu: lly in order, except in a few cases whose detailed dis- 
cussion can be found in (13). 
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Solution of the Boundary Value Problem for 
Fixed Final Mach Number 


In the previous sections general solutions have been de- 
rived for the problem of extremizing 7;. It has been shown 
that the extremal arc is composed of sub-arecs 4 = 0, sub- 
arcs f = max and sub-arcs to be flown with regulated thrust. 

Here the boundary value problem is investigated. It con- 
sists of determining that combination of sub-ares which satis- 
fies a set of prescribed end conditions. Assuming that the 
initial Mach number is M; = 0, the extremal path generally 
starts (11) with an initial sub-are 7A flown with maximum 
engine output (Fig. 4), continues with an intermediate sub- 
are AB flown with variable thrust, and terminates with a final 
sub-arc BF. The latter can either be traveled with a coasting 
flight (Case I) or with maximum engine output (Case II), 
depending upon the boundary conditions of the problem 
(Fig. 4). It is to be noted that, for given initial and final 
values of Mach number, altitude and mass, the end points / 
and F are completely defined in the (M, 7, €) space. The 
problem is, therefore, one with fixed end points. The six 
coordinates of corner points A and B are unknown and must 
be determined. 


Determination of Corner Points in the (ce, M) Plane 


The first step is to locate the end points J and F of the de- 
sired extremal solution in the (e, M) plane (Fig. 5). Assum- 
ing that the initial Mach number is M; = 0, the sub-are p = 
Mmax Starting at J is calculated by forward integration of the 
equations of motion. Analogously, the -two possible. sub- 
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arcs w = Oand w = pmax arriving at F are calculated by back- 
ward integration of the equations of motion. In this way 
several possible couples of corner points A and B are indi- 
vidualized, with a graphical procedure, in the (e, M) plane. 
Notice that there is one couple of corner points A and B for 
each value of the constant C characterizing the variable 
thrust sub-are of the extremal solution (Fig. 5). 


The A-Function 


‘The next step is to account for the altitude condition. The 
below-indicated functional is defined 


f + YredM ....... [61] 
1ABF 


The latter has the form A = A(C), if it is applied to the family 
of trajectories JABF of Fig. 5. 

The final step is to search, among all trajectories JABF of 
Fig. 5, that special path such that A(C) = 0. This trajectory 
evidently satisfies all boundary conditions of the variational 
problem. 


Numerical Example 


To illustrate the theory, a numerical example is carried out 
for the set of conditions” 


=0 n = 0 
M; = 0 nr = 25 m,/m; = 0.6 

[62] 


10 The symbol m, stands for propellant mass; m; denotes initial 
or take-off mass. 


680 


= 


é/y 


variable 


= 

M 

Fig. 5 Combination of sub-arcs in the («/y, M) plane 


After determining all possible couples of corner points A 
and B in the (e, M) plane, the A-function is calculated and 
plotted in Fig. 6. Since the extremal trajectory is to be such 
that A(C) = 0, a fundamental result is obtained: For a 
given set of end conditions there are two extremizing tra- 
jectories, the path /ABF associated with the variable thrust 
sub-are C = —0.06, and the path JABF associated with the 
variable thrust sub-are C = 0.68. 

Both solutions JABF include an initial sub-are 7A flown 
with maximum thrust, an intermediate sub-are AB flown 
with variable thrust and a final sub-are BF flown by coasting 
(Figs. 7 and 8). As shown in (13) the trajectory JABF asso- 
ciated with C = —0.06 maximizes the time (7, = 8.82); the 
trajectory JABF associated with C = +0.68 is the brachi- 
stochronic one, i.e., the path of minimum time (7; = 7.29). 


Solution of the Boundary Value Problem for Free 
Final Mach Number 


If the final Mach number M, is free of choice, two possible 
techniques may be used in solving the boundary value 
problem. 

1 One method is based on the transversality condition of 
the Calculus of Variations, represented by Equation [33]. 
This technique requires that the distribution of Lagrange 
multipliers be calculated at points of the sub-ares w = [max 
and « = 0. For the sub-ares » = 0 the computation of \: 
and A; is materially simplified (13) since the Euler Equation 
[23] admits a first integral 


JET PROPULSION 


M I M 8 
A 
A 
y 
ng 
x 


ts A 
and 
such 
or a 
tra- 
rust 
the 


own 
own 
ting 

the 
chi- 


ree 


Maximum time solution 


(C=-006) 
Minimum time solution 
(C = 068) pe 
4 08 
02 0. 06 
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C; being an integration constant. The calculation of A, and 
\; is ultimately reconduced (13) to a process of simple quad- 
ratures. 

2 Another method of solution consists in attacking the 
fixed end-points problem for several values of the final Mach 
numbers M;; calculating the time 7; necessary to reach each 
final Mlach number M;; plotting the resulting function 7;(M,); 
and finding, by a graphical procedure, the trajectory yielding 
(t;)mine This method bypasses the natural condition result- 
ing into Equation [33], but it is slightly less precise in the 
determination of the optimum trajectory. 

3 The above two methods (which yield consistent results) 
have been developed to a considerable extent in (13). As 
shown by a typical example, the computational effort in- 
volved is approximately the same with both techniques. 
Method 2 is closer perhaps to the immediate intuition of the 
engineer and offers a clear picture of the effect which a 


deviation from the optimum flight program has on the final 


time 7,. 


= 
Numerical Example 


A sample calculation is now carried out for the following 
set of conditions! 


M; =0 nn, =0 =0.9 
M, = 10 “+ = 2 €; = 5.624 
=? 

m,/m; = 0.62 y =0.3 Mmax = © 


[64] 


For the sake of brevity only method 2 is employed, even 
though identical results have been obtained with method 1 
in (13). 

As a first step the fixed end-points problem is solved for 
several final Mach numbers M;. The family of brachisto- 
chronic paths JABF satisfying the above set of conditions 
and leading to different fixed final Mach numbers M, is com- 
puted (Figs. 9 to 12). 

Each brachistochronic path JA BF can be described either in 
terms of the final Mach number M, or in terms of the constant 
C associated with the singular sub-are AB. Notice that each 
trajectory JABF includes: an initial pulse-burning sub-are 
IA, a central sub-are AB flown with variable thrust and a 
- sub-are BF. With regard to the latter, the analysis 
shows: 

1 For M, = 7.7 (C = 0.36) the final point F is reached 
with regulated thrust (Figs. 9 to 12). 


"| Notice that to simplify the analysis the ideal case of an engine 
capable of a pulse burning [umax = ©] is considered. 
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Fig. 11 Family of brachistochronic paths leading to different fixed 
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Final point reached 
by coasting 


Solution of the free 
boundary value problem 
(Mp = 5.91, C =6.25) 


02 


Fig. 13 Time necessary to reach the final point along a brachis- 
tochronic path as a function of the final Mach number 


08 
2 For M,>7.7 (C < 0.36) the final point F is reached with 
a pulse-burning sub-are BF (Figs. 9 and 10). B 
3 For M, < 7.7 (C > 0.36) the final point F is reached 
with a coasting sub-are BF (Figs. 11 and 12). 
Once all possible brachistochronic paths 7A BF are drawn in 04 
the (», M) plane the time necessary to reach the final point Fy) 
is computed from the line integra, (Fig. 13) © 


d 
ixed M @ A 


of 


Clearly the special path minimizing the time for free final 


Mach number is to be consistent with the necessary condition Fig. 14 Particular brachistochronic trajectory for free final Mach 
number 


As it is evident from the integral [65] the condition [66] 
cannot be met by any of the paths JABF ending with a 
pulse-burning (Fig. 9). This equation, on the other hand, 
can be met by a path JABF ending with a coasting flight (Fig. 
— ll). In this connection, Fig. 13 shows that the trajectory 

3 solving the free boundary value problem is the path JABF . | 
corresponding to C = 6.25, M, = 5.91 and terminating with bs 
a coasting flight. Such a path is indicated in Figs. 14 and <“/y : 
15 


Remarks 


For the particular examples previously discussed the 12 
extremal path was shown to be composed of sub-ares of three 
different kinds. There are cases, however, where the com- 
position of the extremal path is quite different from the one 
indicated in Figs. 6 to 15. For instance, the variable thrust 
sub-are may not be present in the solution; or the extremal 
are may reduce to a single sub-are flown with maximum 
engine output. The particular extremal are associated with 
a given problem strongly depends on the characteristics of 
the airframe and of the engine and, in addition, on the 
boundary conditions of the problem. 


Conclusions 

a The brachistochronic burning program for a rocket powered 

— missile moving along a rectilinear, near vertical trajectory 

_ imbedded in an isothermal medium is analyzed. Under O 

suitable assumptions for the drag function, solutions in a 

rod closed form are obtained for the relationships betwe een alti- Fig. 15 Particular brachistochronic trajectory for free final Mach 
tude, time, mass, and Mach number. number 
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The boundary value problem is investigated for the case 

where the final Mach number is prescribed and the case 
where it is free of choice. 

For the fixed end-points problem, a numerical example 
- shows that a given set of end conditions may be met by two 

different extremizing solutions: The brachistochronic solu- 

tion and one which maximizes the time. 

In regard to the free boundary value problem, two pro- 
cedures are developed. One is based on the transversality 
condition; the other one is conceived as a sequence of re- 
peated solutions of the boundary value problem for several 
fixed final Mach numbers. The amount of computational 
work necessary to predict the optimum path is roughly the 
same in both methods. 
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Technical Notes 


A Method for Calculating Impact 
Points of Ballistic Rockets: 


Convenient Representations 


F. SINGER! and R. C. WENTWORTH? 
__ University of Maryland, College Park, Md. 


Nomenclature 


a = semimajor axis of orbit ellipse (a = [Ra + Rp]/2), m 

A; = fraction of the area of the ellipse swept out by r in going 
from the perigee to R,(A, =[cos™ [(Yz? — 1)/e] — 
e(1 — [1 — 

a = angle in the inertial frame between north and the pro- 


jection of the launching direction in the inertial frame 
on the horizontal plane measured clockwise as seen 
from above, deg or rad 


Received May 14, 1958. =, je 
1 Associate Professor of Physics. Member ARS. 
Graduate Research Assistant. 4 


L,), deg or rad 


169, 


= angle in the Earth’s frame between north and the projec- 
tion of the launching direction on the horizontal plane 
measured clockwise as seen from above, deg or rad 

= difference in longitude between Z; and L; (AL = L; — 


‘ = change in longitude in the Earth’s frame, deg or rad 


= total angle r has turned through from r = Ry, on one 
side to r = R, on the other side of the ellipse in the 
flight of the rocket, deg or rad 
= eccentricity of the orbit ellipse (e = (1 — Yy,? cos? 6 
= Newton’s constant (G = 6.670 X 107! 
nt:-m?- kg~?) 
h = peak altitude reached by the rocket above R,, m 
6 = angle between the launching direction and the horizontal 
in the inertial frame, deg or rad 
6’ = angle between the launching direction and the horizontal 


in the Earth’s frame, deg or rad 
L; = longitude of the point of impact at a distance Rz, 
the center of the Earth, deg or rad 


L, = longitude of the Earth where rocket was launched, d 
rad 
A, = latitude of the point of impact at a distance R; from 


center of the Earth, deg or rad 


within two months of the date of receipt. 


Eprror’s Note: The Technical Notes and Technical Comments sections of JET PROPULSION are 0 
developments or offering comments on papers previously published. Such manuscripts are published without editorial review, usually 
Requirements as to style are the same as for regular contributions (see masthead page). 
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\, = latitude of the Earth where rocket was launched, deg or 
rad 

\,/ = impact latitude in the Earth’s frame, deg or rad 

Mz = mass of the Earth (Mg = 5.983 X 10*4 kg) 

P; = impact point 

P, = launching point 

P;' = impact point in the Earth’s frame 

¢, = the value of the polar angle (i.e., geocentric angle) of the 
orbit measured from the perigee to the launching 
point, deg or rad 

Ry = apogee distance, m 

R, = distance of launching point from center of Earth, m 

Rp = perigee distance, m 


i = total time of flight of the rocket above Fz, sec 

t = period of a circular orbit of radius Ry (t- = 2rR,/V-), sec 

T, ~= period of the rocket in the complete ellipse (7. = 27 
(GM —/?a3!2), sec 

V. = velocity of an object in a circular orbit of radius R, 
(Ve = m- 

VY, = launching velocity magnitude in the inertial frame, 
m-sec™} 

VY,’ = launching velocity magnitude in the Earth’s frame, 
m-sec7! 

V,’ = launching velocity vector of the rocket in a reference 
frame fixed with respect to the Earth, m- sec! 

» = angular rotation of the earth (o = 7.272 X 10-5), rad- 
sec! 

Y, = dimensionless launching velocity in the inertial frame 


Y, = Vz/V-) 
(Yr L/ 


Introduction 


\IETHOD for determining impact points of ballistic 

rockets on a rotating Earth has been developed by 
Singer and Wentworth (1).* In the present note, a dimen- 
sionless representation is given which puts the results of (1) 
ina more convenient form. Also, (1) is amended to include 
nonequatorial launchings. 

The method used in (1) to calculate the impact point of a 
ballistic rocket considers the elliptic trajectory of the coast- 
ing rocket in an inertial reference frame fixed with respect 
to the stars rather than using a coordinate system which 
rotates with the Earth. The rotation of the Earth is then 
taken into account by adding a suitable horizontal com- 
ponent to the launching velocity of the rocket, and by con- 
sidering the angle turned through by the earth during the 
time of flight of the rocket. 

The launching parameters are the velocity V,’ of the rocket 
just after final burnout (the subsequent flight of the rocket 
being therefore unpowered) and the latitude \,, longitude 
L,, and radius from the center of the Earth R, where final 
burnout occurs. For simplicity, the impact point is taken 
to be the same distance R, from the center of the Earth, and 
it is also assumed that the final burnout takes place above the 
atmosphere, so air friction can be neglected. 

Specifically, the launching velocity as seen from the Earth’s 
frame V,’ is designated by its magnitude V,’, the angle it 
makes with the horizontal 6’, and the angle its projection 
makes with north, measured clockwise a’. Similar quantities 
as seen from the nonrotating inertial frame, V,, 6, and a, are 
calculated from these, and V, and 6 directly determine the 
orbit ellipse. The orbit ellipse in turn determines the angular 
range Ad, the peak altitude h, and the time of flight ¢ of the 
coasting rocket. Knowledge of these quantities coupled 
with the launching latitude \, and longitude Z, then enables 
one finally to determine, with the aid of some spherical trig- 
onometry, the point of impact. 

However, since the publication of (1), it has been found 
possible to put the important results of that paper into con- 
venient dimensionless form. Thus, the angular range Ad, 
the dimensionless time of flight of the rocket ¢/f,, and the 
dimensionless altitude h/R, are expressed as functions of the 


’ Numbers in parentheses indicate References at end of paper. 
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Fig. 1 The elliptic orbit of a ballistic rocket showing notation 
used : 
. 


Y= Vi/ Ve 
é =(I-ye cof [2 - y2])”* 
ag = 2 sin” (ye sin20/ Ze) 


h/R,= 


N 
/ 
a¢ 
dr 
y 
L- aL*L,-L, 


Fig. 2. The angular trajectory A¢ of a ballistic rocket on the 
surface of the Earth showing notation used 


dimensionless launching velocity Y, and launching angle @, 
in the first part of this note. 

Also, the authors of (1) received a communication from 
Professor E. E. Zukoski* who pointed out that Equations 
[19 and 20] of that paper were in error for nonequatorial 
launchings and indicated (2) as containing the correct solu- 
tion. Therefore, the second part of this note is devoted to a 
derivation of the correct forms of those equations which give 
the impact point as a function of the launching point and the 
launching parameters. 


4 Assistant Professor of Jet Propulsion, California Institute of 
Technology, Pasadena, Calif. 
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Dimensionless Representation 


It is most convenient to work with dimensionless quantities 
in the following expressions, and with that in mind we will 
make use of the so-called circular velocity V,, that is, the 
velocity required for an object to maintain a circular orbit at 
a distance R, from the center of the Earth. Then the di- 
mensionless launching velocity Y, is the ratio of the actual 
launching velocity V; at the launching radius R, to the cir- 
cular velocity at that radius. 


Y, = 


Ve 
Using the definition of V,, Equations [6] of (1) become 


= 2V,? — V,? 
= 2R,V.? 
= R,2V;2 cos? 6 
Substituting [2] into Equations [8] of (1), we have 


Ra = + [1 — (2V22 — cos? 
(2V2 — Vi! 
Rp = Ri — [1 (2 — Y_,?)(¥ cos? (2 — 


Then making use of Equations [1, 2, and 3], Equations /10, 
11, 15 and 17] of (1) reduce, with a few minor changes in 
notation, to 


= (1 — cos? — 
= cos~! cos? — 1]/e) 
= R1/(2 — Y;z?) 
27Az, = cos [(Y,? — 1)/e] — e{1 — (1 — 


Using [4], Equation [18] of (1) becomes 
Aé = 2 ein [sin [5] 

The dimensionless peak altitude h/R, is 
h/Ry = (Ra — = — 14+ (1 + €)/(2 — Yz?). . [6] 


Finally, defining the circular period ¢, of an object in a cir- 
cular orbit of radius R;, about the Earth by 


Equation [12] of (1) may be written 


so that we have for the dimensionless time of flight, making 
use of Equation [13] of (1) 


t/te = [1 — = [1 — 2Az]/(2 — 


Equations [5, 6 and 9] are plotted as functions of Y; with @ 
as a parameter in Figs. 3,4 and 5. These graphs then easily 
permit the determination of the angular range (in radians) 
Ad, the altitude h and the time of flight ¢ as functions of the 
launching velocity V; and angle 0. ieee 


Impact Point 


Suppose that the launching velocity vector as seen from 
the earth V,’ is defined by its magnitude V,’, the angle it 
makes with the horizontal plane 6’ and the angle measured 
clockwise that its projection on the horizontal plane makes 
with north a’ (see Fig. 2). 

To determine the corresponding quantities in a fixed refer- 
ence frame, one must add to V;’ the velocity to the east due 
to the rotation of the Earth wR, cos \,. Then a considera- 
tion of the geometry of the situation yields 


a = tan~! [(Vz' sin a’ cos 6’ + wR, cos cos 6’ cos a’) 
= tan~! [tan a’ + wR, cos A,/Vz' cos cos 
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[2] 


-€ + cos" [(v2-1)/e] 


Also 
tan = sin n = cos cos a’/cos a... [11] 
or 
+ 6 = tan™ [tan 0’ cos a/cos a’]........... [12] 
and 
V, = sin 0’/sin @.......... 


In the special case of vertical launching in the Earth’s frame 
a = + 90 deg a 
6 = [V,'/wR, cos AL] 
6’ = + 90 deg 
= (Vz, + w?R,? 


It has been shown how Ag and t may be determined with 
the aid of [12 and 13] and Figs. 3 and 5. As can be seen 
from Fig. 2, the dependence of \; and L; on Ad, a, and d; is a 
fairly complicated problem in spherical trigonometry. How- 
ever, it has been solved in (2), from which it follows, with 
appropriate changes in notation, that 


= sin=! (cos Ag sin A; + sin Ag cos cos Xz) 
L; — = AL = cot~(—sin A, cot + cos Az cot Qa) 
. [15] 


In [15], the change in longitude AZ between the launching 
longitude L, and the impact longitude L, is positive if the im- 
pact point P; is to the east of the launching point P;. Also, 
hy and \, are positive if they are in the northern hemisphere 
and negative if they are in the southern hemisphere. 

Vinally, the impact point P,’ in the Earth’s reference frame 
is given by 
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Satellite Orbit Perturbations Under a 
Continuous Radial Thrust of 


Small Magnitude 
DOBROWOLSKE 
‘Flight Propulsion Laboratory, General Electric Co., 
Cincinnati, Ohio 
5 
Nomenclature 
4 
F = thrust in radial direction ’ 5 
r = distance from the center of Earth =" 
ro = radius of the initial circular orbit 7 + 
R = radius of Earth " 
F =~ 
~ mg \R 
¢@ = angular displacement 
HE problem of take-off from satellite orbits has been i in- 
vestigated by Tsien (1),? Ehricke (2), Peterson (3) and 
others. In what follows, those investigations are supple- 


mented by discussing the condition of continuous radial 
thrust when 8 \ < 1. Tsien’s paper was based on the con- 
dition of 8\A> 1. 
b- Initially, the spaceship is in a circular orbit of radius ro. 
At time ¢ = 0, a radial thrust in outward direction and such 
that F = const, m = const begins to act. From equations of 
motion 


mr — ro?) = 


¢ is eliminated to obtain 


gk? a 
m 


Integration and substitution of initial condition (7),.,, = 0 
gives 


2F Rt7 1 R? 1 
ft = -1 + -1) 
m To To \r/ro 


where 


and is the period in the original circular orbit 
mg \lR 
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the following results = , Also Equation [8] under similar approximations will give 
4 Integrating Equation [3] It is seen that a continuous radial thrust, such that 8 \ <1, 
Aa ; will cause only a temporary departure from the satellite orbit, 
= oP = f wink. [3a] The vehicle is bound to return, and the maximun radia! dis- 
ts ; V 2 Af(s) placement never surpasses in length the radius of the original 
where 


The graph of f(s) against s (Fig. 1) shows that s oscillates 
between unity and circular orbit. This is seen from Fig. 2. Fig. 3 shows the 
perturbed orbit. 


Vian 


4X 


formin th two transformations Tsien, H. S8., ‘““‘Take-Off from Satellite Orbit,’ Jer Pro» 
SION, vol. 23, July-August 1953, pp. 233-236. 


s=1+2? [4] 2 Ehricke, Krafft A., “Take-Off from Satellite Orbits,’’ Jer 
Proputsion vol. 23, Nov.-Dec. 1953, pp. 372-374. 
The integral is evaluated _ A : 3 Peterson, J. B., ““Note on ‘Take-Off From Satellite Orbi:’, ” 
ee _ Jer Proputsion, vol. 27, Dec. 1957, pp. 1263-1264. 
where F and £ are elliptic integrals of the first and second f(s) . 
kind, respectively. apy b 
From [2] ‘ 
dP 
0 
fe 
4 i 
+ 


Fig.1 The interval within which Equation [3a] can be integrated 


— s + 1)(s — 1) 


oh 
Integration of the above between s = 1 and 8max = 1 + 6? a 
gives 

=4/= b II (b2, 7/2). 

(1 + bx?) — — 2?) 
= =? 

where II is the elliptic integral of the third kind. This last ome 
expression can be shown to be equal to | 


»b 
= pe + — 2/2) (V1 — 8) — (V1 — 54,0) + F (0%, 2/2)F (V1 — 04, 


where cot 6 = b. 

The above gives exact expressions for the perturbed orbit. 
However, some approximations can be executed. From 
Equation [6], by expanding in series 


Then if 8 \ < 1, and hence, 6? ~ 2 X, will transform the pre- 
into 
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Fig. 2 The variation of b?, the measure of the maximum 
radial displacement from the original circular orbit, against 
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Fig. 3 Orbit perturbations under continuous radial thrust 


Flat Plate Laminar Skin Friction and 
Heat Transfer in the Free Molecule to 
Continuum Flow Regimes 


HAROLD MIRELS! 


National Advisory Committee for Aeronautics, 
Cleveland, Ohio 


EVERAL authors have solved Rayleigh’s problem (i.e., 

impulsive start of an infinite plate) using continuum 
equations with slip boundary conditions and have applied 
the results to estimate skin friction (1, 2)? and heat transfer 
(3) on a semi-infinite flat plate in the free molecule to con- 
tinuum flow regimes. The transformation from time, in 
Rayleigh’s problem, to distance downstream of the leading 
edge, in the semi-infinite flat plate problem, according to the 
relation u¢ = Kz, introduces the arbitrary constant K. 
This constant was chosen, in (2), such that the expression for 
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wall shear approached the correct value in the limit of con- 
tinuum flow. The resulting expression overestimated the 
wall shear in the free molecule regime by a factor of 2. In (1 
and 3) the constant K was chosen in such a way that the wall 
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Equation [2b] assumes 7’, is constant. 


f(x) 


shear and heat transfer did not agree with the correct limiting 
values for either the free molecule or the continuum regimes. 

More recently, (4) has studied Rayleigh’s problem at low 
Mach numbers from the point of view of kinetic theory. An 
approximate expression for wall shear was obtained which is 
valid in the entire range from free molecule to continuum 
flow. Surprisingly, this expression is identical—except for 
the numerical constants—with the expression for wall shear 
which results from the solution of Rayleigh’s problem using 
continuum equations with slip boundary conditions. This 
suggests that the latter solution can be used as a model to 
provide the functional dependence of wall shear and heat 
transfer on Mach number and Reynolds number in the oe 
molecule to continuum regimes. In particular, if time is re- 
placed by distance and all numerical constants are arbitrarily 
modified so as to provide agreement in the limits of free 
molecule and continuum flow, then expressions for semi- 
infinite flat plate shear and heat transfer can be obtained 
which should give reasonable estimates of these quantities _ 
over the entire range from free molecule to continuum flow. © 
Such an approach is outlined herein. { 

Consider Rayleigh’s problem for a fluid with constant 
properties. The momentum and energy equations, with slip — 
boundary conditions are, respectively (1,3) 
Ou 


ot oy? 


Ou 
u(t,0) = 1 (=) 


u(t, ©) = u(0, y) = ue 


oT vor 
ol a Oy? 
oT 
T(t, 0) = Ty +m (27) 
Oy J w 


T(t, ©) = T(0,y) = Te 


where / and m are constants defining the velocity slip and 
temperatufe jump at the wall. Equation [1b] neglects dis- 
sipation. The resulting nondimensonal expressions for wall 
shear and heat transfer are (1,2,3) 


(« 
OY] w 


2u = 


» — Tw) 


If t is replaced by 2, 


Equations [2] have the form 


2 
= aft (b x) + (b 
> VT 
(b Vr)< 1. . [3b] 


(bV/x) > 1..[3e] 


where a and b are independent of x. Considering the limits 
of free molecule flow (b+/.x > 0) and continuum flow (6./z > 
«), aand can be expressed as (from Equations [3b, 3c ]) 
a = f(0) 
Vr lim [Vz f(z)] 


give 
al di 7 
Te) : 
L- 
05 
JET 
fi 
FAR Ste 
j 
~ 
j 
o= 
\_ 
Dp 
| | 
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Thus a and b can be evaluated from known wall shear and 
heat transfer results in the free molecule and continuum re- 
gimes. 


since they were matched to continuum solutions which include 
these effects. The wall temperature is assumed constant in 
these equations. The overall friction and heat transfer co- 
efficients for a plate of length L can be obtained by noting 


2 
f(x)dx = Ee (b\/L) — 1 + erfe (b\/L) exp | [9] 


The present approach results in relatively simple analy: ical 


Identifying f (2) with Cy, we get expressions for wall shear and heat transfer on a semi-inf nite 


a flat plate. A major deficiency is that the low-speed Ray!) igh 

@=3(Cyeur _ ee ae [5a] problem model used here does not take into account boun: ary 

e Vary d layer induced modifications of the free stream. For exan ple, 

consider Equation [8c] to represent local shear in ne 

ayy fe continuum flow. Maslen (7) has shown that the first de- 

ee (Crono M. Nry [5b] parture from flat plate continuum theory is associated with a 

VT lim Iva zCy) self-induced pressure gradient which introduces a ive 
= 0.664 <7 term of order (1/b+/z) in the bracket of Equation [8c]. ‘he 


present method therefore tends to underestimate the vall 
shear, at least in the border regime between continuum :ind 
slip flow. With regard to heat transfer, Maslen (8) has ilso 
shown that the departure from continuum first appears ::s a 
term of order (1/b+/z)? in Equation [3c] so that the pre. ent 
model has the correct order of magnitude with regard to the 
heat transfer perturbation in the border region between «on- 
tinuum and slip flow. 


Equation [5a] is the free molecule flow result for diffuse re- 
flection (5). The denominator of Equation [5b] employs the 
Chapman-Rubesin (6) continuum result for wall shear. The 
latter assumes viscosity proportional to temperature accord- 
ing to the relation (u/ue.) = C(T/T.). Other choices for the 
free molecule and continuum shear could have been used. 
Substituting Equations [5] into Equation [8a] gives (for y = 
1.4 and Re = u.2/vC) 
u exp { 0.573 —— M7 Pn |) 1 Schaaf, S. A., “A Note on the Flat Plate Drag Coefficie t 
: “4 Rep. No. HE-150-66, Inst. Eng. Res., Univ. Calif., Feb. 20, 1950. 
2 Mirels, Harold, “Estimate of Slip Effect on Compressible 
Laminar-Boundary-Lay er Skin Friction,’ NACA TN _ 2609, 


M,C; = 0.674 erfe ( 


Identifying f(x) with qu, as defined in Equation [2b], we get 


Equation [7a] is the free molecule result for an accommoda- 
tion coefficient of one (5) while the denominator of Equation 
[7b] employs the Chapman-Rubesin heat transfer result for 
= 0.72. If, in addition, y = 1.4, the expression for local 
heat transfer becomes 


M 
ta 


1 +0.350M,2 — 


7’. 
1 + 0.350M2 — — 


= 0.193[erfe z] exp (2?)...... [8] 


where 


2 = 0.265 


In order for the expansions of Equations [3] to be valid, the 
coefficient of V Re/M., in Equation [8], must be positive. 
Thus it is required that 


4 
: T. > 1+ 0.350 M,.? = F. M. recovery temperature 


or 


< 1+ 0.170 M,? = continuum recovery temperature 


For Equation [8] to give reasonable results, the ratio 
{1 + 0.350 M,? — (Tw/T.)]/[1 + 0.170 — (Tw/T,.)] 


should be of order one. 
Equations [6 and 8], which are the main results of the present 
study, include effects of variable properties and dissipation 
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Relaxation Lags in Gas Flow 


HENRY AROESTE! and MORTON RUDIN? 
Lockheed Missile Systems Division, Palo Alto, Calif. 


T IS the purpose of this note to present a brief discussion 
of methods presently used to estimate lags in flow with 
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consequent determination of the approach to equilibrium. 
One of the authors has recently developed a more precise 
treatment in which the detailed couplings of the flow equations 
are considered (1).* It is interesting to compare his results 
with those from more approximate considerations (2, 3). 

In the case of vibrational temperature lag it is possible to 
separate a comparison factor J, which one may multiply with 
the more approximate result to achieve the more exact cri- 
terion. Following this procedure the vibrational lag in tem- 
perature for a diatomic molecule may be written as 


where 7 is the usual vibrational relaxation time, and DT/Dt 
represents the usual cooling rate. An outline of a derivation 
will be given later to show that 


7(7 — 5K) 
—1)(7 + V2/RT) 


Hee, V is the flow velocity, R is the gas constant divided by 
the molecular weight and K is given as 
_ 7/2 + (0/T)*/{exp (8/7) — 1)? 


wer 


where is the “vibrational temperature.” 

The factor J has been computed under various conditions 
for the molecule Oz, and the results are displayed in Fig. 1. 
From Equations [2 and 3] it may be observed that in the 
limit of infinite temperature J approaches 1. This yields 
the approximate formula‘ 


This formula when used at finite temperatures is a conserva- 
tive estimate of the temperature lag; i.e., it gives a consider- 
ably larger value of the lag than when computations are 
made with the appropriate values of J for the temperature 
and velocity of interest. A temperature lag computed from 
[4] involves the following assumptions: (a) Complete vibra- 
tional partition (vibrational energy = RT) is attained, a 
condition only true at infinite temperature; and (b) there is a 
lag only in temperature, whereas actually there is a coupling 
in order to conserve energy such that there is also a lag in 
velocity. When the latter effect is considered, one obtains 
a velocity dependence as shown in Equation [2]. This de- 
pendence, which is small at the conditions for which the com- 
putations were made, vanishes at infinite temperature. 

Equations [1 to 3] may be derived with the use of Equations 
[2C and 16] in (1), given below in modified form as Equa- 
tions [5 to 7] 


a. AT/T = [5] 
where 
E, = (V?/RT) (1 — K)/((V2/RT) — K)........ [6] 
and 
(1/A)dA/dl = AT/JTVTE,............. [7] 


where / is the length along a streamline. Equation [1] is 
seen to follow directly’ from Equations [5 to 7] with the 
auxiliary use of the relationship 


[8] 


Numbers in parentheses indicate References at end of paper. 
The formula used in (3), which neglected 2/7, results in an 
even more conservative estimate, as if all the energy change 
during flow came from the vibrational mode. 
° The above derivation fails at the choking point because dA = 
0; nevertheless, the same Equations [1 to 3] may be shown to 
follow from the equations valid at choking (1). 
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Fig. 1 J, the comparison factor, as a function of velocity in 
ft/sec for various absolute temperatures 


which is true under steady state conditions. 

It should be noted that according to (1) AT gives the value 
that the temperature lag approaches under the condition that 
d(AT/T)/dl = 0. A better estimate for AT may be ob- 
tained from a determination of d(AT/T)/dl by the methods 
described in (1). It is therefore significant to state that if a 
careful measurement of rates, whether internal or chemical, 
is considered through the measurement of lags, one requires 
either a measurement which will yield the derivative along 
the streamline or measurements at two points so that the 
derivative may be computed. 

In the case of chemical lag it has not been possible to sepa- 
rate a similar comparison factor, although it is possible to re- 
duce the results in (1) to the form, AT = zDT/Dt, which is 
the same form of the result in (2). Here z is the so-called 
reaction time. It would be likewise interesting to compare the 
results of calculations of z by the two different methods for a 
simple case as dissociation. In computing the reaction time 
by the method of (2) one neglects the following three effects, 
which are considered in (1). First, and probably most im- 
portant, is that the temperature lag as defined in (2) would 
result if all the energy change during flow came from dissoci- 
ation. This is analogous to the estimate for vibration given 
in (3).6 Second, an analogous velocity effect would be pre- 
sumed to be of the same order of magnitude as in the vibra- 
tional comparison given here. Third, there would be an 
added effect of density variation which does not appear in 


the case of vibration and is also probably small. 
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Optical tracking recording system (2,- 
836,894). H.A. Wagner, Thousand Oaks, 
Calif., assignor to H. A. Wagner Co. 

Tracking sighting means responsive to 
the movements of a control stick. Signals 
generated by the rate of motion of the 
control stick are recorded. 

Fuel control for a gas turbine power plant 
(2,836,957). S. S. Fox, W. Hartford, 
Conn., assignor to United Aircraft Corp. 

Metering valve in the conduit upstream 
of a chamber through which the fuel flows, 
feeds maximum fuel into the chamber as a 
function of turbine speed, compressor air 
inlet temperature and compressor dis- 
charge pressure. 

Jet power plant with unobstructed rotating 
combustion chamber (2,836,958). J. A. 
Ward III, Duluth, Minn. 

Rotary jet engine with front and rear 
supports having axially aligned bearings 
in which hubs rotate. Gases flow through 
the engine in a straight path, and a sta- 
tionary fuel supply means extends into 
the combustion chamber through one of 
the hubs. 

Gas turbine engine supporting frame 
(2,836,959). C. J. McDowell and O. V. 


Montieth, Indianapolis, Ind., assignors to 


General Motors Corp. 

Rotary engine with a cylindrical support 

secured to a torque ring, and extending 
rearward. A central shaft has a rear bear- 
ing connected by rear struts to the cylin- 
drical support. 
Wing mounted jet propulsion system 
(2,838,256). I. M. Davidson and N. A. 
Dimmock, Woking, England, assignors 
to Power Jets (Research and Develop- 
ment), Ltd. 

Hollow wing with an air inlet extending 

along the leading edge. A hollow main 
spar is formed with inlet and outlet 
apertures. Jet engines are mounted with- 
in the wing supplied by air from air supply 
passages and discharging in narrow ex- 
haust slots along the wing trailing edge. 
Jet sustained aircraft with enclosed 
compressor rotor (2,838,257). M. Wi- 
bault, New York, N. Y., assignor to 
Vibrane Corp. 

Centrifugal blower in the body of an 
aircraft. The downstream end’ of the 


outlet passage through which air is dis- . 


charged from the blower discharges into 
the ambient atmosphere in which the air- 
craft flies, and different parts of the down- 
stream end are at different locations with 
respect to the c.g. A control varies the 
amount of air discharged at the different 
locations. 


Device for demonstrating mathematical 
probability (2,838,851). R. Lusser, 8. 
Pasadena, Calif. 

Group of balls each movable with respect 
to all remaining balls, and composed of 
electrically conductive material, except 
that a predetermined area (constituting a 
minor fraction of the surface of some of 
the balls) is nonconducting. An elec- 
trical circuit in each ball is completed by 
contact with the conducting surface of 


Epritor’s Note: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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each ball. Movement results in a random 
rearrangement of the balls, breaking the 
circuit when contact with any one or more 
balls lies within a nonconducting surface 
area. 

Means for the deflection of fluid jets 
(2,838,909). H. L. P. Meulien, Courbe- 
voie, France, assignor to SNECMA. 

Pair of movable members mounted for 
angular displacement in the exhaust 
nozzle. The members are streamlined 
when closed, allowing undisturbed flow 
of the jet discharge. When opened, they 
expose two intercepting surfaces to deflect 
the flow. 

Sun follower (2,839,689). G. M. Trinite 
Jr., Baltimore, Md., assignor to Aircraft 
Armaments, Inc. 

Body containing coarse-eye photo- 
electric cells responsive to light rays over 
a wide sector or space to produce an out- 
put voltage. A plurality of fine-eye 


cells responsive to rays within a relatively 
narrow sector of space within the wide 
sector, to produce an output voltage. 
The fine-eye cells alone provide output 
while the sun follower is closely tracking 
the sun. 


Aircraft design (180,297). E. Schmued 
and W. E. Gasich, Los Angeles, Calif. as- 
signors to Northrup Aircraft, Inc. 

Delta wing jet airplane with air intake 
scoop below fuselage, aft of cockpit. The 
tail is at the same level as the wing. 
Variable area jet nozzles for gas turbine 
units (2,839,890). F. W. W. Morley, 
Castle Conington, England, assignor to 
Rolls-Royce, Ltd. 

Jet pipe with an air inlet leading to the 
space between inner and outer skin mem- 
bers. An outlet in the nozzle induces the 
flow of cooling air through the air pas- 
sage. Nozzles may be adjusted to vary the 
orifice area. 

Thrust reversers for jet engines (2,839,- 
891). G.T. Drakeley, Seattle, Wash., as- 
signor to Boeing Airplane Co. 

Deflector flaps mounted for movement 
between an inoperative position along the 
tailpipe and an operative position to the 
rear of (and directed outwardly from) the 
tailpipe axis. In operative position, gases 
are directed forward as they follow the in- 
terior surface of the flaps. 
Gas Turbine cycle (2,839,892). H. Rosen- 
dahl, Yonkers, N. Y. 

Process of making valuable products of 
incomplete oxidation of a fluid-fuel by 
employing secondary fuel as a coolant and 


George F. McLaughlin, Contributor 


utilizing the turbine exhaust gases in 
chemical reactions. Turbine gases are 
passed into a chamber in the presence of a 
catalyst, using the sensible heat of the 
turbine exhaust for initiating catal\tic 
reaction. 

Multiple vortex tube generator cooling 
unit (2,839,898). F. H. reen, ‘os 
Angeles, Calif., assignor to The Garrett 
Corp. 

System for cooling and pressurizing an 

instrumentality in an aircraft having a jet 
engine compressor and a ram air dict. 
Flow of high pressure air contro'led 
through vortex tubes is limited when sonic 
velocity of air is reached. 
Regenerative vortex cooling systems 2,- 
839,900). Green, Palos Verdes 
Estates, Calif., assignor to The Garrett 
Corp. 

System for supplying cooling air to an 
aircraft. Means for passing ram air in 
heat exchange relation to first and second 
vortex tubes. 

Position indicator (2,839,921). A. Wood, 
Albuquerque, N. M. assignor to Engincer- 
ing Development Co., Inc. 

Rack representing the equator, in en- 
gagement with arms on a pivot represent- 
ing the elevated pole. A track guide rep- 
resents departure and destination posi- 
tions, and a movable guide corresponds to 
estimated ground speed. A gyro compass 
having latitude correction is connected to 
a latitude gear engaged to the rack and 
arms. 

Shaped charges (2,839,997). J. H. Church 
and G. J. Kessenich, Madison, Wisc., as- 
signors to the U.S. Army. 

Pair of identical six-pointed star-shaped 
sheets of material having uniformly 
spaced arms. Sheets are connected paral- 
lel axially, and a strip of material extends 
along the edges of the sheets to form a 
star-shaped container which is filled with 
an explosive charge. 

Air intake for aircraft power plants (2,- 
840,322). A.A. Griffith, Derby, England, 
assignor to Rolls-Royce, Ltd. 

Center body dividing an air intake into 
two passages. The body has a wedge 
shaped so that at the design Mach number 
the shock wave formations extend outside 
or onto the leading edge lips. Intake 
passages may be varied in effective area. 
Hinged airfoil with jet nozzle for aircraft 
flight control (2,840,323). G. K. Hood 
Jr. and J. A. O’Malley Jr., N. Tonawanda 
N. Y., assignors to Bell Aircraft Corp. 

Attitude control consisting of a pair of 
telescoping cylindrical members rotatable 
about the longitudinal axis. One member 
is mounted on a fixed airfoil and the other 
on a movable airfoil. A duct conveys 
pressurized gas into the inner member so 
that upon swinging the movable airfoil 
away from its neutral position, parts align 
to release the gas into the atmosphere, 
rotating the airplane about a flight control 
axis. 

Aircraft structure (2,840,325). A. A. 
xrifith, Derby, England, assignor to 
Rolls-Royce, Ltd. 

Pilot’s windshield with an external sur- 
face lying in the surface of a cone inter- 
secting the external skin, and disposed 
with its apex forward. At supersonic 
speed a shock wave formation is generated 
which is of conical surface form with its 
apex coincident with the apex of the wind- 
shield. 
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Emulsion Types: 7ype / has the highest order of 
light sensitivity provided by modern emulsion- 
ae making methods, medium contrast (gamma 

somewhat greater than 1), and resolving power 


n TE = of about 60 lines per millimeter for a test pattern 

€ Contrast we’ of 20:1 contrast. With this as a basis, the chart 
a a Resolving Power = shows how Types I, 103, II, III, 1V, and V differ. 
3 In addition there are: 
Type 649—Extreme resolving POWEE- 

a—Simil: Il and 103 

ypes Ila and 103a—Similar to an 
; Tuvecessssesescccsooaverserercestessss respectively, but much higher in sensitivity for 
t low-intensity, long-time exposure. 
n 
t 
This chart shows the approximate relationship between the principal emulsion 
d types with respect to light sensitivity, inherent contrast, and resolving power. 
= 


z \V) 
1030 
IV) 
1 IV Vv) 1030 
t 
‘iL. N ZZ (1#, IV, Vv) 
F (1, Ha, IV, V, 103, 1030, 649) 
) E 11, Ha, IV, 103, 103a) 
Ai, Wl, IV, Vv, 103, 103a, 649) 
B (1, IV, Vv, 103, 103) 
| We ho e D (1, 11, Ha, IV, V, 103, 103a) 
p T (1, HE, IV, V, 103, 103a) 
ZA(\, IV, V, 103, 1030, 649) 
H Wt, IV, V, 103, 103a) 
these charts | « 
J (1, 11, Ha, IV, V, 103, 1030) 
. (1, 11, Hla, IV, 103, 1030, 649) 
im ress 300 400 500 600 700 800 900 1000 200 
= TOTAL SENSITIVITY SPECTRAL REGION FOR WHICH CLASS IS ESPECIALLY VALUABLE 
+The I-N film is available as ‘‘Kodak High Speed Infra Red Film."' 
*Class O represents unmodified spectral sensitivity of silver halide. 


This chart shows the spectral sensitivity regions of the various classes of sensitizing and the emul- 
sion types in which they can be furnished. 


In the more advanced aspects of photography, more ideas emulsions that has expanded and greatly improved over 
are floating around today than can be prudently enu- the years. Today there are 111 of these “Kodak Spectro- 
merated in a freely circulating journal. Suddenly the poten- _scopic’’ emulsions, and practically every astronomer knows 
tialities of the photographic emulsion have become highly to order by a combination of a numeral and a letter, as 
exciting. Suddenly many people whose interests hitherto explained in the above charts, e.g. I-C, IV-L, 103a-Z. 
have been concentrated on other fields have to learn the Maybe one of these combinations just neatly fits one of 
most they can about photography in the least time. the brave new ideas now being born. Maybe something 
First off, they should familiarize themselves with the two is required that isn’t in the charts yet. Maybe it can be 
charts above. The charts illustrate how specifically we can made, and maybe it can’t. What needs to be impressed on 
“tune” a photographic emulsion to its purpose. Many years __ the authors of the ideas and their principal helpers is that 
ago, the astronomers (and others who have to define pho- very early in their photographic thinking they ought to 
tographic characteristics with high precision) came and address a few questions to: 
told us their needs. We responded with a large family of 
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Book Reviews__ 


Elements of Pure and Applied Mathe- 
matics, by Harry Lass, McGraw-Hill, 
New York, 1957, 591 + ix pp. $7.50. 
Reviewed by Patrick BENSON 
TEMPO, General Electric Company 


This book is a mathematical compen- 
dium comparable to those of Reddick and 
Miller, Margenau and Murphy, Pipes, 
Wilie, etc.; and, like those, it has certain 
merits and defects. In this case, how- 
ever, it appears that the former far out- 
weigh the latter. The book is well named. 
Dr. Lass sticks to the elements of mathe- 
matics; at the same time, however, he 
provides a number of practical examples 
which will keep engineers and physicists 
interested in the material. Furthermore, 
his choice of examples is timely. For 
example, he includes practical discussions 
of pursuit problems; topics from game 
theory, group theory; applications of the 
Laplace transform, nonlinear mechanics, 
and so on, all of which are becoming 
more and more important in today’s 
technology. 

An advantage of the book is the number 
of examples worked cut, all of which are 
wisely chosen. While, additionally, a 
set of problems is given at the end of each 
chapter, answers to these sets are not 
provided. The present reviewer feels 
that one of the more important functions 
of a good compendium of this type is its 
use as a reference text or a self-study text. 
In either case a list of answers is an im- 
portant asset. 

To take the work in chapter order, the 
first deals with determinants, linear 
equations and matrices. In the second 
chapter, the essentials of vector analysis 
are covered. Here Dr. Lass includes the 
discussion of the kinematics of a rigid 
body with one fixed point, and the dis- 
cussion of pursuit problems. Chapter 3, 
Tensor Analysis, starts with the concept 
of generalizing the notion of a vector and 
leads into a standard discussion of tensors. 
Some of the applications include Rieman- 
nian geometry, Hamilton’s equations of 
motion and the Navier-Stokes equation 
of hydrodynamics. The second and third 
chapters are unusually well presented. 

In Chapter 4, those elements of complex 
variable theory are included which would 
normally be encountered in an elementary 
senior or first year graduate course, but 
there is greater emphasis on rigorous proof 
than is usually found in textbooks of this 
type. Many of the problems are im- 
portant in engineering. In particular the 
application of the Schwartz-Christoffel 
transformation to electrostatic problems 
is included. 

In Chapter 5 topics of ordinary and 
partial differential equations are treated 
and some of the typical solution functions 
(e.g., hypergeometric, Laguerre and Her- 
mite polynomials) are introduced. Un- 
fortunately, a more complete discussion 
of Eigenvalues and Eigenfunctions, which 
could have been related to the solutions of 
boundary value problems, is omitted. 
They are mentioned only briefly at the 
beginning of the book. Chapter 6 in- 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


cludes the elements of orthogonal poly- 
nomials, Fourier series and Fourier In- 
tegrals. The author points out that these, 
combined with the discussions under 
complex variables and the calculus of 
variations, usually provide sufficient back- 
ground for the student to understand 
Wiener’s results concerning linear pre- 
dictors and filters. Some discussion of 
nonlinear differential equations involving 
use of Fourier series is also included. 

In Chapter 7 he deals with the Stieltjes 
Integral, the Laplace transform and the 
calculus of variations. The treatment is 
brief but there seems to be enough material 
to enable the student to understand 
applications to the problems of physics and 
engineering. 

In Chapter 8 two of the more funda- 
mental, and consequently much more 
difficult, topics of mathematics are in- 
cluded, namely group theory and algebraic 
equations. 

Chapter 9 contains discussions of the 
elementary concepts of probability theory, 
the chi-squared distributions of statistics, 
and a brief discussion of game theory and 
the Monte Carlo method. Since it is 
becoming more apparent than ever that the 
engineer must have a thorough knowl- 
edge of probability theory if he is to keep 
abreast of such fields as information theory 
and filter theory, this is an important 
contribution in a textbook. 

The final chapter of the text deals with 
many basic aspects of real variable theory. 
It is over 100 pages long and, taken by 
itself, could comprise a text on real 
variable theory. 

In general the text is written in such a 
way that each chapter is independent. of 
the others, although there are several 
cross references. The reader who is 
interested chiefly in obtaining a cursory 
knowledge of the material will probably 
find that the text contains most of the 
elementary topics found in existing applied 
mathematics texts. On the other hand a 
sufficient degree of rigor is introduced to 
make the treatment more palatable to 
the pure mathematician. 


Chemistry Problems in Jet Propulsion, 
by S. S. Penner, Pergamon Press, 
New York, 1957, 394 + xiv pp. $12.50. 

Reviewed by JoHn B. FENN 
Project squip, Princeton University 


It is » long way from the principles which 
put an electron into an orbit o! the Bohr 
atom to the practice which pushed a 
satellite into an orbit of Earth: Pro- 
fessor Penner has written an excellent 
guidebook for this ambitious journey in 
chemical science and technology. The 
trip has been broken into three stages. 
The first is brisk passage through the 
development of modern structural chem- 
istry. It begins with spectroscopic ob- 
servations. By way of the elements of 
quantum theory it arrives at a description 
of chemical bonds and their implications 
as energy sources and sinks. This is home 
territory for the author. His familiarity 
with the terrain is apparent. The trav- 


eler’s road is smooth and clearly murked 
except for the last few steps. Compared to 
the highway which precedes it, the sk«-tchy 
account in Chapter VII of chemical vom- 
pounds as propellants is rather like a 
detour. 

The second stage of the journey «om- 
prises a tour through classical thermo- 
dynamics and thermochemistry. It -gins 
with a briefing on the mathematics of 
partial differentiation and the prop: rties 
of differentials. It includes instruciions 
for computing equilibrium tempera‘ ures 
and compositions under both static and 
flow conditions. It concludes wiih a 
treatment of the calculation of ther- 
modynamic properties of gases from =! ruc- 
tural and spectroscopic data, where «vain, 
with evident relish, the author is i: his 
element. 

The last lap of the journey leads thr \ugh 
the relatively unknown territory of rate 
processes. After an introduction to «las- 
sical chemical reaction kinetics, consvrva 
tion principles in reacting flow systems : 
transport processes in gases, a numbet 
problems are treated in detail. Tl 
include diffusion flames, premixed lami 
flames and reactions during expand 
flow. Following a discussion of simil: 
and dimensional analysis, the aut 
finally considers some of the real and «« 
plicated technical problems in liqu 
propellant rockets which have not su 
cumbed to exact analysis. There 
short stopovers at combustion instabilit 
scaling, transient phenomena, car! 
formation, combustion of sprays, liqu 
monopropellant burning and combust 
volume requirements. 

Thanks to his terse and lucid styl 
Professor Penner has been able to cral 
a surprising amount of information into 
this relatively small guidebook. The 
numerous and extensive tabulations of use- 
ful data are extremely welcome. The 
cited references are adequate charts fot 
excursions from the main route for more 
detailed views of the scenery along the 
way. Be he an eager student or jaded 
expert, the traveler cannot fail to be 
stimulated into the realization that there 
is still much country to be settled and 
many new trails to blaze. 

On the other hand, the chosen route 
bypassed many attractions. This tourist 
was disappointed at the absence of an 
account of many of the interesting prob- 
lems relating to solid propellants. He re- 
grets the omission of any treatment of pre- 
mixed turbulent flames or homogeneous 
combustion as in the stirred reactor con- 
cept. The discussions of detonation and 
heterogeneous kinetics were tantalizingl) 
brief. The philosophy which justified 
rather extensive review of quantum theory 
could surely have allowed a section on the 
elements of the kinetic theory of gases. 
To keep a book within reasonable size or 
a course within reasonable scope, an autho! 
and teacher is forced to be selective. His 
choice of topics is a matter of taste to 
which there can be little objection. [low- 
ever, the text contains so much which 
is so well done that one cannot help but 
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ENGINEERS 
& SCIENTISTS 


Here is your opportunity to grow 
with a young, expanding subsidiary 
of the Ford Motor Company. Out- 
standing career opportunities are 
open in Aeronutronic’s new RE- 
SEARCH CENTER, overlooking the 
Pacific at Newport Beach, and the 
facility in Glendale, California. 
You will have all the advantages of a 
stimulating mental environment, 
working with advanced equipment 
in a new facility, located where you 
can enjoy California living at its 
finest. 


PhD and MS RESEARCH SPECIAL- 
ISTS with 5 to 7 years’ experience in heat 
transfer, fluid mechanics, thermodynamics, 
combustion and chemical kinetics, and 
thermoelasticity. To work on theoretical 
and experimental programs related to re- 
entry technology and advanced rocket 
propulsion. Specific assignments are open 
in re-entry body design, high temperature 
materials studies, boundary layer heat trans- 
fer with chemical reaction, thermal stress 
analysis, and high temperature thermo- 


dynamics. 

PROPULSION ENGINEERS with 5 years’ 
experience in liquid and solid rocket design 
and test. Familiarity with heat transfer 
problems in engines desirable. To work on 
program of wide scope in R & D of advanced 
concepts in rocket engine components, and 
for missile project work. 

ADVANCED AERODYNAMIC FACIL- 
ITY DESIGNER. Advanced degree de- 
sired. To supervise work in design and in 
instru t of d aerodynamic 
test facilities such as shock tubes, shock 
tunnels, plasma-jets, and hyper-velocity 


STRI 
TRUCTURAL ANALYSIS SECTION 
SUPERVISOR with 8 to 10 years’ experi- 
ence, including supervision, in the missile 
field. Graduate degree for design and 
analysis required. Will be required to 
apply knowledge of high temperature ma- 
terials and methods, thermal stress, 
etc. to advan 
e-entry bodies, and space hicles 
FLIGHT TEST & INSTRUMENTATION 
ENGINEERS with 5 to 10 years’ experi- 
ence in laboratory and flight test instru- 
mentation techniques. Will develop tech- 
niques utilizing advanced instrumentation 
associated with space vehicles. 
THEORETICAL AERODYNAMICIST. 
Advanced degree and at least 5 years’ ex- 
rience in high- aerodynamics 
ledge of viscid and inviscid gas flows 
. To work on program leading to 
advanced missile configurations. Work in- 
volves analysis of the re-entry of hypersonic 
missiles and space craft for determining 
optimum configuration. 
DYNAMICIST. Advanced degree, ap- 
plied mathematics background, and ex- 
perience in missile stability analysis de- 
sirable. Work involves re-entry dynamics 
of advanced vehicles and dynamic analysis 
of space craft. 
ENGINEER or PHYSICIST. With ex- 
perience in the use of scientific instru- 
ments for making ysical measurement. 
Work related to flight test and facility 
instru degree de- 
sired with minimum of 3 years of related 
experience. 


Qualified applicants are invited to send res- 
umes and inquiries to Mr. L. R. Stapel. 


AERONUTRONIC SYSTEMS, 
INC. 


a subsidiary of Ford Motor Company 


1234 Air Way, Bidg. 19, Glendale, Calif. 
CHapman 5-6651 


wish there ‘ai been more. yey 

Of more concern than omissions are the 
implications of this book in terms of a 
growing problem in engineering education. 
It is food for thought that incumbent 
upon a graduate course in aeronautical 
engineering is the need to cover elementary 
quantum mechanics, thermodynamics and 
even some basic mathematics. At least 
half of the book, were it separated from 
the title page and preface, could never 
be identified by a casual observer as 
jet propulsion course material. Famili- 
arity with quantum theory and _ther- 
modynamics is fully as important to other 
kinds of engineers and scientists. One 
‘an reasonably ask whether a student who 
does not have a working acquaintance 
with these basic disciplines is in any posi- 
tion to have already decided to concentrate 
on jet propulsion or any other specialized 
technology. By the same token it would 
seem a pity that Professor Penner’s 
admirable presentation of this funda- 
mental material should be available 
only to engineering students majoring in 
jet propulsion. 


Vibration and Impact, by Ralph Burton, 
Addison-Wesley Publishing Co., Inc., 
Mass., 1958, 310 + x pp. $8.50. 

Reviewed by D. EK. Hupson 
California Institute of Technology 

The two classic engineering vibration 
books of Timoshenko (1928) and Den 
Hartog (1934) established a pattern which 
has well served the student and the 
teacher. Through successive revisions 
these books have increased their usefulness 
as standard texts and references in the 
field. Recent developments, however, in 
both the subject matter itself and in its 
place in engineering education have pro- 
vided incentives for a number of new 
books with significantly different ap- 
proaches. The present book is an example 
of one type of response to these new re- 
quirements. 

The major principles behind this book 
appear to be: (1) To provide an elemen- 
tary treatment of the subject suitable for 
undergraduate students with a minimum 
background of dynamics and mathematics; 
(2) to emphasize the physical aspects 
rather than the mathematical, and the ap- 
plications rather than the basic theory; 
(3) to include brief introductory treat- 
ments of a number of topics of current in- 
terest, such as shock problems, random 
loads, control systems, instrumentation 
and acoustics. 

The book is designed as a textbook; 
many of the explanations would require a 
considerable expansion by the teacher, and 
a number of the problems would need 
much class discussion. These are of 
course advantageous features from the 
teaching standpoint, provided that close 
class supervision is available. 

The laudable desire to introduce many 
topics of current interest sometimes leads 
the author to oversimplifications. The 
treatment of random loads, for example, 
seems inadequate even as an introduction, 
and the amount of attention given to the 
impact problem hardly justifies the prom- 
inence of the word in the title. In some 
cases, as in the group of formulas for drag 
forces, the restrictions on the use of the 


information have not been fully stated, 


and some of the examples of damping cal. 
culations seem to imply an ease and exaet 
ness of calculation that might not 
realized in practice. The definition 
equivalent viscous friction used by 
author at one point seems to be a differey 
one than customarily used, and for hyst 
esis damping it is not clear that the resv 
obtained by the two approaches are simil 

A commendable feature of the bouk 
the introduction and frequent use o! 
concept of dynamic stability, and this 
ables the author to show the unity bet wee 
vibration problems and control syst 
problems. Other favorable feature 
clude the interesting brief historical int 
duction and the annotated reference |is 

Mechanical Vibration courses are n 
quite common in our undergraduate « 
gineering curricula, and this book may 
find a useful place there for the bric! in- 
troductory course. The student who 
wishes ultimately to fully develoy) his 
knowledge of the field, however, would 
probably be better advised in his under- 
graduate program to concentrate on |asic 
dynamics and mathematics, whicl: he 
could then use as a foundation for ad- 
vanced studies. 


Book Notices 


Nuclear Engineering, by C. F. Bonilla 
McGraw-Hill, New York, 1957, 850 py 
$12.50. This excellent text was pre 
pared to give the engineer the details o 
nuclear powerplant design. Because of 
the vastness of the subject matter, differ- 
ent outstanding specialists wrote about 
their respective fields. The authors and 
their subjects are: 

Theodore Baumeister—Power Gener- 

ation 
Charles F. Bonilla—The Flow of Fluids; 
Heat Removal 

John R. Dunning—Introduction 

Gioacchino Failla—Basic Concepts of 
Radiation Protection 

Alfred M. Freudenthal—Thermal-Stre 
Analysis and Mechanical Design 

William W. Havens Jr.—Nuclear Pa 

ticles; Nuclear Physics 

John W. Hoopes Jr.—Instrumentatic 

and Control 

R. Wayne Houston—Introduction; Ele 

mentary Reactor Physics 

George L. Kehl—Metallurgy of Ura 

nium and Uranium Alloys 

John Landis — Nuclear-Reactor Types 

John G. Palfrey—Legal Aspects of Nu- 

clear Power 

David C. Peaslee—Shielding of Power 

Reactors 

Chien-Shiung Wu—Particle Detection 

The First One Hundred and Fifty Years, 
by The House of Wiley, John Wiley & 
Sons, 1957, 242 pp. $7.50. An interest- 
ing book relating the history of the Joht 
Wiley & Sons, Inc., publishing house be- 
tween the years 1807-1957 and written by) 
many leaders from industrial firmsand edt 
cational institutions. 

The Cosmic Radiation, by J. E. Hooper 
and M.Scharff, Wiley, New York, 1958, 172 
pp. This little volume should be useful 
to those who wish to get a bird’s eye 
view of the cosmic radiation techniques 
for measurements and physical effects. It 
is not a popular book, but requires 4 
good background in modern physics 
including some nuclear physics. 
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e labs of Annapolis to the launching 


the early spac 
pads at Cape Canaveral — setting the standards for 


the control of high pressure fluids...the regulators 
which served as the critical control on World War II 
flame throwers and torpedoes...and made missile 
development work possible... Grove High Pressure 


Regulators... frequently imitated but never equalled 


PROVE VALVE and REGULATOR COMPANY 


829 Hollis St., Oakland 8, California - 2559 W. Olympic Blvd., Los Angeles 6, California 
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SELECT 
ASSIGNMENTS IN 
SPACE PROPULSION 


Chance Vought offers opportunity to join 
small-group research and development 
of propulsion systems for space and for 
utilization of nuclear power. Required 
for each opening: advanced Engineering 
or Physics Degree, at least 6 years 
applicable experience. 


Senior Propulsion Specialist. To conduct 
research and development in the fields 
of advanced heat transfer and space 
environmental systems; pursue technical 
studies of specific cooling methods for 
hypersonic vehicles; interpret advanced 
hypersonic cooling data and advise 
Chief of Propulsion. 


Senior Propulsion Engineer. To report 
to Chief of Propulsion in staff capacity; 
conduct special technical studies to 
optimize energy sources for space 
propulsion systems (both primary and 
auxiliary power); supervise applied 
research and development in the space 
engine and auxiliary power and control 
fields; monitor advanced studies of 
propulsion requirements and methods. 


S. J. Townsend, 
Chief of Propulsion, 
Dept. JP-6 


AIRORAFT 


MARQUARDT 


Professional Personnel 
Requisition 


Test Facilities 
ANALYSIS ENGINEER 


To analyze problems associated with 
engine test equipment at Marquardt 
Jet Laboratory, the West's largest jet 
engine testing facility. 


Will perform highly analytical calcu- 
lations involving thermodynamic, aero- 
dynamic and stress analysis problems 
as related to such equipment as high 
pressure storage tanks; compressors; 
cooling towers; heat exchangers; water 
systems; temperature, pressure and 
flow instrumentation and controls. 


Excellent opportunity. Requires gradu- 
ate engineer with at least 5 years’ 
experience. 


For information, contact: 
Jim Dale, Manager 
Professional Personnel 
_ Marquardt Aircraft Co. 
He 16556 Saticoy Street 


Technical Literature Digest 
M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Heat Transfer and Fluid 
Flow 


On Transition Experiments at Moder- 
ate Supersonic Speeds, by M. V. Mor- 
kovin, J. Aeron. Sci., vol. 24, July 1957, 
pp. 480-486. 

Laminar Jet Mixing of Two Compres- 
sible Fluids with Heat Release, by S. I. 
Pai, J. Aeron. Sci., vol. 23, Nov. 1956, p. 
1012. 

Normal Conduction Effects on Heat- 
transfer Data During Transient Heating 
of Thin-skin Models, by Morton Cooper 
and Edward E. Mayo, J. Aeron. Sci., 
vol. 24, June 1957, pp. 461-462. 

Characteristics and Sound Speed in 
Nonisentropic Gas Flows with Nonequilib- 
rium Thermodynamic States, by E. L. 
Resler Jr., J. Aeron. Sciences, vol. “24, 
no. 11, Nov. 1957, pp. 785-790. 

Experiments on Interaction between a 
Traveling Shock Wave and a Turbulent 
Jet, by Darshan S. Dosanjh, J. Aeron. 
Sciences, vol. 24, no. 11, Nov. 1957, pp. 
838-844. 

A Method of Heating Matter of Low 

Ep1tor’s Note: Contributions from E. R. 
G. Eckert, J. P. Hartnett, T. F. Irvine Jr. 
and P. J. Schneider of the Heat Transfer 
Laboratory, University of Minnesota, are 
gratefully acknowledged. 


Density to Temperatures in the Range 
10° to 10°K, by F. B. Knox, Australian J, 
Phys., vol. 10, 1957, p. 221. 

Heat Transfer Media for Use at Ele- 
vated Temperatures, by P. A. Rotten- 
burg, T’rans., Inst. Chem. Engrs., vo'. 35, 
1957, p. 21. 

Theoretical Analysis of Viscous Flow in 
Narrowly Spaced Radial Diffuser, by 
H. W. Wollard, Trans. ASME (J. - ‘ppl. 
Mech.), vol. 24, 1957, p. 9. 

Power from Solar Energy, by J. I. 
Yellott, Trans. ASME, vol. 79, ‘957, 
p. 1349. 

Prediction of Heat Transfer Burnout, 
by Louis Bernath, Chem. Engng. Pvogr., 
vol. 53, 1957, p. 226. 

Recorrelation of Data for Convective 
Heat Transfer Between Bases and Single 
Cylinders with Large Temperature Dif- 
ferences, by W. J. M. Douglas and +. W. 
Churchill, Chem. Engng. Progr., vol. 53, 
1957, p. 226. 

Pressure Drop during Forced-circula- 
tion Boiling, by Max Jakob, George Lep- 
pert and J. B. Reynolds, Chem. Engng. 
Progr., vol. 53, 1957, p. 226. 

Measurement and Prediction of Den- 
sity Transients in a Volume-heated Boil- 
ing System, by R. P. Lipkis, C. Liu and 
N. Zumer, Chem. Engng. Progr., vol. 53, 
1957. 


PRECISION 
ACHINE PARTS 


Small parts requiring very close tolerances and numerous machine 
operations have been our specialty during the past half-century! 
We invite your quotation requests and suggest that you write for 
our illustrated brochure to see how our facilities can be put to 


MACHINE WORKS 


4635 WEST LAKE ST., CHICAGO, ILLINOIS _ 
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Better Things for Better Living... through Chemistry 


Lino-Writ 4...fastest to use, fastest to process 


__Lino-Writ 4 is the fastest, toughest, thinnest 
paper you can use. Now with Du Pont’s new 
low-odor Rapid Processing Chemicals Kit, it 
offers the additional advantages of fastest proc- 
essing without sacrificing fine detail. 

With this combination you can now use proc- 
essing temperatures up to 120° F. and get better 
trace densities, cleaner records, higher photo- 
graphic speeds and faster processing throughout. 


REG. U.S, PAT. OFF. 


No other photorecording paper-chemical team 
can offer these across-the-board “plus” features. 
Ask your Du Pont Technical Representative 
about Lino-Writ 4, the new low-odor Rapid 
Chemicals Kit and the entire Lino-Writ line. 


E. I. du Pont de Nemours & Co. (Inc.), Photo 
Products Department, Wilmington 98, Delaware. 
In Canada: Du Pont Company of Canada (1956) 


Limited, Toronto. 
ively by Phototypography 
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«BUILDERS OF POWER FOR OUTER SPACE 


CONTROL ENGINEERS 


(Electric » Servo « Valves) 


Move into Large Rocket Engineering 
and put yourself way ahead in your field — 


Help us to automate millions of horsepower designed into a 
4 jet-size package—the High-Thrust Rocket Engine. Here are the 
fields : 

_ The Electrical System includes Ground Support and Check- 

out Equipment which must be operable by military personnel. 

Aboard the missile, engine controls must be carefully isolated 

_ from other missile systems. Miniaturization is striven for, but 

- never at the expense of reliability in extremes of temperature, 

_ vibration and acceleration. You'll cover all aspects of circuitry, 

deal with every branch of weapons systems. 

_ Servo-mechanisms offer a broad spectrum—electronic, pneu- 
matic, mechanical, hydraulic. Your analytical ability will be at 
“3 a premium here, to evaluate methods of Mixture Control, Thrust 


efficiency, heat, cold, G, vibration etc. tad ll be free of routine 


a 4 details, able to apply your training and experience toward a high 


a level of professional growth. 
_ Valves run to 6” diam. and up, with very high pressures and 
_ flow rates, extremely rapid action, temperatures down to —300°F. 
This is where the real advanced work in controls is being done. 
_ Join the trailbreakers. Write, giving your background: Mr. 
A. K. Jamieson, Rocketdyne Engineering Personnel, 6633 
Avenue, Canoga California. 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 


eeeese 


Flow of Liquid Helium I through Nar- 
row Channels, by D. C. Champeney, 
Proc. Physics Soc. vol. 69, 1956, p. 1112. 

Flow of Compressible Fluid through 
Sudden Enlargement in Pipe, by W. B. 
Hall and E. M. Orme, Proc. Inst. Mech, 
Engrs., vol. 169, 1955, p. 1007. 

Velocities and Temperatures on Axis 
of Downward Heated Jet from 4-inch 
Long-radius ASME Nozzle, by R. Kiaak, 
Heating & Air Conditioning, vol. 29, 1957, 
p. 137. 

Pressure Drop and Flow Characteris- 
tics of Short Capillary Tubes at Low Reyn- 
olds Numbers, by F. Kreith and R, 
Eisenstadt, Trans. ASME, vol. 79, i957, 
p. 1070. 

An Empirical Correlation for Velocity 
Distribution of Turbulent Fluid Flow, by 
B. F. Ruth and H. H. Lang., J., Al‘ hE, 
vol. 3, 1957, pp. 117-120. 

Determination of Burnout — of 
Polyphenyl Coolants, by K. Sato, ojel- 
General Corp., (US Atomic Plage m.), 
AGC-AE-32, Feb. 1957. 

Measurement of Local Heat Transfer 
Coefficients with Sodium-potassium Eu- 
tectic in Turbulent Flow, by K. D. Kuvzen 
and T. R. Bump, Nuclear Sci. and Enyng., 
vol. 2, 1957, p. 181. 

The Heat-transfer and High-tempera- 
ture Properties of Liquid Alkali Metals, by 
I. I. Novikov, J. Nuclear Energy, vol. 4, 
1957, p. 387. 

Mercury Vapour Transfer Studies, 
R. W. Maxwell and J. Anderson Storro 
Chem. Engng. Sci, vol. 6, 1957, p. 204. 

Thermal Conductivity of Some Re- 
fractory Materials, by J. F. Clements a 
J. Vyse, Trans., Brit. Ceramics Soc., v 
56, 1957, p. 296. 

Vapor Pressure Data More Common 
Elements, by R. E. Honig, RCA Rev., 
vol. 18, 1957, p. 195. 

Determination of Thermal Conduc- 
tivities of Metals by Measuring Transient 
Temperatures in Semi-infinite Solids, by 
8S. T. Hsu, Trans. ASME, vol. 79, 1957, 
p. 1197. 

Thermal Relaxation in Carbon Dioxide 
as a Function of Temperature, by F 
Douglas Shields, J. Acoustical S 
America, vol. 29, 1957, p. 450. 

Radiative Transfer in Discontinuous 
Media, by R. G. Giovanelli, Australian 
Phys., vol. 10, no. 2, June 1957. 

Some Transition Forms between Rapid 
and Tranquil Flow, by N. V. Feodoroff 
and B. A. Bakhmeteff, T'rans., N. 
Acad. Sci., vol. 19, 1957, p. 518. 

Rotational Waves in Turbulent Liquid, 
by R. V. L. Hartley, J. Acoustical Soe. 
America, vol. 29, 1957, p. 195. 

Material Transport in Turbulent Gas 
Streams: Radial Diffusion in Circular 
Conduit, by S. Lynn, W. H. Corcoran and 
B. H. Sage, J. AIChE., vol. 3, 1957, p. 11. 

Numerical Solution of Equations of the 
Diffusion Type with Diffusivity Concen- 
tration-dependent II, by J. R. Philip, 
Australian J. Phys., vol. 10, 1957, p. 29. 

Mass Transfer between Co-current 
Fluid Streams and Boundary Layer Solu- 
tions, by Owen E. Potter, Chem. Engrg. 
Sci., vol. 6, 1957, p. 170. 

Improving Gas Turbine Performance, 
3: Cooling Turbine Blades, by F. J. 
me and N. Bell, Engng., vol. 304, 1:57, 

». 183. 


The Reflection of Shock Waves from an 
Orifice at the End of Duct, by George 
Rudinger, Project Squid. Tech. lep. 
CAL-73-P (ASTIA AD 142826), Sept. 
1957, 28 pp. (available only on mi:ro- 
card). 
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ACOUSTIC 
NOISE 


When huge rocket engines lift a 

ballistic missile or satellite vehicle 

into the air, the thundering noise of 

the engines violently shakes intricate 
electronic components and the missile 
structure itself. Acoustic noise is 

one of the many environments that must 
be thoroughly explored and understood 
by the missile designer. 

But it is expensive and impractical 

to conduct acoustic noise tests by 
operating rocket engines on the ground. 
Noise—of the same frequencies and 
amplitudes as rotket engine noise—must 
be produced in the laboratory. 


Engineers at Avco’s Research and 
Advanced Development Division have met 
this problem by designing a versatile 
high-intensity acoustic noise generato: 
Based on electromagnetic principles, 

it creates noise environments of th¢ 

same frequency spectra as rocket 

or jet engines. 

This successful effort is now part o1 

the division’s complete environmenta 
testing laboratory, where it is producing 
information for Avco scientists and 
engineers responsible for ICBM nose cones. 


Acoustic noise work is one of the many 
development programs that proceed 
simultaneously with theoretical and 
applied research at RAD. The mutua! 
stimulation of scientists and engineers 
in many disciplines occurs naturally 

at Aveo, where the search for new 
knowledge and its application 

go forward under one roof. 


Research and Advanced Development is 
more than a descriptive title at Avco. 
It is a concept that promotes creativity. 


For information on unusual career ».. 
opportunities for exceptionally 

qualified scientists and engineers, 

write to: Dr. R. W. Johnston, 

Scientific and Technical Relations, 

Avco Research and Advanced Development Di 
201 Lowell Street, Wilmington, Mass. 


J 
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If you are interested in new 
approaches to fins, airfoil sections, 
access doors, or bulkheads, 


information awaits your request. 


Such high strength aluminum parts | 


for missiles are accurately cast with 
the unusual foundry techniques of 


Morris Bean & Company, 
Yellow Springs 6, Ohio. a 


A Method for Calculating Turbulent 
Boundary Layer Development in Super- 
sonic and Hypersonic Nozzles Including 
the Effects of Heat Transfer, by Jerome 
Persh and Roland Lee, NAVORD Rep. 
4200 (Aeroballistic Res. Rep. 320), June 
1956, 10 pp., 31 fig. 

A Note on Turbulent Boundary Layer 


| Allowances in Supersonic Nozzle Design, 


by E. W. E. Rogers and B. M. Davis, 
Gt. Brit., Aeron. Res. Council, Current 
Paper 333, (A.R.C. Tech. Rep. 18490), 
1957, 14 pp. 11 fig. 

Converging Cylindrical Shocks in Mag- 
netohydrodynamics, by G. B. Whitham, 
N. Y. Univ., Inst. Math. Sci. Rep. IMM- 
NYU 237, Jan. 1957, 13 pp., 2 fig. 

The Penetration of a Shock Wave into a 
Magnetic Field, by J. M. Burgers, Univ. 
Maryland, Inst. Fluid Dynamics Appl. 
Math., Tech. Note BN-82 (AFOSR-TN- 
56-486; ASTIA AD 136511), June 1957, 
84 pp. 

On Exact Solutions of One Dimensional 
Flow Equations of Magnetogasdynamics, 
by S. S. Pai, Univ. Maryland, Inst. Fluid 
Dynamics Appl. Math., TN BN-82 
(AFOSR-T N-56-486; ASTIA AD 97370), 
Sept. 1956, 12 pp. 

On the Dissociation Rate of a Diatomic 
Gas, by Gianni Jarre, Torino, Poliiechnico, 
Lab. di Meccanica Applicata, TN 

(AFOSR-TN-57-528; ASTIA AD 1365- 

3), July 1957, 12 pp. 

Research on the Design of Hypersonic 
Nozzles and Diffusers at High Stagnation 
Temperatures, Part I: Design and Con- 
struction Problems of a Hypersonic Facil- 
ity and Preliminary Investigation of Liquid 
Film Cooling, by R. Hermann, H. Leitinger 
and W. L. Melnik, Wright Air Dev. 
Center, TN 55-507, Univ. Minn., Rose- 
mount Aeron. Labs., Res. Rep. 127), 
March 1955, 182 pp. 

Vibrational Relaxation Phenomena in a 
Shock Tunnel, by W. J. Lick and T. Y. Li, 
Rensselaer Polytech. Inst., Dept. Aeron. 


| Engng., TR AE 5704 (AFOSR-TN-57- 


346; ASTIA AD 132419), June 1957, 
21 pp. 7 fig. 

Gas Properties Behind Shocks at Hy- 
personic Velocities, II: Introduction to 
General Thermodynamics of a Real Gas, 
by Adolf R. Hochstim, Convair, San Diego, 
Physics Group, Rep. ZPh- 003, May 1957, 
50 pp. 

Approach to Equilibrium Ionization Be- 
hind Strong Shock Waves in Argon, by 
Harry Petschek and Stanley Byron, 
Cornell Univ. Grad. School Aeron. Engng., 
1955, 49 pp., 16 fig. 

Thermodynamic Properties and One 
Dimensional Flow of a Partially Ionized 
Monatomic Gas, by J. Pomerantz, E. 
Winkler and A. E. Seigel, NAVORD Rep. 
4222, (Aeroballistic Res. Rep. 328), Feb. 
1956, 19 pp., 7 fig. 

Thrust Reversers for Jet Aircraft, by 
F. J. Stimler and J. F. McDermott, SAE 
Prepr. 112, Apr. 1957, 7 pp., 14 fig. 

Experimental Velocity and Temperature 
Profiles for Air in Turbulent Pipe Flow, by 
C. A. Sleicher Jr., ASME Paper 57-HT-9 
(ASME-AIChE "Heal Transfer Confer- 
ence, University Park, Pa., Aug. 11-15, 
1957), Aug. 1957, 10 pp. 

Heat Transfer in Liquid Metals, by 
P. S. Lykoudis and T. S. Touloukian, 
ASME Paper 57-HT-16 (ASME- AIChE 
Heat Transfer Conference, University Park, 
Pa., Aug. 11-15, 1957), 11 pp. 

Heat Transfer from Superheated Vapors 
to a Horizontal Tube, by G. Balekjji jian and 

. L. Katz, ASME Paper 57-HT-27 
(ASME-AIChE Heat Transfer Conference, 
University Park, Pa., Aug. 11-15, 1957), 
Aug. 1957, 18 pp. 


No ordinary 


SYSTEMS 
ANALYST 


will do for this job... 


The engineer we need has superic 
creative ability and an analytic: 
mind. He now has senior status— 
Mechanical, Chemical, Electric: 
Engineer or Engineering Physici: 
(preferably with advanced degree 
who is versed in classical vibratio 
analysis, as well as feedback analysi 
for control of systems composed o 
Heat Transfer, Fluid Mechanics an: 


“Thermodynamic processes. 


The Dynamical system of th 
High-Thrust Rocket Engine is one 0 
extraordinary interest, exceptiona 
performance. No matter what you 
achievements have been, you'll fin 
new interests at Rocketdyne. You wil 
be confronted with the analysis o 
design and operational problems ot! 
the rocket engine as a dynamic sys 
tem. You must develop valid mathe 
matical models of both systems and 
components, using advanced physica 
concepts, and empirical data. These 
must be combined using digital ccm 
putation and analog simulation. 

You'll work with the leading pro 
ducer in the nation’s fastest growing 
industry. Rocketdyne builds the high- 
thrust rocket propulsion systems fo: 
America’s major missiles. 

We know we can show you, in a 
personal discussion, all the opportu- 
nity you can wish for. Write to: Mr. 
A.K. Jamieson, Rocketdyne Engineer- 
ing Personnel, 6633 Canoga Avenue, 
Canoga Park, California. 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION, IN 


BUILDERS OF POWER FOR OUTER SPACE 
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Moving upstairs’ through the application of the 
same sound principles which have made 
S@érvomechanisms, Inc. a leading supplier of avionic 
equipment for modern-day manned aircraft 
3M1 is today successfully meeting the challenges 
nthe equally complex field of missile 
systems. Major design and production projects 
by 5 MI hove resulted in a family of 
a@vanced products for both airborne and 
gfound support missile requirements. 


Te same high standards of reliability and 
mpintainability inherent in SM] 
Mechatronic Design Phi 

mOcuies céveioped 

nave 

for missile 

We woulc welcome the 

opportunity fo work with you 

on your missile problems. 


. 
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LITIES 


Engine Controls 


Instrumentation 


Ground Checkout Equipment 
Infrared Devices 

Guidance and Autopilot Subsystems 
Production Test Equipment ae 


Complete Qualification 
and Environmental Facilities 


Rotating and Magnetic Components 


Vacuum Film Techniques 


COMPONENTS 


GENERAL OFFICES: A, ON BOULEVARD. HAWTHORNE, CALIFORNIA 


of SMI are availoble in Conoda and throughout 


Advanced Weapons 


broad creative mission at Chance Vought 


MILITARY ASTRONAUTICS: Vought studies toward space research 
vehicles and manned spacecraft include multistaging, space 
communications, nuclear and ionic propulsion, celestial navi- 
gation, and membership on Boeing’s Dyna Soar space glider 
development team. 


ASW RESPONSIBILITY: Under the Office of Naval Research, 
Vought antisubmarine specialists are conducting extensive 
research and development work in ASW detection and clas- 
sification phases. Goal is to bring detection abreast of 
destruction capabilities. 


Advanced Weapons Staff Engineer. Ph.D. preferred, with at 
least 10 years background in guidance or navigation and control 
systems. To develop completely new concepts in guidance, naviga- 
tion, or control systems. 


Electro-mechanical Systems Engineer or Specialist. A.E., E.E., 
or M.E., (advanced degree preferred) with at least 7 years expe- 
rience in autopilot, flight control, stability systems and inertial guid- 
ance systems and design work. To make high-level technical studies 
of various control and stabilization systems for advanced weapons. 


ASW Detection Specialist. Physicist or Electronics Engineer with 
Sonar or electromagnetic detection experience. Familiarity with 
submarine tactics, equipment highly desirable. To devise new 
methods for submarine detection, conduct necessary preliminary 
analyses, and prepare information leading to hardware design for 
laboratory testing. 


Lead Instrumentation Engineer. B.S. or M.S. in E.E., plus 5 
years electronics experience in telemeter applications of flight test 
instrumentation data. To guide a Research and Development group 
in the design and application of electronic gathering and reduction 
systems. 


Hydrodynamicist. B.S. or M.S. in Engineering plus 5-7 years 
experience in hydrodynamics and cavitation of torpedoes and other 
fully submerged vehicles. 

Qualified engineers and scientists who would like to join Vought’s 
development of advanced weapons are invited to inquire. 


A. L. JanRETT, Manager, Advanced 
Weapons Engineering, Dept. JP-5 


INCORPORATED DALLAS, TEXAS 


CHANCE 


Wine 


The Effect of a Variable Evaporation 
Rate on the Ballistics of Droplets, by C.¢. 
Miesse, J. Franklin Inst., vol. 264, no. 5 
Nov. 1957, p. 391. 

Mixing in Laminar-flow Systems, }, 
W. D. Mohr, R. L. Saxton and C. i. 
Jepson, Ind. and Engng. Chem., vol. 49, 
vol. 49, no. 11, Nov. 1957, pp. 1855 -1856. 

Film Boiling of Flowing Subcooled 
Liquids, by E. I. Motte, Ind. and !ngng. 
Chem., vol. 49, no. 11, Nov. 1957, pp. 
1921-1928. 

Ionization in the Trail of High-v-locity 
Pellets, by William S. Partridge «id [, 
Dale Harris, J. Appl. Phys., vol. -'8, no. 
11, Nov. 1957, pp. 1269-1271. 

Solution of the Bloch Equations for 
Determination of Relaxation Times in 
Liquids, by P. S. Hubbard Jr. ano T. J. 
Rowland, J. Appl. Phys., vol. 28, io. 11, 
Nov. 1957, pp. 1275-1281. 

Time Lag between High-speed Pellets 
and the Ionization in Their Trai's, by 
R. A. Davidson and W. S. Par: ridge, 
J. Appl. Phys., vol. 28, no. 11, Nov 1957, 
pp. 1304-1308. 

Ionization by Ultra-speed Pelleis, by 
Charles D. Hendricks Jr., J. Appl. hys., 
vol. 28, no. 11, Nov. 1957, pp. 1339 1341. 

Static Pressure Distribution in the Free 
Turbulent Jet, by David R. Miller and 
Edward W. Comings, J. Fluid \ech., 
vol. 3, Part 1, Oct. 1957, pp. 1-16. 

On the Instability of Small Gas Bubbles 
Moving Uniformly in Various Liquids, by 
R. A. Hartunian and W. R. Sears, J. 
Fluid Mech., vol. 3, Part 1, Oct. 1987, 
pp. 27-47. 

The Laminar and Turbulent Mixing of 
Jets of Compressible Fluid, Part II: 
The Mixing of Two Semi-infinite Streams, 
by L. J. Crane, J. Fluid Mech., vol. 3, 
Part I, Oct. 1957, pp. 81-92. 


Combustion, Fuels and 


Propellants 


Effects of Additives on Flame Propaga- 
tion in Acetylene IV, Univ. Texas, TN, 34 
(AFOSR-T N-56-598; ASTIA AD 1150- 
25), Nov. 1956, 6 pp. 

Experimental Investigation of Tail Pipe 
Burner Design Variables, by W. A. 
Fleming, E. William Conrad and A. W. 
Young, NACA RM. E50K22, March 1951, 
75 pp. (Declassified from Confidential by 
authority of NACA Res. Abs. 122, p. 7, 
12/3/57.) 

Preflight and Flight-test Investigation 
of a 50-percent-magnesium 50-percent 
JP-4 Slurry Fuel in a Twin-engine Ram- 
jet Vehicle, by Otto F. Trout Jr. and 
Thomas L. Kennedy, NACA RM L56C06, 
May 1956, 27 pp., diagrs., photos. (De- 
classified from Confidential by authority 
of NACA Res. Abs. 122, p. 14, 12/3/57.) 


Flight Investigation of a Ram Jet 
Burning Magnesium Slurry Fuel and Hav- 
ing a Conical Shock Inlet Designed for a 
Mach Number of 4.1, by Walter A. 
Bartlett Jr. and Charles F. Merlet, 
NACA RM L56124a, Jan. 1957, 23 pp., 
diagrs., photos. (Declassified from Con- 
fidential by authority of NACA Res 
Abs. 122, p. 14, 12/3/57.) 

Flight and Preflight Tests of a Ram Jet 
Burning Magnesium Slurry Fuel and 
Utilizing a Solid-propellant Gas Generator 
for Fuel Expulsion, by Walter A. Bartlett 
Jr. and William K. Hagginbothom Jr., 
NACA RM 155A24, Apr. 1955, 35 pp. 
diagrs., photos. (Declassified from on- 
fidential by authority of NACA Res. Abs. 
122, p. 13, 12/3/57.) 

Sources of Error in the Determination of 
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Exclusive Produce 


Only Becco can supply you 98% Hydrogen Peroxide 
incommercial quantities—in tank cars or drums. For 
'V.-34 |use both as a monopropellant and an oxidant, super- 
1150- Foure 98% HO, is a stable liquid over a wide range 
of temperatures (B.P.: 299.2°F.; F.P.: 27.5°F.). Its 
high density (Sp. Grav.: 1.44) reduces space require- 
ments per unit weight. Its ease of catalyticand thermal 
al by decomposition makes for simple and reliable systems. 


Yet, these are only a few of the advantages of 98% ae a 
ation #4,0. in propulsion. Becco, as the pioneer in the field, = = =— 
Ram- §'8 especially well equipped to provide you with com- oe 
_2n4 Dolete information on all phases of the use and handling 
(De- fof Concentrated Hydrogen Peroxide. Write today for 
' J ree samples and our newest Technical Bulletin (No. 

Jet $93), “Anhydrous Hydrogen Peroxide as a Propellant”. 
Hav- f Address: Dept. JP-25. 


r A, R 

erlet, 

in Broygere B 

‘Res. 
BECCO CHEMICAL DIVISION Food and Chemical Corporation 

Pe Station lo 7, New York 


ial Sulphur © OHIO-APEX Plosticizers and Chemicals © FAIRFIELD Pesticide Compounds and Organic Ch 


con 


tne FMC CHEMICALS INCLUDE: BECCO Peroxygen Chemicals © WESTVACO Phosphates, Barium and Magnesium 
Chemicols © WESTVACO Alkalis, Chlorinated Chemicals and Carbon Bisulfide * NIAGARA Insecticides, Fungicides and 
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the Characteristics Velocity of Solid Pro. 
liants, by Ellis M. Landsbaum, Calif. 


4 H nst. Tech., Jet Propulsion Lab., M 
Provide instant 
A General Method of Calculating the 


relief with 
eee Performance of Propellant Systems Con- 


EMERGENCY EYE and Ore Other Halogen, by C. A. Se fede 
FACE WASH FOUNTAINS 


Calif. Inst. Tech., Jet Propulsion Lab, 
Mem. 20-144 July 1957, 29 pp. 

A Single Generalization Chart of De- 

tures, by Robert G. Dunn and Bern:::d 

Eye and face hazards are inevitable in Wolfson, Wright Air Devel. Cente, TN 

industry ...so take positive steps to reduce the 57-263 (ASTIA AD 130906), Aug. 1957, 
extent of injury and minimize insurance claims. 
HAWS Emergency Eye and Face Wash Fountains 

flood contaminated areas with water — instantly 

ridding face and eyes of caustics, chemicals and dan- 

gerous particles. This instant relief may well pre- 


My eyes! My face! 


11 pp. 

Study of Some Burner Cross-s«ction 
Changes That Increase Space-heating 
Rates, by Donald R. Boldman and erry 
L. Blackshear Jr.,. NACA TN 4162, Noy. 
1957, 38 pp. 


Model 8933: 
Face spray ring acts simultaneously 
with eye-wash, sending controlled 
streams of water from specially de- 
signed twin fountain heads. 


MARQUARDT 


Professional Personnel 
Requisition 


vent permanent damage — bridging the gap until 
medical aid arrives. Install extra safety —HAWS 
Emergency Fountains and Drench Showers. 


Write today for detailed infor- 
mation on the complete HAWS 
line of emergency facilities. HAWS 
DRINKING FAUCET COMPANY, 1443 
Fourth Street, Berkeley 10, Calif. 
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AERODYNAMICIST 


Internal Flow Systems 


To formulate and test inlet systems and 
exit nozzles for hypersonic air breath- 
ing propulsion devices. Create new 
designs for advanced and unique en- 
gine configurations. Supervise aero- 
dynamic design of hypersonic free jet 
and wind tunnel research facilities. 


Five years’ experience in internal flow 
and fluid dynamics analysis. Re- 
quires Mechanical or Aeronautical 
Engineering degree with appropriate 
advanced study. 


Small work group offers unusual oppor- 
tunity for professional recognition 
and advancement. 
Contact: Jim Dale, Manager 

Professional Personnel Department 


Marquardt Aircraft Co. 
16556 Saticoy Street 


Van Nuys, California ‘ 


COMPLETE 


TESTING 


FACILITIES 


* Qualification Tests 
* Evaluation Tests 

* Performance Tests 
* Environment Tests 


Manufacturers of Metal Boss 
Seals to Military Specifications 


TESTING CO. 
1812 Fleet St., Baltimore 31, Md. 
ORleans 5-8337 ORleans 5-2222 


Final Report on Studies of Turbulent 
Premixed Flames and Methods of Sup- 
pressing Oscillations, by Abbott A. Put- 
nam, Gail M. Clough and Milton J. 
Kenworthy, Wright Air Devel. C:ner., 
Tech. Rep. 56-583 (ASTIA AD 13072), 
Dec. 1956, 125 pp. 

Theoretical Rocket Performance of JP 4 
Fuel with Several Fluorine-oxygen Mix- 
tures Assuming Frozen Composition, by 
Sanford Gordon and Kenneth S. Drellishak 
NACA RM E57G16a, Nov. 1957, 62 pp. 

Estimation of Specific Surface of Finely 
Divided Magnesium, by Murray L. Pinns, 
NACA RM E56D10a, Aug. 1956, 28 pp. 
(Declassified from Confidential by author- 
ity of NACA Res. Abs. 121, p. 18, 11/4/57.) 

Design Considerations of a Condensing 
System for Magnesium, by 
Walter R. Witzke, George M. Prok and 
Thomas A. Keller, NACA RM E55120, 
Nov. 1955, 19 pp. (Declassified from 
Confidential by authority of NACA Res 
Abs. 121, p. 18, 11/4/57.) 

Quarterly Periodic Status Report, .\/ass. 
Inst. Tech., Hydrogen Peroxide Labs., 
Sept. 1957, 18 pp. 

Restrictive Specifications Will Increase 
Jet Fuel Costs, by H. R. Porter, G. D. 
Gould and M. L. Wolden, SAE Prep., 
Aug. 1956, 8 pp., 1 tab., 3 fig. 

Unstable Chemical Species: Free Radi- 
cals, Ions, and Excited Molecules, edited 
by Henry C. Thacher Jr., Wright Air 
Devel. Center, Tech. Rep. 57--217 (ASTIA 
AD 118181), May 1957, 219 pp. (Con- 
ference of New York Academy of Sci- 
ences, 1956). 

The Luminous Mantle of Fuel-rich 
Oxyacetylene Flames, I: Spectroscopic 
Temperature Measurements, by (G. 
Marr, Canadian J. Phys., vol. 35, no. 11, 
Nov. 1957, pp. 1265-1274. 

The Luminous Mantle of Fuel-rich 
Oxyacetylene Flames, II: Free Radical 
and Continuum Intensities, and Their 
Influences on C; Emissions, by G. V. 
Marr, Canadian J. Physics, vol. 35, no. 11, 
Nov. 1957, pp. 1275-1283. 

The Thermal Decomposition of Am- 
monium Nitrate: The Roles of Nitric 
Acid and Water, by R. D. Smith, 7'rans. 
Faraday Soc., vol. 53, pt. 10, no. 418, 
Oct. 1957, pp. 1341-1345. 

Thermal Drift of Floated Gyroscopes, 
by L. E. Goodman and A. R. Robinson, 
J. Appl. Mech., vol. 24, no. 4, Dec. 1°57. 
pp. 506-508. 

Hot Wire Anemometer for Measure- 
ments in Unsteady Flow, by O. Wehr- 
mann and R. Wille, Berlin-Charlotten! 79, 
Tech. Univ., Hermann Fottinger Inst. fur 
Stromungstechnik, TN, (AFOSR-T -56- 
26), 1956, 18 pp., 18 fig. 
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Think 


Mechanical brains for missiles must be as 
tough and tiny as possible...a design 
problem that calls for experts skilled in both 
electronic computers and miniaturization. 

Arna's computer group has shrunk a 
digital computer module until it’s the size 
shown above...a feat comparable to 
squeezing the contents of a steamer trunk 
into a cigarette package. 

Right now, in fact, through new techniques 


small 


of solid state circuitry, systematic design and 
compatibility testing, a@as#a@ is producing a 
family of airborne digital computers that are 
operational under the most severe conditions 
of vibration, temperature, noise, acceleration 
and deceleration, and nuclear radiation. 
For information on our fully transistorized, 
airborne digital computers, contact 4A@A#s4, 
Garden City, N. Y. A division of American 
Bosch Arma Corporation. 
5836 
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Builders of more large, thin wall, 


high strength solid propellant 
rocket engine cases and nozzles 
for development purposes than 


any other company in America. 


RV-A-10 SERGEANT 
NIKE-HERCULES JUPITER JR. 
JUPITER SR. POLARIS O 
POLARIS A POLARIS A-1 
AIR FORCE — RE-ENTRY -— X-17 
POLARIS — RE-ENTRY X-36 
MINUTEMAN PERSHING 


AND MANY OTHER CLASSIFIED PROJECTS 


A small experienced organization geared to handle 


your development and prototype requirements for 


static and flight tests in the shortest possible time. 


MILL STREET *® PHONE 101 
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A Multiple Range Viscometer for 
iguids, by E. F. Johnson, P. E. Parisot 
gd D. R. MacRae, Project Squid, Tech. 

PR-72-M (ASTIA AD 142429), 
. 1957, 11 pp. (available only on micro- 


rd) 


he Influence of Gradient Temperature 
ids on Thermocouple Measurements, 
N.R. Johnson, A. S. Weinstein and F. 
lerle, ASME Paper 57-HT-18 (ASM E- 
iE Heat Transfer Conference, Univer- 

Park, Pa., Aug. 11-15, 1957), Aug. 

, 20 pp. 
The Effectiveness of Stub Fins as De- 
mined by the Teledeltos Paper Analog, 
l’. Landis and T. Zupnik, ASME 
wr 57-H5-20 (ASME-AIChE Heat 
msfer Conference, University Park, 
Pa., Aug. 11-15, 1957), Aug. 1957, 12 pp. 
Report on Countermeasures, by Philip 
1 Klass, Aviation Week, vol. 67, no. 21, 

Nov. 25, 1957, pp. 102-109. 

Avionics Fights Its Own Silent War, 
wy Philip J. Klass, Aviation Week, vol. 
7, no. 20, Nov. 18, 1957, pp. 50-63. 

Hybrid Analog-digital Computing Sys- 
fms, by Max Palevsky, Instruments and 
futomation, vol. 30, no. 10, Oct. 1957, 
gp. 1877-1880. 

High-temperature Measurements, by 
i. M. Beede, Instruments and Automation, 
ml. 30, no. 10, Oct. 1957, pp. 1896-1898. 
Tube Developments for Guided Missile 
Applications, by R. W. Slinkman, Sy/- 
tania Technologist, vol. 10, no. 4, Oct., 
1957, pp. 102-105. 

Application of Statistical Theory to 
Beam-rider Guidance in the Presence of 
Noise, I1: Modified Wiener Filter The- 
wy, by Elwood C. Stewart, NACA RM 
A55Ella, Aug. 1955, 48 pp. (Declassified 
fom Confidential by authority of NACA 
Res. Abs. 121, p. 15, 11/4/57.) 

Analysis of a Flight Investigation at 
Supersonic Speeds of a Simple Homing 
System, by Robert A. Gardiner, Clarence 
L. Gillis and G. B. Graves Jr., NACA 
RM L.55J28, Jan. 1956, 55 pp. (Declassi- 
fed from Confidential by authority of 
VACA Res. Abs. 121, pp. 27, 11/4/57.) 

Theory of Aircraft Structural Models 
Subject to Aerodynamic Heating and Ex- 
ternal Loads, by W. J. O’Sullivan Jr., 
VACA TN 4115, Sept. 1957, 48 pp. 

Simulating Kinetic Heating in High- 
Speed Aircraft, Engng., vol. 184, July 26, 
1957, p. 111. 

A Comparison of the Calculated Maxi- 
mum Maneuver-response Characteristics 
of Three Air-to-air, Beam-rider Guided 
Missiles Having Different Lift Ratios, by 
Howard F. Matthews and Elwood C. 
Stewart, NACA RM A51F18. Sept. 1951, 
pp. (Declassified from Confidential 
by authority of NACA Res. Abs. 121, p. 
13, 11/4/57.) 

Development of a Piston-compres- 
sor Type Light-gas Gun for the Launching 
of Free-flight Models at High Velocity 
by A. C. Charters, B. Pat Denardo and 
Vernon J. Rossow, NACA TN 4143, Nov. 
1957, 95 pp. (supersedes NACA RM 
A55G11). 

Russian Missiles, Aviation Week, vol. 
67, no. 20, Nov. 18, 1957, pp. 26, 29, 30. 

Hypervelocity Missile Launcher, by 
B. M. Shepard, Ind. & Engng. Chem., vol. 
49, Dec. 1957, pp. 1967-1968. 

Missiles up to Date, by John W. R. 
Ta) lor, Air Power, vol. 5, Oct. 1957, pp. 
(-16. 

Upper Atmosphere to Power Air Force 
Research Vehicle (Project HARE), by 
William Beller, American Aviation, vol. 
21, Dec. 2, 1957, pp. 38-39. 
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TYPICAL ‘APPLICATIONS of 


“Pyrocore’ Ipniter. You can 


string, thread, weave it through or around the propellant. 


New DuPont Pyrocore’ Ignite 


... provides faster, more uniform ignition 


TRADEMARK 


of rocket propellants 


“‘Pyrocore” is a small diameter flexible 
metal tubing, available in practically 
any length you desire, containing a 
detonating-ignition composition as its 
core. 

Propellant ignition may be effected 
along its length at velocities ranging 
from 12,000 to 21,000 feet per second, 
depending upon the type and amount 
of core composition you specify. 

In laboratory tests, ‘““Pyrocore” cut 
ignition time 99.5% when compared 
with a standard primer. Total ignition 
of a 2114” length of cannon primer was 
achieved in less than 14 millisecond, as 
against 50 milliseconds required with 
a squib primer. 


Greater Design Freedom 


The stringlike form of ‘““Pyrocore” gives 
you new freedom in designing rocket 
motors. For example, you can thread 
“Pyrocore” through jelly roll or basket 
type igniter assemblies, or string it 
around the propellant grains. Or you 
can weave coils or loops of ‘“‘Pyrocore” 
directly into the propellant grain either 
during or after grain formation. 
Further, “Pyrocore” has very low 


DUPONT 


5. Pat. OFF 


brisance (shattering effect), and can 
convey a non-electric stimulus safely 
past expensive equipment and _ bulk- 
heads. 

Initiator assemblies may be located 
outside the propellant chamber, elimi- 
nating a major source of residue in the 
reaction zone. The ‘“Pyrocore’’ itself 
leaves very little residue. The portion 
of “Pyrocore” not actually used for 
ignition may be easily insulated with 
any one of a variety of plastic com- 
pounds. 


Sample Kit Available 


Sample kits, con- 
taining several 
variants of the 
“Pyrocore” Ignit- 
er and various 
primers and ac- 
cessories, are 
available from 
DuPont. Kit 
prices and _ tech- 
nical bulletins can 
be obtained by writing to E. I. du Pont 
de Nemours & Co. (Inc.), Explosives 
Department, Wilmington 98, Del. 


“PYROCORE” IGNITER 


TRADEMARK 


PRODUCT OF DU PONT RESEARCH 


Better Things for Better Living ... through Chemistry 
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“MONOBALL” 


Self-Aligning Bearings 


ROD END 
TYPES 


CHARACTERISTICS 


ANALYSIS RECOMMENDED USE 


Stainless Steel For types operating under high temper- 
Ball and Race (800-1200 F.) 


Chrome Alloy { For types operating under shih radial 
Steel Ball and Race ultimate loads (3000-893,000 Ibs.). 


Bronze Race and { For types operating under normal loads 
Chrome Stee! Ball with minimum friction requirements. 

Thousands in use. Backed by years of service life. Wide variety 
of Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
similar size range with externally or internally threaded shanks. 
Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
have been used to meet specific needs. Write for revised Engi- 
neering Manual describing complete line. Dept. JP-58. 


SOUTHWEST PRODUCTS CO. 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


RESEARCH 
ENGINEERS 


Explosives 


The expanding Explosives Research Raskin of 
Armour Research Foundation has immediate 
openings and opportunities for Research Engi- 
neers and Scientists interested in research con- 
cerned with properties and behavior of explosives, 
explosives components, explosives trains, fuzes 
and warheads. 


These positions offer excellent employee benefits, 
tuition free graduate study and good salaries. 


Men experienced in this field who desire to work 
for a progressive organization with some of the 
leading scientists in this field, please send resume 
to: 

E.P.Bloch 

ARMOUR RESEARCH FOUNDATION 
of 

Illinois Institute of Technology 
oa 10 West 35th Street 
Chicago 16, Illinois 
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